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1. Introduction and summary 


QCD displays an interesting phase structure as the number of flavors Nf and number 
of colors N c are varied. It is natural to discuss the phase diagram in the Veneziano 
limit [1,2]: 

N f 

N c — y oo , Nf —* oo , x = — fixed , g 2 N c hxed , (1.1) 

as a function of the variable x which has become continuous in this limit. The 
“standard” expectation for the diagram is shown in Fig. 1 (at zero temperature and 
quark mass). We restrict to the interval 0<x<ll/2 = xbz where the theory is 
asymptotically free. Various regimes can be identified: 

• The QCD regime 0 < x < x c where the infrared (IR) dynamics is similar to 
ordinary QCD (having N c = 3 and a few light quarks), with confinement and 
chiral symmetry breaking. 

• The walking regime with 0 < x < x c and x c — x -C 1 where the coupling 
constant of the theory varies very slowly, i.e., “walks”, over a large range of 
energies. 

• The conformal window x c < x < Xbz where the theory runs to an IR hxed 
point (IRFP), and chiral symmetry is intact. 

The existence of the conformal window is solid in the Banks-Zaks (BZ) limit x —> 
Xbz, because the value of the coupling is parametrically small and perturbation theory 
is trustable [3]. It is also credible that dynamics at small x < 1 is similar to ordinary 
QCD - the ~ 1 /N% corrections arising in the Veneziano limit are not expected to 
change the picture qualitatively. But the nature of the “conformal transition 1 ' at 
x = x c , and the behavior of the theories near the transition, is an open question. 

The phase diagram of Fig. 1 with a walking regime is obtained in the Dyson- 
Schwinger approach [4]. The transition is then of the BKT type [5], associated 
with the so-called Miransky scaling law [6]. The existence of the walking regime is 
important, since theories in this region may have properties, which are desirable for 
technicolor candidates [7,8] (see also the reviews [9]). The location of the transition 
has been estimated by using different approaches [10]. However, it has also been 
suggested that the transition is discontinuous, such that the dynamics “jumps” and 
walking is absent [11]. There is an ongoing effort to clarify these issues by using first- 
principles lattice simulations [12,13], but as it turns out, obtaining reliable results in 
the transition region is difficult. 

The holographic V-QCD models [14] also have the phase diagram of Fig. 1, 
including a BKT transition an walking. The “V” in V-QCD refers to the Veneziano 
limit of (1.1), where the models are defined. V-QCD is based on two building blocks, 
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Figure 1: The phases of QCD in the Veneziano limit (at zero temperature and quark 
mass) as a function of x = Nf/N c . 

the improved holographic QCD (IHQCD) [15] for the gluon sector, and a method for 
adding flavor by inserting a pair of space filling DA — DA branes [16,17]. This means 
that the gluon sector is described by a five dimensional Einstein-dilaton gravity, 
and the flavor sector is described by a tachyonic Dirac-Born-Infeld action. The 
sectors are fully backreacted in the Veneziano limit. The structure of V-QCD at 
finite temperature and chemical potential has been studied in [18-21], Two-point 
correlators and bound state masses were analyzed in [22,23]. 

There have also been numerous other holographic models addressing the phenom¬ 
ena expected in QCD at finite values of x. Walking gauge theories have been modeled 
by using as a starting point the traditional bottom-up models for example in [24-26]. 
Walking within the top-down framework has been found and studied in [27]. The 
conformal transition [28] was studied in a top-down setup in [29], and by using a 
tachyon-Dirac-Born-Infeld action in [30]. IHQCD has been used to study walking 
and IR conformal theories by modifying the holographic RG flow “by hand”, without 
inclusion of dynamical fermionic degrees of freedom [31,32], Walking dynamics and 
the conformal transition have also been studied in Dynamic AdS/QCD [33,34], which 
has partially similar ingredients as V-QCD. 

Including effects due to finite quark masses in holographic models for QCD is 
important for various reasons. First of all ordinary QCD, which describes the strong 
interactions observed in nature, has finite quark masses. Therefore, proper under¬ 
standing of the dependence of the model on the quark masses should lead to better 
models for ordinary QCD. 

Knowledge on the dependence of quark masses is also useful when analyzing 
the data from lattice simulations, which are often carried out at unphysically large 
quark masses due to technical reasons, and therefore extrapolation in the quark 
mass is needed. For lattice studies which aim at uncovering the phase diagram of 
QCD as a function of Nf, extrapolations in (typically flavor independent) quark 
mass are particularly important in the probably most interesting region near the 
conformal transition where lattice simulations are demanding. For example, the so- 
called “hyperscaling relations” [35] in the conformal window have turned out to be 
useful close to the conformal transition [36], Holography may further help to analyze 
the physics of the transition, because it can provide a unified description of the phases 
in the vicinity of x = x c where most important observables can be computed rather 
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easily for all values of the quark mass. 

Finally, it is also important to ensure that the behavior of the holographic model 
is realistic for all values of the quark mass. Studying limiting behavior (such as the 
limit of large quark mass) may lead to nontrivial constraints on the model. Fixing 
the behavior of the model to agree with QCD in limiting cases where field theory 
computations are tractable, is also expected to improve the model for all values of 
the quark mass and other parameters. 

In this article we study in detail the dependence of holographic QCD on a flavor- 
independent quark mass m q in the Veneziano limit. We argue that many of the basic 
results, such as the dependence of the energy scales on x and m q , are essentially uni¬ 
versal, i.e., independent of the details of the holographic model (given some natural 
assumptions which will be specified below). Additional scaling results for meson and 
glueball masses, the decay constants, the chiral condensate, the S-parameter, and the 
critical temperatures are derived for V-QCD analytically at large and small rn q and 
in the different regimes for x. Explicit results are computed numerically in V-QCD 
and are shown to agree with the expectations from the analytic studies. 

The numerical analysis of [14,18, 23] was mostly done at m q = 0, but some 
observables were already computed for a few values of m q and for a limited range 
of x. Here we extend these results to cover all relevant regimes (with 0 < x < x-qz) 
on the ( x , mQ-plane. Importantly, this extension does not require tuning the models 
or adding any new terms in the V-QCD action, as the value of the quark mass is 
determined through the boundary conditions of the model. Actually, we will use 
the choice for the V-QCD action with various potential terms in the action defined 
exactly as in [23]. 

We carry out a particularly detailed analysis of the dependence of the chiral 
condensate on rn q , again comparing analytic formulas with explicit numerical results 
in V-QCD. It is then shown how this analysis is used to study the deformation of 
QCD by a four-fermion operator ~ (qq) 2 - Understanding the effect of four-fermion 
operators is interesting for purely theoretical reasons, but also because such terms 
naturally arise in technicolor models due to the so-called extended technicolor inter¬ 
actions [9]. 

This article is organized as follows. In the remaining part of introduction we 
summarize our main results, and discuss the status and future of the exploration 
of V-QCD. In Sec. 2, we give a brief review of the V-QCD models. In Sec. 3, we 
discuss the universal scaling results for the energy scales of (holographic) QCD in the 
Veneziano limit, without explicitly referring to V-QCD. These results are confirmed 
numerically for V-QCD in Sec. 4, and extended to include the m^-dependence of the 
bound state masses and decay constants. In Sec. 5, we analyze the mass dependence 
of the chiral condensate, and prove the Gell-Mann-Oakes-Renner (GOR) relation for 
the pion masses for the fully backreacted case. The results of Sec. 5 are then used to 
study four-fermion deformations in Sec. 6. The mass dependence of the S-parameter 
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and the pion decay constant is analyzed in Sec. 7. Finally, a finite temperature is 
switched on in Sec. 8, where the scaling laws for critical temperatures are derived 
and compared to numerical results. The Appendix provides technical details for the 
derivation of the results. 

The article is long but many of its sections are largely independent, so a reader 
only interested in a specific topic may want to jump directly to the corresponding 
section. There are, however, exceptions: Sec. 4 requires studying Sec. 3 first, and 
Sec. 6 requires Sec. 5. 

1.1 Summary of results 

Let us then summarize the main results of this article. We identify three different 
regimes on the (x,m q ) -plane, shown schematically in Fig. 2, where the dependence 
of the model on the quark mass is qualitatively different. In regime A, the quark 
mass is a small perturbation. Regime B is the “hyperscaling” region where the cou¬ 
plings constant walks, and the amount of walking is controlled by the quark mass 1 . 
Regime C refers to the limit of large quark mass (in units of the scale of the ultravi¬ 
olet (UV) RG flow A uv ). The white regions in the plot indicate crossovers between 
the regimes A, B, and C - there are no phase transitions at finite quark mass and 
zero temperature. 

In this article, we discuss in detail how the diagram arises quite in general in 
holographic 2 models. We show how scaling results for the energy and mass scales 
in the various regimes can be obtained by using relatively simple assumptions and 
straightforward analysis (explicit results are given in Eqs. (3.9), (3.10), and (3.14)). 
The most important findings are: 

• The physics at small m qi including regimes A and B of Fig. 2, is universal, i.e., 
qualitative features are independent of the details of the holographic model. 
That is, the results are proven for all V-QCD models with such “regular” po¬ 
tentials that there is a BKT transition and therefore the model displays the 
structure of Fig. 2, but they also hold in other models which involve a BKT 
transition triggered by the same mechanism as in V-QCD (details on this mech¬ 
anism are given in Sec. 3 and in Appendix A). 

• At large m q (in regime C) the results are model dependent, and can directly 
be compared to field theory results which are tractable in the limit of large m q . 
We demonstrate that V-QCD reproduces most important features such as the 
decoupling of the massive quarks. Remarkably, we find that a specific subclass 
of the V-QCD models, which was found to be closest to QCD by analyzing the 

1 In this article the term “walking regime” always refers to the region with x c — x <C 1 where 
walking is found at zero quark mass, even though the coupling constant also walks in regime B. 

2 As such, the structure of Fig. 2 is not surprising, and one can argue how it arises from QCD 
by using arguments directly based on field theory [35,37]. 
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asymptotic meson spectra in [23], also works best at large quark mass. These 
are the models with “potentials I” below. 

We compute numerically the meson and glueball masses in V-QCD, and show that 
they agree with the generic scaling results in all regimes. 

Notice that Fig. 2 only shows changes in the dependence of the quark mass. There 
are also other features: In regime A, the theories close to x = x c are walking and 
therefore much different from the theories at low values of x, reflecting the behavior of 
the m q = 0 backgrounds, which are perturbed by the quark mass. Within regime B, 
the scaling laws change in a nonanalytic manner exactly at x = x c (see Eqs. (3.9) 
and (3.10)). 

We also analyze in detail how the chiral condensate depends on the quark mass. 
The main results are the following: 

• In the QCD and walking regimes, there is an interesting spiral structure, which 
we call the Efimov spiral. We give an asymptotic formula [38] of the spiral at 
small m q in Eq. (5.14), and compare this to data in Fig. 9. We demonstrate 
how the spiral is related to the subleading Efimov vacua 3 of the theory and to 
the Miransky scaling law in the walking regime. 

• The standard holographic proof of the GOR relation is extended to the case of 
V-QCD, which requires handling of the full backreaction and the logarithmic 
corrections which appear in near the UV boundary due to the RG flow. The 
relation is given in Eq. (5.25) and checked numerically in Fig. 11. 

We use the results for the chiral condensate to check how Witten’s method [40] 
for adding multi-trace deformations in holographic models by modifying the UV 
boundary conditions works for V-QCD. The resulting phase diagram for the case of 
a double-trace (four-fermion) deformation oc g 2 {Q(l) 2 at zero quark mass is shown in 
Fig. 13 in the (x, g 2 )-plane. For positive coupling g 2 , we see no change with respect 
to g 2 = 0, whereas negative g 2 induces an instability. 

The mass dependence of the S-parameter, the pion decay constant, and related 
quantities is studied in detail. 

• As any finite m q is turned on in the conformal window, the S-parameter discon- 
tinuously jumps from zero to a O ( NfN c ) number. Except for this discontinuity, 
the dependence on x and m q is weak. 

• We demonstrate a novel power-law scaling of the subleading terms in the S- 
parameter at small quark mass in the conformal window and in the walking 
regimes (see Figs. 19 and 20). It is argued that the power can be expressed in 
terms of the dimension of TrF 2 at the IRFP as in Eqs. (7.28) and (7.30). 

3 The terminology refers to the analogous Efimov effect in the formation of three-body bound 
states of identical bosons [39]. 
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• We analyze the ^-dependence of the S-parameter by writing it as a series over 
the contributions from the low-lying vector and axial vector meson poles (see 
Eq. (7.31)). It is shown that the increase of the S-parameter with x in the QCD 
and walking regimes can be attributed to slower convergence of the series. 

Finally we analyze the dependence of the critical temperatures of the deconfine¬ 
ment transition and various crossovers on m q and x. The scaling laws of (8.5) - (8.8) 
and (8.10) - (8.11) are demonstrated numerically in Fig. 22. The second order chiral 
transition, which is found at zero quark mass in the walking regime [18], transforms 
into a crossover as finite m q is turned on. We also discuss how the thermodynamics 
approaches that of the YM limit (x —> 0) as the quark mass is taken to infinity so 
that the quarks are decoupled. 

1.2 Outlook 

This article is part of an ongoing program [14,18-23] for studying the properties 
of the V-QCD models, and in more general the structure of QCD in the Veneziano 
limit. There are several possible future directions to explore. As a continuation of 
this study, one could consider the case of flavor dependent quark masses, which would 
be interesting in order to construct more realistic models for ordinary QCD, where all 
quark masses are unequal. In holography, this means that the background solutions 
are nontrivial dependence on the flavor indices, and consequently are described in 
terms of a non-Abclian Dirac-Born-Infeld action. The precise definition of such a 
non-Abelian action is not known, so requiring the physics to be correct might lead 
to interesting constraints for it. 

The study of the CP-odd terms of the V-QCD action is in progress at the mo¬ 
ment [41]. These terms govern the physics of the axial anomaly and the theta angle 
of QCD. Similarly to the analysis at zero values of theta, requiring the regularity of 
the solutions in the IR, and among other things the correct behavior of the asymp¬ 
totic meson spectra at finite values of theta result in constraints for the potentials 
appearing in the CP-odd terms, in addition to the constraints obtained in the YM 
limit [15,42,43]. After analyzing these constraints, we may make physically reason¬ 
able choices for these potentials, and in particular for their asymptotic behavior near 
the UV boundary and deep in the IR. 

So far the potentials of the V-QCD have not been tuned to fit any nonperturba- 
tive QCD data, and therefore all predictions of the model are qualitative. But once 
the CP-odd sector of the action has been fixed, one can start fitting the parameters 
of the potentials both to experimental results for QCD (such as meson masses) and 
to lattice data (both for QCD at finite Nf/N c [12] and for YM [44,45]). The hope is 
that the overall fit, together with the other constraints for the potentials, fixes the 
predictions of the model to a good accuracy for all relevant values of the parameters 
(such as x, the quark mass, temperature, and chemical potential). Consequently the 


model would give a effective description of QCD with real predictive power, rather 
than being just a toy model. 

2. V-QCD 

In this section, we will briefly introduce a class of bottom-up models for QCD, which 
we call V-QCD [14]. The V in the name refers to the fact that the models are defined 
in the Veneziano limit: 

N f 

N c — y oo and Nf —> oo , with x = — and g 2 N c fixed . (2.1) 

*’ C 

V-QCD is based on two “building blocks”. The first block (the glue sector of 
V-QCD) is IHQCD [15] which is a bottom-up model for Yang-Mills (YM) theory 
inspired by five-dimensional noncritical string theory. The second block (the flavor 
sector of V-QCD) is a framework for adding flavor via tachyonic Dirac-Born-Infcld 
actions [16,17]. This framework has been tested previously [46] in the probe (or 
’t Hooft) limit, i.e., without including the baekreaction of the flavor branes to the 
background. However in V-QCD, and more generally in the Veneziano limit, the 
flavor and glue sectors are fully backreacted. 

2.1 The V-QCD action 

Let us then discuss briefly the dictionary of V-QCD. The most relevant fields are 

• The dilaton 0. The exponential e~^ is dual to the operator TrT 2 . We will 
denote A = below. As this notation indicates, its background value is 
identified as the’t Hooft coupling on the field theory side. 

• The tachyon T %3 which is a Nf x Nf matrix in flavor space. The combination 
T + T' is dual to the operator q^qj whereas T — T' is dual to qc^q r For the 
background solutions considered here we will take T t3 = r(r)<V/ so the flavor 
structure does not appear explicitly. 

• The left- and right-handed gauge fields A R / R which are dual to <?7 M (1 ±75 )q. 
They are also matrices in flavor space but we have hidden the flavor indices 
here. These fields evaluate to zero for the backgrounds considered in this article. 

In addition, the scale factor A of the metric 

ds 2 = e 2A{r \-dt 2 /f(r) + dx 2 + f{r)dr 2 ) ( 2 . 2 ) 

of the vacuum solution is identified as the logarithm of the energy scale on the held 
theory side. The blackening factor / in the metric may be either identically equal to 
one or a nontrivial function of r. Our convention will be that the UV boundary lies at 
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r = 0, and the bulk coordinate therefore runs from zero to infinity (when f(r) = 1) 
or up to a horizon at a finite value of r (when f(r) has a nontrivial prohle). The 
metric will be close to the AdS metric in the UV: A ~ — log (r/£), where l is the 
(UV) AdS radius. In the UV, r is therefore identified roughly as the inverse of the 
energy scale of the dual held theory. 

The action for the V-QCD model consists of three terms: 

S = S g + Sf + S a (2.3) 

where S g , Sf, and S a are the actions for the glue, flavor and CP-odd sectors, respec¬ 
tively. Explicit expressions for the Erst two terms will be given below. As discussed 
in [14], only these terms contribute in the vacuum structure of the theory if the phases 
of the quark mass matrix and the theta angle vanish. The last term S a is important 
for the realization of the theta angle and the axial anomaly of QCD [17]. This term 
has been written down explicitly in [23], and it will be zero for all configurations 
discussed in this article. 

The glue action is that of IHQCD [15]. It includes five dimensional Einstein 
gravity and the dilaton A = eA 

Sg = M 3 N 2 c J d 5 x \J — det g (r - + V fl (A)) . (2.4) 

The flavor action is the generalized Sen’s action [17,47] (see also [48]), 

Sf = -]r,M 3 N c Tr f <fx (v)(A, T*T) y/~ det A t + V,(\ 7T*) v 7 - det A s ) , 

(2.5) 

where the quantities inside the square roots are defined as 
A lmn = g M N + w( A, T)F$ n + [{D m T)\D n T) + ( D n T)\D m T)\ , 

Armn = 9mn + w( A, T)F™ + [( J D M T)( J D JV T) t + (J D N T)(D M T) ] ] , (2.6) 

with the covariant derivative 

DmT = 8mT + iTA l m — iA^T . (2-7) 

The trace Tr is over the flavor indices - recall that the fields Ar, Ar as well as T 

are Nf x Nf matrices in the flavor space. 

It is not known, in general, how the determinants over the Lorentz indices in (2.5) 
should be defined when the arguments (2.6) contain non-Abelian matrices in flavor 
space. However, for our purposes such definition is not required: our background 
solution will be proportional to the unit matrix I/v f , as the quarks will be all massless 
or all have the same mass m q . In such a case, the fluctuations of the Lagrangian are 
unambiguous up to quadratic order. 
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The form of the tachyon potential that we will use for the derivation of the 
spectra is 

Vf(X, 7T t ) = l/ / 0 (A)e- a(A)Trt . (2.8) 

This is the string theory tachyon potential where the constants have been allowed to 
depend on the dilaton A. For the vacuum solutions (with flavor independent quark 
mass) we will take T = r(r)I Nf where r(r) is real, so that 

V>(A, T) = . (2.9) 

The coupling functions , T) and w(\,T) are allowed in general to depend on T, 
through such combinations that the expressions ( 2 . 6 ) transform covariantly under 
flavor symmetry. In this article, we will take them to be independent of T, emulating 
the known string theory results. Under these assumptions, and for the vacuum 
solutions (so that the gauge fields also vanish) the flavor action simplifies to 

Sf = - M 3 N c N f J d 5 x v 7 - det g U / 0 (A) e~ aWr '' *y/l + g rr k{\) (t') 2 . ( 2 . 10 ) 

2.2 Potentials and the holographic RG flow 

In order to fully fix the action, the potentials V g (X), Vf 0 ( A), a(A), k(A), and w( A) 
need to be specified. It turns out [14] that V g (A) must satisfy the same constraints 
as in IHQCD [15]. The other potentials will be subject to analogous constraints. We 
review the main idea here, and the details can be found in [14,23]. 

First, identification of the field A as the’t Hooft coupling and the scale factor A as 
the logarithm of the field theory energy scale defines the holographic renormalization 
group (RG) flow and the holographic beta function for the coupling as in IHQCD, 

A ' (r) 

W) = ^ (2- 11 ) 

with the understanding that the fields are evaluated on the r-dependent background 
solution. In IHQCD, the dilaton potential V g can be directly mapped to the holo¬ 
graphic beta function [15] at any value of r, that is, at any energy scale. In V-QCD, 
there is an additional field, the tachyon, whose background value is linked to the run¬ 
ning quark mass. Therefore one may define a holographic gamma function, which 
controls the holographic RG flow of the quark mass [14], The mapping between the 
beta and gamma functions and the potentials is, in general, more complicated than 
in IHQCD, but simplifies in the UV. 

The behavior of the potentials in the UV (where A —> 0) is then restricted by 
requiring that the holographic beta and gamma functions match with their QCD 
counterparts in the UV. For the UV structure to be consistent, all potentials are 
chosen to be analytic at A = 0, and the series coefficients can be related to those of 
the perturbative beta and gamma functions in QCD. It turns out that the dilaton 
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potential V g is consequently mapped to the perturbative beta function of YM theory 
as in IHQCD. Due to the backreaction, the beta function of QCD (in the Veneziano 
limit) is mapped to the combination V g — xVfo. This mapping leaves one unde¬ 
termined parameter, the UV normalization of V/o, which we call Wq. The gamma 
function of QCD fixes the UV behavior of the ratio cl/k. Here we will match the 
expansions of the potentials around A = 0 up to two loops for the beta function and 
up to one loop for the gamma function of QCD. 

Using perturbative QCD to determine the UV behavior of the potentials may 
be surprising, since holography is in general not expected to work at small values of 
the coupling A. The idea is, however, that by using this procedure correct bound¬ 
ary conditions for the more interesting IR dynamics are obtained. Notice also that 
the procedure can be seen as a rather mild generalization of what is usually done 
in bottom-up holography. For example, the bulk mass of the tachyon is typically 
required to satisfy the relation —m 2 £ 2 = A(4 — A) (at least in the UV), where £ is 
the UV AdS radius and A is either the dimension of the quark mass or the chiral 
condensate. Here this relation is effectively generalized to include loop effects, i.e., 
the perturbative anomalous dimension of the quark mass, which is roughly mapped 
to the first few coefficients in the expansion of the bulk mass of the tachyon at A = 0. 

In this article, we will carry out one more check which demonstrates that our 
UV boundary conditions make sense. Namely, we prove that the GOR relation holds 
even in the backreacted case and that the UV RG flow of the quark mass and the 
condensate, imposed by the matching to perturbative QCD, cancels in this relation, 
as it should (see Sec. 5.3). 

The choice of the potentials in the IR is more relevant since it affects how the 
nonperturbative physics of QCD is modeled. Since we are working with bottom- 
up models, there is a lot of freedom in choosing the potentials, and it is important 
to choose them such that the IR physics resembles that of QCD. The mapping to 
the beta functions is not useful at large values of the coupling because of scheme 
dependence. The asymptotic behavior of the potentials at large A can however be 
constrained heavily by comparing to several different observables, most importantly 
the asymptotics of the spectra at large mass. 

The IR behavior of the dilaton potential V g , i.e., its asymptotics as A —* oo, 
can be fixed by requiring (among other things) confinement and correct asymptotic 
behavior of the glueball spectrum at large excitation numbers. The remaining param¬ 
eters have been fitted to YM data [43]. Similarly, asymptotics of the meson spectra 
sets strict constraints on the large A asymptotics of the other potentials in the V- 
QCD action [23]. The remaining degrees of freedom will be fitted to experimental 
and lattice data in future studies. 

In this article, we will be using the choices “potentials I” and “potentials II” 
which are exactly the ones given in [22,23], and are defined explicitly in Appendix G. 
The potentials I reproduce more accurately qualitative features of QCD. This choice 
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has the UV parameter W 0 set to a constant value 3/11, and assumes w = k. The 
potentials II are included in order to study the model dependence of the results. This 
choice has the UV parameter Wq set to a value which guarantees “automatically” the 
Stefan-Boltzmann normalization of pressure at high temperatures [18], and assumes 
w — 1. Most of the results derived in this article are qualitative (e.g., scaling laws of 
various observables) and therefore insensitive to the details of the potentials. 

2.3 Background solutions 

Let us discuss some general features of the background solutions of the V-QCD 
models, first restricting to the standard case, which has a phase diagram similar to 
what is usually expected to arise in QCD. Such a phase diagram is obtained if the 
potentials are chosen as discussed above. 

2.3.1 Zero temperature 

In this article we will mostly discuss solutions at zero temperature. In this case, the 
blackening factor / in (2.2) is trivial, f = 1. To find the background, we consider 
r-dependent Ansatze for A, and A. As pointed out above, we assume that the quark 
mass is flavor independent, and therefore take T = r(r) 1^. We also set all other 
fields to zero, and look for solutions to the equations of motion (EoMs). The models 
are expected to have two classes of (zero temperature) vacuum solutions [14]: 

1 . Backgrounds with nontrivial A(r), A(r) and with zero tachyon r(r) = 0. These 
solutions have zero quark mass and intact chiral symmetry. 

2 . Backgrounds with nontrivial A(r), A(r) and r(r). These solution have broken 
chiral symmetry. As usual, the quark mass m q and the chiral condensate are 
identified as the coefficients of the normalizable and non-normalizable tachyon 
modes in the UV. 

In the first case, the EoMs can be integrated analytically into a single first order 
equation, which can easily be solved numerically. The regular solution ends on an 
IRFP, where the dilaton approaches a constant value, and the geometry is asymp¬ 
totically AdSs. In the second case, one needs to solve a set of coupled differential 
equations numerically. The regular solution ends in a “good” IR singularity [49], 
where both the dilaton and the tachyon diverge. This kind of singularity supports 
extension to finite temperature and is repulsive: perturbations around the regular 
solution develop a nonanalyticity before reaching the singularity, which signals the 
fact that IR boundary conditions are uniquely fixed. 

Let us first recall what happens at zero quark mass. The ratio x = Nf/N c is 
constrained to the range 0 <x<ll/2 = x-qz where the upper bound was normalized 
to the Banks-Zaks (BZ) value in QCD, where the leading coefficient of the /3-function 
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turns positive. The standard 4 phase diagram at zero quark mass has two phases 
separated by a phase transition at some x = x c within this range. 

• When x c < x < x-qz, chiral symmetry is intact. The dominant vacuum solu¬ 
tion is in the first class with the tachyon vanishing identically and an IRFP. 
Therefore the geometry is asymptotically AdS ,5 in the IR. 

• When 0 < x < x c , chiral symmetry is broken. The dominant vacuum therefore 
is in the second class with nonzero tachyon even though the quark mass is zero. 
The geometry ends at a good IR singularity in the IR. 

For potentials satisfying some reasonable requirements [14,23], the phase tran¬ 
sition at x = x c is due to an instability of the tachyon at the IRFP [28,50]. That is, 
the solution with vanishing tachyon is unstable if the bulk mass of the tachyon at the 
IRFP —mill, where £* is the IR AdS radius, violates the Breitenlohner-Freedman 
(BF) bound [51]: 

-m 2 Jl = A*(4-A*)<4, (2.12) 

where A* is the dimension of the quark mass at the fixed point. Notice that when 
the bound is violated, A* becomes complex. 

As a consequence of violating the BF bound, the phase transition at x = x c 
(which is only present at zero quark mass) involves BKT [5] or Miransky [6] scaling, 
for values of x right below the critical one. The order parameter for the transition, 
the chiral condensate cr ~ (gg) vanishes exponentially, 

cr ~ exp (- ) (2-13) 

V y/x c -xj 

as x —> x c from below. Here the constant K is positive. 

At finite quark mass, the BKT transition disappears: the background is always 
in the second class with finite tachyon, and the dominant vacua at all values of x are 
smoothly connected. Therefore the geometry ends in a “good” singularity in the IR. 
In particular, the IR geometry changes in a discontinuous manner (from an IRFP to 
the good singularity) when a small quark mass is turned on in the conformal window 
(x c < x < Xbz)- This discontinuity causes interesting behavior of observables which 
will be discussed in the following sections. 

The models may also have unstable subdominant vacua when 0 < x < x c . We 
will discuss such vacua in more detail in Sec. 5. 

2.3.2 Finite temperature 

At finite temperature, one can first identify two types of background geometries [18]: 

4 For some choices of potentials, also different structure can appear. In particular, there is the 
possibility that chiral symmetry breaking is absent in the regime of very small x [18]. 
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1. The thermal gas solutions which have the same (and therefore temperature 
independent) r-dependence as the zero temperature solutions, and in particular 
/ = !• 

2. The back hole solutions which have a nontrivial f(r) and end on a horizon in 
the IR. 

As usual the temperature is given by (the inverse of) the length of the compactihed 
time direction, and is equal to the Hawking temperature of the black hole for the 
second type of solutions. Both of these geometries further split into two classes, one 
having zero and the other having nonzero tachyon. Therefore there can be up to four 
qualitatively different competing saddle points. 

The finite temperature phase diagram has been studied in [18]. The following 
structure was found at zero quark mass: 

• For x c < x < Xbz, the finite temperature phase is the tachyonless black hole. 
A tachyonless thermal gas solution which has the same r-dependence as the 
zero temperature solution also exists at all temperatures but is subdominant. 
Therefore if the system is first prepared at zero temperature, any amount of 
heating makes the system jump to a different phase immediately. At tempera¬ 
tures which are much smaller than the characteristic scale of the RG flow be¬ 
tween the two fixed points, the finite temperature backgrounds are obtained by 
deforming the zero temperature backgrounds, which are asymptotically AdSs, 
only very close to the IR end. Consequently, small temperature thermody¬ 
namics is that of AdS, p oc T 4 , and the zero temperature transition is of 4th 
order. 

• For 0 < x < x c , the low temperature phase is the tachyonic thermal gas 
phase which is smoothly connected to the zero temperature solution and breaks 
chiral symmetry. The high temperature phase is the tachyonless black hole 
phase. There is always a first order (“deconfinement”) transition separating the 
phases, but it is also possible that an intermediate, chirally broken tachyonic 
black hole phase appears. If this is the case, chiral symmetry is restored at 
a separate second order transition, which has higher critical temperature than 
the deconfinement transition. As x — > x c from below, the critical temperatures 
go to zero following Miransky scaling. 

At finite quark mass, all phases are tachyonic clue to the UV boundary conditions, 
and chiral symmetry is broken. At low temperatures the dominant vacuum is the 
thermal gas phase. When the system is heated it undergoes a first order transition 
to the black hole phase, which is interpreted as the deconfinement transition (since 
there is no transition linked to chiral symmetry breaking). The phase structure is 
therefore similar for all 0 < x < Xbz- The dependence of the critical temperature 
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on m q and x will be discussed in Sec. 8. In addition to the phase transition, we find 
several crossovers which are linked to changes in the zero temperature geometry. 

3. Energy scales of (holographic) QCD in the Veneziano limit 

In this section we shall discuss the dependence of observables (at zero temperature) 
on the quark mass and x = Nf/N c for holographic models of QCD in the Veneziano 
limit and in general, i.e., not necessarily only for V-QCD. The results will not be 
proven rigorously, but we will sketch how they arise from rather natural assumptions. 
Comparison to quantitative results from V-QCD will be carried out in Sec. 4. 

It is useful to discuss separately the regions with small (possibly zero) and large 
quark mass. 

• At small quark mass, we shall demonstrate that under certain natural assump¬ 
tions, a universal picture arises from holography, which is not dependent on 
the details of the model. That is, the results of this section will apply to the 
V-QCD models with quite generic choices of potentials but also to holographic 
models with more general actions, if they satisfy some natural requirements 
which will be given below. 

• At large quark mass, the m^-dependence of energy scales and some other ob¬ 
servables (such as meson mass gaps) can be found via arguments based on 
QCD. Large quark mass probes the UV structure in holography, where predic¬ 
tions may not be reliable. We shall see that the results from holography are 
model dependent in this region, and discuss what is needed to match with the 
known results in QCD. 

In this section we will concentrate on the behavior of the various energy scales for a 
generic holographic model. In the remaining of this article we will analyze concrete 
observables such as bound state masses, the chiral condensate, and the S-parameter 
for V-QCD. Also the behavior of these observables at small quark mass is to large 
extent independent of the model details, but it is just much easier to work with the 
explicitly fixed V-QCD action. 

We will exclude the BZ limit (i.e., the limit x —> Xbz) from our analysis for the 
moment. In this region the theory is fully under perturbative control and holographic 
approach may not be useful. 

First we need to give rough definitions for various energy scales at zero temper¬ 
ature, assuming a generic holographic model. There must be a field dual to the qq 
operator, which we call the tachyon, and the geometry must be asymptotically AdS 
in the UV. We take the UV boundary to he at r = 0 as above. Because we want 
to keep the discussion generic, the definitions for the energy scales will be rather 
sketchy. Precise definitions in the case of V-QCD will be given in Sec. 4. 
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• Ayv is the scale of the UV RG flow in QCD. In holography, it can be identified 
as the scale of the UV expansions of the vacuum solution. Since V-QCD also 
implements the UV RG flow of the coupling constant through the flow of the 
dilaton field, it will be most natural to define Auv in terms of the UV expansion 
of the dilaton (see (4.1) in the next section). In the generic treatment of this 
section (which assumes small quark mass, i.e., m q f A uv -C 1), this scale is most 
conveniently defined in terms of the UV behavior of the tachyon instead: it is 
identified as the (inverse of the) boundary of the interval where the standard 
UV expression (3.1) for the tachyon holds as a good approximation. 

• Air is the soft IR scale which governs the IR expansions (or is set by an IR 
cutoff), in close analogy to the definition of Auv in the UV. Precise definition 
for V-QCD will be given in (4.3) and (4.4) in Sec. 4. 

• The quark mass m q is defined as the source for the tachyon. 

• A r is the scale of chiral symmetry breaking (whenever it is broken). In holog¬ 
raphy it can be identified as the energy scale (inverse of r) where the tachyon 
grows large (or becomes 0( 1), more precisely). As we shall see, its dependence 
on m q is similar to that of the constituent quark mass in nonrelativistic quark 
models. 

Notice that we only discuss the dynamical scales appearing through the background 
solutions, and for example the UV AdS radius is assumed to be fixed. More precisely, 
the ratios of the various energy scales are determined by the dynamics of the model, 
whereas one of the energy scales (most conveniently Auv) can be taken as a fixed 
reference scale. 

3.1 Small quark mass 

By small quark mass we mean here that m q /Auv -C 1, and m q may also be zero. 
In this case our results will be essentially universal, i.e., independent of the details 
of the holographic model. Naturally, several assumptions must be made on the 
model, in order to ensure that the physics resembles that of QCD in the Veneziano 
limit. Most importantly, there must be an IRFP and a conformal window, and 
the conformal transition at x = x c is assumed to be a BKT transition, which is 
naturally implemented in holography by the tachyon hitting its BF bound at the 
IRFP [28, 50] (implemented through a geometry which is asymptotically AdS in 
the IR). These assumptions are met by the V-QCD models with potentials which 
fulfill reasonable constraints in the UV and in the IR, and do not have any peculiar 
structure at intermediate energy scales [14,23]. But these assumptions can also hold 
more generally in holographic theories the actions of which cannot be written in the 
V-QCD form. Indeed similar results which we will present here have been found in 
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related top-down [29] bottom-up [30,33,34] models. For simplicity we also assume 
that the upper edge of the conformal window lies at the QCD value, x — 11/2 = xbz- 
The model assumptions can be made explicit in terms of the tachyon background 
in various ranges of the bulk coordinate r. We will only need the solutions for small 
tachyon, so that the tachyon EoM is approximately linear (see Appendix A for more 
detailed derivation of these solutions). In the UV (where r —> 0), we assume the 
standard behavior determined by the dimension of the quark mass: 

r ~ m q r + or 3 , (r<C-—) > (3-1) 

V Au v/ 

where a is proportional to the chiral condensate. Here logarithmic running of the 
quark mass (and the condensate) could also be included, but it is not important 
because it will not affect the leading scaling behavior. In the vicinity of the fixed 
point when the tachyon is small and the BF bound is satisfied, corresponding to 
x c < x < xbz in QCD, we find that 

T ~ C m (rA uv ) A * + C' cr (rA uv ) 4_A * , < r < , (3.2) 

where the precise x-dependence of the anomalous dimension of the quark mass A* 
depends on the model. When the BF bound is violated, and x c — x is small 

t -C w (rA uv ) 2 sin [u log(rA uv ) + (f>] , ( < r < ) , (3.3) 

\ A UV A t/ 


where v = ImA* ~ n^/(x c — x)/K as x —> x c from below. The coefficient K depends 
on the model, and it satisfies 

K = l = ■ (3,4) 

v / A[A.(A.- 4 )] I , Ie 

Now it is straightforward to fix the integration constants (a, C m , C a , C w and (j)) 
and compute the mass dependence of the various energy scales by using the following 
recipes: 


• Both the normalizable and nonnormalizable terms of the tachyon solution 
are separately continuous 5 . Therefore an approximate tachyon solution for 
all r <C 1 /A t can be found by gluing together the results from the different 
regimes (3.1), (3.2), and (3.3). Because the bulk mass of the tachyon varies 
smoothly, the terms of the tachyon can only jump by a factor O (1) as we move 
from one regime to another. 

5 Near the IRFP in (3.3) the normalizable and nonnormalizable terms cannot be separated, but 
the strategy is then interpreted as requiring the continuity of the tachyon and its first derivative. 
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• The IR behavior of the tachyon in the regime r 1/A r is nontrivial and qual¬ 
itatively different from the formulas given above. Therefore, the IR boundary 
conditions obtained by requiring continuity at r ~ 1/A r are taken to be generic, 
i.e., it is assumed that the conditions are not fine tuned to pick any specific 
solution for r <C 1/A T . 

3.1.1 Zero quark mass 

Let us first recall how the energy scales depend on x at zero quark mass. The results 
for x < x c are obtained by using the above assumptions and recipes. Therefore they 
are also independent of the details of the model. The computation can be found in 
Sec. 10 of [14] and the scaling results will also be reproduced by the analysis of Sec. 5 
below, where the chiral condensate is studied in detail (see also [29]). We will only 
list the results here. 

When x < x c , A IR ~ A r (i.e., the tachyon becomes sizeable where the asymptotic 
IR geometry starts), and A r is not defined in the conformal window (x c < x < x R z)- 
Depending on the value of x = Nf/N c , there are three regions with qualitatively 
different behavior: 


• In the QCD regime, meaning that 0 < x < x c and x c — x > 1. there is only one 
scale, Auv ~ A IR , which corresponds to Aq CD . 


• In the walking regime, which is found when x < x c and x c — x -C 1, the IR and 
UV scales are related through Miransky scaling: 


Auv 

- rs_/ 

Air 

The chiral condensate satisfies 

° ^ Air 

A uv ~ A^ 


exp 


K 

y/x c ~ X 


~ exp 


2 K 

V x c - X 


(3.5) 


(3.6) 


• In the conformal window (x > x c ) the tachyon is zero as chiral symmetry is 
intact. In this case one expects that the scales of the UV and the IR expansions 
are the same, Auv ~ Air (as there is no obvious mechanism which would 
separate the scales). 


3.1.2 Small but finite quark mass 

Let us then generalize these results to nonzero quark mass, first assuming a small 
mass (m q /Axj\/ <C 1). The detailed analysis uses the strategy formulated above, and 
can be found in Appendix A. We will only discuss the results here. 

One can identify, for m q /Auv •C 1, two regimes on the (x,m q )~ plane, shown 
schematically in Fig. 2: 


19 













m q /A uv 



0 1 2 3 4 5 


Figure 2: A schematic diagram of the different scaling regions as functions of x and m q . 


A The regime where the quark mass is a small perturbation. This is possible 
for 0 < x < x c as the m q = 0 solution is continuously connected to the 
solution with finite rn q . In order to determine the extent of regime A, we 
will perturb the solution at vanishing quark mass by adding a small m q . It 
is a small perturbation so long as the nonnormalizable term remains small for 
r <C 1/Auv- This is equivalent to m q <C | a(m q = 0)/A^j V |, or 


m q 

Auv 


< 1 


(3.7) 


in the QCD regime, and 


m q 

< exp 

Auv 


2 K 

\Jx c -x 


(3.8) 


in the walking regime. Apart from the perturbation caused by the quark mass, 
the energy scales behave as in the m q = 0 case discussed above. In particular, 
the Miransky scaling law (3.5) applies for i c - i C 1. 


B The “scaling” regime which involves walking of the coupling, and the amount 
of walking is determined by the quark mass. Continuity of the tachyon implies 
for x <x c that (see Appendix A for details) 


m g _ 

AUV Ayy A^jy 


fx < x c and 


exp 


2 K 

\Jx c - x 




m g 

Auv 
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whereas in the conformal window 


m q 

f a ir A 

Auv ~ 

V A UV J 

<T 

f AiR_y 4 ~ A: 


^Auv J 



< x < xbz and 


m q 

vW 



(3.10) 


Notice that (3.9) is obtained from (3.10) by setting A* = 2, which is indeed 
the value of the anomalous dimension as x —» x c from above, and the scal¬ 
ing behavior is therefore continuous and both expressions are valid at x — x c . 
Moreover, (3.10) gives what is termed the “hyperscaling” relation for the chi¬ 
ral condensate (see, e.g., [35]) and often written in terms of the anomalous 
dimensions 7 * = A* — 1 , 


/ \ 3 ~ 7 * 
a ( m q \ 1 + t * 

V : y Auv/ 


(3.11) 


Such scalings have been studied recently in a specific holographic model [34]. 
The behavior of A T is the same as for m q = 0 and in the regime A, i.e., 


A T ~ Ai R . 


(3.12) 


Notice that we have not discussed the m 9 -dependence of bound state masses and 
decay constants, because they are difficult to analyze without specifying the details 
of the model. In the regime A and at small x, however, it is clear that the lowest 
bound state masses (except for the light pions) and decay constants must be O (Air) 
because this is the only available scale (apart from the small perturbation due to 
the quark mass). The pions are expected to obey the GOR relation since this arises 
in holography quite in general. In the regime B, as well as in the regime A when 
x c — x -C 1 , the most natural expectation is that the masses and decay constants are 
still given by the soft IR scale O (Air). This will be demonstrated for V-QCD below. 

3.2 Large quark mass 

Let us then discuss the limit m q /Auv 1 (regime C in Fig. 2). Since the quark mass 
is large, the tachyon grows large very close to the UV boundary where the dilaton is 
still small. Therefore regime C probes the limit of small dilaton and large tachyon in 
the holographic model. Due to the smallness of the dilaton it is not obvious that the 
holographic description is reliable for all observables in this regime. But in analogy 
to how the UV structure of V-QCD is fixed in order to guarantee correct boundary 
conditions for the IR dynamics (as explained in Sec. 2), it is important that the 
holographic model is as close to QCD as possible also at large quark mass, in order 
to have the best possible boundary conditions for the physics at small and 0(1) 
quark masses. We will therefore analyze what are the possibilities in this limit. 
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Deep in the UV, for r <C 1 /m q , the tachyon will have the standard form of (3.1) 
which already implies scaling results for a and A T , with similar assumptions as above. 
At r ~ 1 /m q it will grow large, after which nonlinear effects are important and a 
generic solution cannot be written down. In particular, the RG flow is not expected 
to approach the IRFP at any point, and the IR structure of the tachyon solution 
shown above is not relevant. 

Instead, additional scaling results can be derived assuming that the holographic 
model implements some key features of QCD (as will be the case in V-QCD). These 
are the RG flow of the’t Hooft coupling in QCD, which also gives the proper definition 
of the UV scale Auy in regime C, and the decoupling of the massive quarks at energy 
scales much smaller than m q . Such decoupling is automatically implemented in V- 
QCD by the Sen-like exponential tachyon potential of the flavor action, as we will 
explain in Sec. 4. Due to the decoupling, the RG flow is that of full QCD for 
r < l/m q and that of YM theory for r 1 jm q . Continuity of the ’t Hooft coupling 
(rather than the tachyon) is required. 

The collected results for the regime C in Fig. 2 are the following (details can 
again be found in Appendix A): 

C The regime of large quark mass (m q /Auv 3> 1). The form of the tachyon 
solution implies that 


er ~ m q and A r ~ m q . 


(3.13) 


Here cr will be connected to the properly renormalized chiral condensate (see 
Appendix C for details). The renormalization also involves scheme dependence 
which is important at large quark mass (see [46] for a discussion in the context 
of holography). The condensate is proportional to m q as in (3.13) for generic 
schemes. Further, analysis of the RG flow leads to 


A 


uv 


A 


IR 


m n 


A 


uv 


bo/b™-l 


m n 


A 


IR 


( m <l \ 

Uuv; 


bo/bF* 


where 6 0 (6 q" m ) is the leading coefficient of the beta function for 
Veneziano limit (YM theory at large N c ). That is, 


b 0 _ , 2x 
b™ ~ 11 


(3.14) 
QCD in the 

(3.15) 


Let us then briefly comment on the size of the bound state masses in regime C. 
First, as the quarks are decoupled for energy scales smaller than the quark mass, 
the glueballs are expected to decouple from the mesons, and have a mass gap of 
O (Ajr). Recall that the meson states in QCD become nonrelativistic at large m q , 
and therefore their mass gap is ~ 2 m q , and the mass splitting of the low-lying states 
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is much smaller than the gap. The example of V-QCD, which we shall discuss below, 
shows that obtaining such a mass gap in holography is nontrivial (see the analysis 
for V-QCD in Appendix B), and the gap can be either O ( m q ) or O (Air) for actions 
which produce reasonable results at small m q . If an IR cutoff is placed at the point 
where the tachyon grows large (as in the dynamic AdS/QCD models [33]), mass gap 
oc m q is obtained. In the presence of such a cutoff, A T ~ m q is the only scale in 
the system. But in this case the splitting between the bound state masses is also 
expected to be O ( m q ). 

Finally let us comment on the structure of Fig. 2 near the BZ region. Namely, 
the boundary of the regime C bends toward higher m q j Auv as x grows. Similarly 
the upper boundary of the regime B bends down in the BZ region. This is a generic 
feature due to RG flow which becomes slower and slower in the BZ limit x —* xbz- 
Due to slowness of the flow the separation of the energy scales m q and Auv needs 
to be larger for the scaling results to apply. This phenomenon is discussed in more 
detail in Appendix A. 

Actually the BZ regime could be identified as an additional scaling regime. Here 
this is not done, however, since from the point of view of holography this regime is not 
the most interesting one. The coupling constant is restricted from above by its small 
value at the IRFP, the theory is perturbatively soluble, and therefore holographic 
description is not expected to be useful. 

4. Scaling of bound state masses in V-QCD 

While we discussed above the generic behavior of (holographic) QCD in the Veneziano 
limit, we shall now derive explicit predictions for the V-QCD models defined in Sec. 2, 
and demonstrate that they agree with the results of the previous section. 

4.1 Energy scales in V-QCD 

The various energy scales can be defined explicitly in terms of the background solu¬ 
tions. For definiteness we will write down the definitions here. We will use the UV 
and IR expansions which can be found in Appendix D of [14] and in Appendix D 
of [23]. The UV expansions can also be found in Appendix C.3. 

• First, Auv is the scale of the UV expansions. The precise definition is most 
conveniently given in terms of the UV expansion of the dilaton A: 

V _ 1 861 log [- log(rAuv)] / 1 A 

6 0 log(rA uv ) 96 q log(r A uv ) 2 V lo g( rA uv) 3 / 

Recall that the coefficients of the potentials in the V-QCD action were matched 
to those of the QCD beta function in the UV. We used this mapping to write 
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the coefficients in the expansion (4.1) in terms of the coefficients fe* of the QCD 
beta function in the Veneziano limit, 


/3(A) 


dX 

d\ og/i 


—6qA 2 + b\X 3 + • • • 


(4.2) 


• The IR scale A TR can be defined analogously in term of the IR expansions. 
For the standard geometry in V-QCD with an IR singularity the definition of 
Air = 1/R is given in limit r —)■ oo as [14] 

l0g A = \lp + ° ^ = \ A ^ rT ’ 2 + ° ( r °) ' ( 4 ‘ 3 ) 

In the presence of an IRFP the definition is modified to 

A ~ A* - = A* - (rAnt) -5 (4.4) 

where 6 = App — 4 is the anomalous dimension of the TrF 2 operator at the 
fixed point. 


• The quark mass is determined 6 in terms of the UV asymptotics of the tachyon 
(r -+ 0), 


r(r) 


= m q r{- log(rAuv)) p 


1 + 0 


log(rAuy) 


(4.5) 


where £ is the UV AdS radius and p = 7o/A) is the ratio of the leading coeffi¬ 
cients /?o and 70 of the beta function and the anomalous dimension of the quark 
mass in QCD, respectively. Notice that in V-QCD, the logarithmic running of 
the quark mass is therefore included, and is matched to agree with that of 
QCD. 


• The scale A r where the tachyon grows large is defined simply by 7 


r 


= 1 


.A—log A r 


(4.6) 


where r and A are to be evaluated on the background solution. 

6 It is well known that the quark mass can be defined only up to a constant. This constant is set 
to one for notational simplicity. 

7 For potentials II slightly modified definition is used: the number on the right hand side is set 
to 1/10. This is necessary because the transition to the IR region where nonlinear terms in the 
tachyon are important turns out to happen when the tachyon is still smaller than one, and in the 
nonlinear region the tachyon grows relatively slow. Therefore using exactly (4.6) would lead to an 
energy scale which does not precisely reflect the change in the dynamics. 
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The background EoMs for V-QCD are invariant under the transformation 

r —> Ar , A —>■ A — log A . (4.7) 

This transformation changes all energy scales defined above by the same number, and 
reflects the choice of the units of energy on the field theory side. Consequently, only 
the ratios of the above energy scales are a priori well-defined. Notice also that Air 
and A r do not have direct counterparts in field theory, whereas A uv is indeed mapped 
to the scale of the RG flow in held theory and m q is identified (up to a constant) as 
the quark mass on the held theory side. Therefore the values of A uv and rn q (up to 
the ambiguities mentioned above) can be inferred from QCD data only, whereas the 
determination of Air and A r also requires pinning down the holographic action. 

There is a small issue with the above definitions which will be visible in the 
explicit results below. Namely, the definition of Auv which is natural at generic 
values of x is not optimal at high x and in particular in the BZ limit. In this limit 
one would expect Auv to match the scale of the UV RG how, but this turns out 
not be be the case. Also at zero quark mass there should be only a single scale 
and therefore Auv should equal Air, but as it turns out, actually Auv becomes 
exponentially suppressed with respect to Air. Because the RG how of the coupling 
is controlled by the two-loop beta function in the BZ limit, a single scale A can be 
defined which has better behavior in this limit (see Appendix B of [52]): 

(yx “ 0 exp (yx)" ( Xr ) “ (48) 

This scale is related to Auv hi the BZ limit by 

b l 

A — ° A uv ; (4-9) 

where bi/b^ ~ (xrz — x)~ 2 grows large in the BZ limit (see Appendix A for some 
more details). Notice that Auv will anyhow be used as a reference scale in the nu¬ 
merical analysis below - since the BZ region is not very interesting from holographic 
viewpoint, we have chosen to use the same definition of Auv as earlier literature, 
even though it has unnatural behavior in this region. 

4.2 Flavor nonsinglet masses and decay constants: scaling results 

Before going to the numerical results let us argue how the scaling laws for the masses 
and decay constants can be derived in V-QCD. The meson masses in each sector (for 
vectors, axial vectors, scalars, and pseudoscalars) may be defined into two classes: 
the flavor singlet and nonsinglet states. The former are singlets under the vectorial 
SU(Nf ) transformation, whereas the latter are the other fluctuation modes, which 
correspond to quark bilinear operators involving the Hermitean traceless generators t a 
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of SU(Nf). We restrict analytic considerations to the nonsinglet states, because the 
singlet mesons mix with the glueball states, which would lead to complications [23, 
53]. Only the main points are summarized here and the details can be found in 
Appendix B. 

The fluctuation equations for the scalar nonsinglet mesons can transformed in 
the Schrodinger form as detailed in Appendices A and B of [23]. The masses for 
each sector are then given as eigenvalues of a Schrodinger equation with a certain 
potential term and the Schrodinger coordinate running from u = 0 (UV) to u = oo 
(IR). In order to find the behavior of the mass gaps and splittings one then needs to 
study the Schrodinger potentials Vs(u). 

Let us first take finite but small quark mass (m q /Auv -C 1), so that the regimes A 
and B are covered. We use here the V-QCD action, but results in this region are 
not sensitive to the details of the action. The Schrodinger potential for u <C 1/Air, 
is given by the background with small or walking dilaton and small tachyon, so that 
the geometry is close to AdS, and consequently Vs(u) ~ const /u 2 . For u 1/Air 
the diverging tachyon creates the confining potential ( V s {u ) ~ Aj R w 2 if the excitation 
spectrum is linear, m 2 ~ n with n being the excitation number). Therefore the 
only relevant scale is Air and the (lowest lying) meson masses as well as the mass 
splittings between the modes are given by this scale: 


(4.10) 


m n ~ Air , 


which is also consistent with the results from Dyson-Schwinger and Bethe-Salpeter 
approaches in regime A [54], In units of Auv, by using the scaling laws from Sec. 3 we 
therefore obtain the same results as at vanishing quark mass [22,23] within regime A: 
m n ~ Auv hi the QCD regime and the masses go to zero obeying Miransky scaling, 
m n /Auv ~ exp(— K/y/x c — x ), in the walking regime. In regime B, we find the 
“hyperscaling relations” [34,35] for the low-lying meson masses: 



(4.11) 


The sole exception to these scaling results is the pion mode (for x < x c ), which 
is massless at m q = 0 and obeys the GOR relation in regime A: 



(4.12) 


so that it is lighter than the other meson states. The GOR relation will be discussed 
in more detail in Sec. 5. At the level of the Schrodinger formalism the absence of the 
mass gap in the pseudoscalar sector is reflected in the negativity of the Schrodinger 














potential in the UV region. In regime B, the pseudoscalar masses will also obey the 
scaling (4.10). 

The various decay constants are slightly more difficult to analyze. In Appendix B 
it is argued that they are similarly of the order of Air when the quark mass is small. 

At large quark mass (regime C) the meson masses and decay constants depends 
more on the details of the action. Actually only the form of the flavor action Sf in 
the limit of large tachyon and small dilaton is relevant, as is shown in Appendix B. In 
this limit the functions a, k, and V/o become almost constants and the action of (2.5) 
takes the form of the standard DBI action with the exponential Sen potential. 

The exponential dependence Vf oc exp(—a(A)r 2 ) of the flavor potential on the 
(squared) tachyon naturally implements the decoupling of the massive quarks. For 
large quark mass the tachyon is roughly proportional to m q r in the UV (for r <C 
l/m q ). When r grows larger than 1 /m q , that is, at energies lower than m q on the 
field theory side, the tachyon grows sizeable (see Appendix B.2 for details) and the 
exponential factor in Vf decays rapidly. Consequently, the flavor part of the V-QCD 
action becomes suppressed with respect to the glue part, and therefore the dynamics 
at energies below m q is governed by the gluons as expected. The decoupling of flavors 
will also be well visible in the numerical results for the scalar singlet states below 
which involve both glueball and qq components. 

Interestingly, as shown in Appendix B, the choice Vy(A,r) = V/ 0 (A) exp(— a 0 r 2 ) 
of potentials I, where the function a(A) is set to a constant value a 0 , is a special 
case. This choice was motivated by the asymptotics of the meson trajectories at 
zero quark mass [14,23], and noticeably a(A) is also constant in tachyon potentials 
obtained from string theory [16,17,55]. In Appendix B we show that this choice is 
essentially the only one which produces physically reasonable mass gap and splitting, 
if the value of a 0 is slightly modified from the value of potentials I which reproduces 
the correct dimension of the quark mass and condensate in the UV. 

We will anyhow discuss the results for potentials I without this modification of 
do, because such a modification was not introduced for our numerical studies. It is 
found that the mass gap for the (flavor nonsinglet) mesons in all sectors (vectors, 
axials, scalars and pseudoscalars) is given by 


m 


gap 


ml 


3q=o 

4£ 2 



(4.13) 


where i is the UV AdS radius, and £ = 1 would be required to match with the 
field theory result for nonrelativistic bound states 8 . We stress that this formula was 
obtained by fixing the value of do to produce the UV dimension of the quark mass 
and the condensate, as was done for the potentials I used in the numerics. The last 
expression in (4.13) suggests that £ = 1 could be obtained by tuning the value of Wo. 

8 One can have both £ = 1 and the correct UV dimension for the quark mass if extra terms (e.g. 
a (1 + #t 2 )) are added in the tachyon potential. 
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This is, however, problematic because the contribution involving Wo should vanish 
in the probe limit x — > 0, and also because W 0 would negative, which causes Vf to 
have a node [14] at which the tachyon background equation becomes singular. 

The mass splitting is suppressed with respect to the gap as is the case in real 
quarkonia. For potentials I it scales as the inverse of the mass gap. The decay 
constants of the lowest meson states are exponentially suppressed in the quark mass 

fn ~ exp(-#mf) (4.14) 

and they are therefore decoupled. This is found for all mesons with masses below 
O ( m q ). The decay constants of the states with masses ~ m q are of the order of Air. 

4.3 Numerical results 

Computing the energy scales and masses numerically is straightforward (but tedious) 
after the background has been constructed numerically [14, 23]. The potentials I 
(given explicitly in Appendix G) are used here unless stated otherwise. We choose 
three reference values x = 1, 4, and 4.5, which lie in the QCD regime, walking 
regime, and conformal window, respectively, and plot the observables as functions of 
the quark mass. We also show plots where m q / Arv is fixed to 10 -6 and x is varied 
over the whole parameter space. These choices cover the most interesting structures 
of Fig. 2. 

Notice that when plotting dimensionful parameters as a function of x we are 
comparing different theories and a choice for the reference scale must be made. A 
natural choice would be A uv , which roughly corresponds to the scale where the ’t 
Hooft coupling takes some fixed tiny value very close to the UV boundary. In many of 
the plots below, however, A IR is chosen as the reference scale instead, simply because 
this makes the plots more easily readable. 

4.3.1 Energy scales 

In Fig. 3 the dependence of the energy scales on x and m q is demonstrated. We have 
chosen to show the scales Arv and A r in the units of Air since this makes the details 
visible. Some of the bound state masses are also shown as thin lines. 

The top-left plot shows the dependence of the scales on a; at a tiny quark mass 
(m q / A uv = 10 -6 ). The data extend only up to x = 5.1 because in the BZ region it 
is difficult to do reliable numerics. The top-right plot shows the mass dependence in 
the running regime (x = 1). The bottom-left plot is in the walking regime (x = 4, 
which is close to x c ~ 4.0830). The bottom-right plot is in the conformal window 
(x = 4.5). 

The thick solid blue curves show the ratio Arv/Air in each plot. In the top-left 
plot the crossover from the QCD-like regime to conformal window is clearly visible: 
As x — y x c from below, the ratio first grows according to the Miransky scaling 



Figure 3: Dependence of the energy scales on x and m q . Also the masses of the rho meson 
m p , the lowest singlet scalar m ss , and the pions m n are shown for reference. The blue solid 
curves is Auv, the red dashed curve is A r , the thin dotted magenta curve is the rho mass, 
the thin dotdashed green curve is the pions mass, and the thin long-dashed brown curve is 
the mass of the lowest singlet scalar in each of the plots. See text for more details. 


law (3.5) in regime A until the condition (3.8) no longer holds. Thereafter the ratio 
saturates to roughly a/A uv / m q in regime B as predicted by (3.9), and then one 
moves out of regime B at even higher x. 

The dependence of Auv/A ir on also follows the predicted scaling laws. In 
the regime A (low m q /Auv hi top-right and bottom-left plots) the ratio is constant 
as the quark mass is a small perturbation. In the regime C, i.e., at large m q / Auv, 
the ratio follows the power law of (3.14) as best seen in the top-right plot. In the 
regime B the ratio obeys the other power law of (3.9) and (3.10) as best seen at 
intermediate m q in the bottom-left plot. Notice that much of the solid blue curve 
was left out in the bottom row plots in order to make the details of the other curves 
better visible. 

The thick dashed red curves show the ratio A t /Air. At small quark mass, includ¬ 
ing regimes A and B, the ratio is close to one as expected (except in the BZ limit). 
For the explanation of the divergence of A t /Air in the BZ limit see Appendix A, 
Eq. (A. 15). At large quark mass (regime C), we End A r ~ m q as predicted in (3.13). 
At high x (plots in the bottom row) the convergence toward this scaling law is quite 
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Figure 4: Dependence of the mass ratios between the lowest excitations of each tower 
on x and m q . The choices of x and m q /Auv are as in Fig. 3. Blue solid curve is the 
lowest vector (p meson) mass, red dashed curve is the nonsinglet axial vector mass, dotted 
magenta curve is the lowest nonsinglet scalar mass, and the dotdashed green curve is the 
nonsinglet pseudoscalar (pion) mass in each plot. In addition, the long-dashed brown curve 
is the lowest singlet scalar mass. 

slow due to the slow running of the quark mass. It could be demonstrated by con¬ 
tinuing the plots up to much larger m q / A uv , but we have chosen not to do so in 
order to show the other details in the plots more clearly. Actually all variables vary 
slower and slower as functions of m q /Auv when x is grows, and therefore we have 
substantially increased the range of m q /Auv hi the plots with higher x, but this is 
still not enough to demonstrate the large m q scaling convincingly. The change in 
the m^-dependence is due to the RG flow (see the end of Appendix A for some more 
details). 

The dependence of a on m q will be discussed in Sec. 5. 

4.3.2 Flavor nonsinglet masses and decay constants 

The flavor nonsinglet spectra can be computed as explained in [23]. The results for 
different values of x and m q are shown in Figs. 3, 4, and 5. 

First, the p and n masses are shown in Fig. 3 with thin dotted magenta and thin 
dotdashed green curves, respectively. As expected, the p mass is O (Air) at small 
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Figure 5: Masses of the lowest four vector states as a function of m q at x = 1. Left: 
masses in units of Air. Right: masses normalized to the mass of the lowest vector meson 
(the p- meson). 


quark mass (regimes A and B) and obeys the power law (4.13) in regime C. The pion 
mass is close to the p mass in regimes B and C, but in regime A it obeys the GOR 
relation instead, as best visible from the top-right and bottom-left plots at small m q . 
In the top-left plot, m n / A IR is O ^^/m 9 /A uv j at small x but then increases with x 
in the walking regime (actually obeying the Miransky scaling law) until the ratio is 
O (1) at the point of the crossover near x = x c . 

Fig. 4 shows the masses of the lowest meson states (i.e., the mass gaps) in each 
sector normalized to the p mass. The choices of x and m q are the same as in Fig. 3. 
The masses of the lowest vector, axial, scalar, and pseudoscalar states are given by 
the solid blue, dashed red, dotted magenta, and dotdashed green curves, respectively 
(the brown curves give the scalar singlet mass gap which will be discussed below). 
These ratios are mostly constant and close to one as predicted by the above scaling 
arguments. The exception is the pion mass (dotdashed green curves) which obeys the 
GOR relation in regime A. Notice that in regime C all meson mass gaps should ap¬ 
proach the same number (roughly 2 m q ) as expected for nonrclativistic bound states, 
but we find instead that the axial and pseudoscalar gaps are larger than those of the 
vectors and scalars. The reasons for this are analyzed in Appendix B. Notice also 
that the lowest scalar states are lighter than vectors even at small values of x , which 
seems to be in conflict with QCD. Such details are, however, sensitive to the choice 
of the potentials in the V-QCD action, and can be changed by tuning the potentials. 

Finally we plot the masses of the four lowest vector states in Fig. 5 in order to 
demonstrate the dependence of the mass splittings in the spectra on m q . The masses 
are given in units of Air (left hand plot) and normalized to the lowest mass, i.e., the 
p mass (right hand plot). The splittings decrease with increasing m q in regime C 
which is in qualitative agreement with the bound states becoming nonrelativistic, 
but the power laws are not exactly correct (see Appendix B for details). 
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Figure 6: The masses of the five lowest scalar singlet states as functions of m q for various 
choices of x. 

4.3.3 Scalar singlet masses 

The singlet sector is qualitatively different from the nonsinglet sector because it also 
contains glueball states which mix nontrivially with the singlet meson states. Such 
mixing takes place in the scalar and pseudoscalar sectors. The scalar sector will be 
discussed in detail here while the pseudoscalar sector will be analyzed in a future 
publication [41]. 

Before going to the numerical results, let us discuss the generic features of the 
spectrum. The singlet mesons masses are expected to show similar 9 m q and x de¬ 
pendence as the nonsinglet mesons above, given by (4.10) and (4.13). In the singlet 
case there is, however, nontrivial mixing of the mesons with the glueball states, the 
masses of which should be independent of m q and therefore always characterized by 
Air. In particular in regime C the masses of the mesons become much larger than 
the glueballs which suggest that the meson and glueball states decouple. 

Let us then demonstrate these features numerically for the scalar singlet sector 
in V-QCD. The mass gap of the scalar singlets is shown by the long-dashed brown 
curves in Figs. 3 and 4. The mass gap is that of the glueballs and therefore O (A IR ) 

9 There is an exception in the singlet pseudoscalar sector: the lowest state, the rf meson, is 
anomalously light at small x in the Veneziano limit [2,56]. This is also reproduced by V-QCD [41]. 
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for all values of m q and x as seen from Fig. 3. The fact that the lowest glueball 
mass is suppressed with respect to the rho mass m p in regime C can also be seen in 
Fig. 4: the ratios m ss /m p given by the brown curves, where m ss is the mass gap for 
the scalar singlet states, decrease with m q at large m q . 

The decoupling of the meson and glueball states is most clearly demonstrated 
by Fig. 6, where the masses of the five lowest scalar singlet states are plotted as a 
function of m q in the QCD regime (x = 1, top-left), in the walking regime (x = 4, 
top-right) and in the conformal window (x = 4.5, bottom). When m q j Auv -C 1 the 
spectrum has both glueball and meson states which are nontrivially mixed. As m q 
increases, the meson masses increase while the glueball masses stay constant, which 
leads to the crossing structure seen in the plots. This would be expected to happen 
at m q /Auv ~ 1, but as was discussed above, A uv is not exactly the scale of the UV 
RG flow when x is large, and therefore the crossing structure shifts to higher values 
of m q as x increases. At large m g /Auv only the glueballs are left. Their decoupling 
from the mesons is demonstrated by the fact that the limiting values of their masses 
as m q —> oo are independent of x, and in fact it can be checked that they match with 
the glueball masses obtained in the YM limit (x —* 0) of V-QCD. 

It was shown in [22,23] that V-QCD does not have a light “technidilaton” mode [8] 
(which would be the lightest scalar singlet state) as x —> x c . Both the singlet and 
nonsinglet scalars fluctuations do have critical [28] behavior in the near conformal 
region (see Appendix I in [23]) and the Schrodinger potential for the nonsinglet 
scalars is negative (as was also found in [32]), but it was shown numerically that 
this is not enough for a technidilaton to appear. The negative result is also seen in 
Figs. 4 and 6: the lowest singlet scalar does not become light with respect to the 
other states in the walking regime (x = 4). We do see, however, that it is lighter in 
regime A than in regime B. 

4.3.4 Masses near the conformal transition 

To conclude this section, let us add a few comments on the scaling of the bound state 
masses near the conformal transition. We plot the masses of the rho meson, the pion, 
and the lowest scalar singlet state in Fig. 7. The left hand plot shows the masses in 
units of Air, and the right hand plot shows the masses in units of Auv- Notice that in 
units of Air the pion mass deviates from the masses of the other states at extremely 
small m q (regime A) where it obeys the GOR relation, but in units of Auv the pion 
mass obeys the same power law ~ \/m q /Auv hi regimes A and B. 

Knowledge of the dependence of the meson masses on m q near the walking regime 
is important for the lattice studies which take place at finite quark mass, and aim to 
locate the conformal transition at m q = 0 [13,36]. As we have pointed out above, the 
regime B extends even to x < x c (see Fig. 2), and the crossover between regimes A 
and B moves to lower x as m q is increased. This suggests that studies at finite m q 
lead to an underestimate for x c . One should recall, however, that the scaling in 
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Figure 7: The masses of the lowest vector (rho meson), pseudoscalar (pion), and singlet 
scalar states as a function of m q and in units of Air (left) and Auv (right). The rho meson, 
pion, and singlet scalar masses are shown as the solid blue, dashed red, and doted magenta 
curves, respectively. 
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Figure 8: The (real part of the) anomalous dimension at the fixed point as a function of x 
for potentials I with Wo = 3/11 (solid blue curve) and for potentials II with SB normalized 
Wq (dashed red curve). The kink lies at x = x c for both potentials. 


regime B involves A* = 7 * + 1 which depends strongly on x. We show Rey*, which 
controls the scaling exponents of the masses, as a f un ction of x for both potentials I 
and II in Fig. 8 . The kinks in the plots are located exactly at x — x c , and 7 * drops 
rapidly right above the kinks. This supports the idea that x = x c can be located 
by extracting 7 * from the meson masses on the lattice. In fact, recent lattice results 
for 7 * in the conformal window report very low values [36,57] that are in apparent 
contradiction of the curves in Fig. 8 . Recall however that the model has not been 
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tuned to fit any QCD data yet 10 . It is actually not difficult to construct potentials 
for which 7 * drops much more rapidly when x > x c . 


5. Quark mass and the chiral condensate 


Let us then discuss the mass dependence of the chiral condensate. Recall that the 
tachyon solution in the UV reads 


r(r) 

i 


= m q r (— log rA uv ) p 
+ or 3 (- log rAuv) p 


1 + 0 
1 + 0 


log rA uv 
1 

log rA uv 


(5.1) 


where a can be identified as the chiral condensate * 11 (the exact identification is studied 
in Appendix C.l), and p can be expressed in terms of the leading coefficients of the 
beta and gamma functions as p = 70 /fro = 9/(22 — Ax). 

Notice that the analysis of previous sections was restricted to the standard, 
dominant vacuum. It is known [14,23], however, that there are subdominant “Efimov” 
vacua in the QCD and walking regimes {x < x c ) which quite in general appear in 
connection to the BKT transition (see, e.g., [29,30,38]). These vacua are mapped 
to different values of the chiral condensate on the (m q , cr)-plane: all possible regular 
vacua form a spiral structure, which will be called the “Efimov spiral” below. The 
results for the spiral structure will be used to analyze four-fermion deformations of 
QCD in Sec. 6. 


5.1 Efimov spirals 

Let us first review the structure of the subdominant vacua. Including the solutions 
with finite quark mass [14]: 

• When x c < x < Xbz, only one vacuum exists, even at finite quark mass. 

• When 0 < x < x c and the quark mass is zero, there is an infinite tower of 
(unstable) Efimov vacua in addition to the standard, dominant solution 12 . 

• When 0 < x < x c and the quark mass is nonzero, there is an even number 
(possibly zero) of Efimov vacua. The number of vacua increases with decreasing 
quark mass for fixed x. 

10 One should also keep in mind that we are working in a bottom-up model, which is defined in 
the Veneziano limit whereas the lattice data was computed at finite Nf and N c . 

11 When m q is finite, the UV expansions in practice only define a up to a linear term in m q . This 
issue will be discussed below and in Appendix C.l. 

12 It is assumed here for simplicity that there is an IRFP for all positive values of x, and the BF 
bound is violated at the fixed point all the way down to x = 0. 
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The infinite tower of Efimov vacua, which appears at zero quark mass, admits a 
natural enumeration n = 1,2,3,... where n is the number of tachyon nodes of the 
background solution as we shall demonstrate below (see also Sec. 10 and Appendix H 
in [14]). A generic feature of these backgrounds is, that they “walk” more than the 
dominant, standard vacuum, so that the scales Auv and Air become well separated 
for all 0 < x < x c when n is large enough. It is possible to show that 



for any 0 < x < x c . The coefficient v will be given below in Eq. (5.8) (see also 
Appendix F in [14]). In the walking regime, one finds that 


Auv {K( n + 1)\ 

- -exp , 

Air V V • x c ~ x J 


{x ->• x c ) 


(5.3) 


for any value of n. In particular, n = 0 corresponds to the standard solution discussed 
in the previous sections, and the relation (5.3) gives the standard Miransky scaling, 
whereas for n > 0 the scaling is even faster. We also found a similar scaling result 
for the free energies of the solutions as x —> x c in [14], therefore proving that the 
Efimov vacua are indeed subdominant in this limit, and verified this numerically 
for all 0 < x < x c . In [23] it was shown that the Efimov vacua are perturbatively 
unstable (again analytically as x —>■ x c ~, and numerically for all 0 < x < x c ). 

When the BF bound is violated at the IRFP, the quark mass and the condensate 
are known to show an oscillating behavior for the (chirally broken) backgrounds 
where the coupling flows very close to the fixed point [14]. Let us first discuss how 
these oscillations arise from the tachyon EoM. 

First, take a background at zero tachyon which reaches the fixed point as r —> oo. 
Then consider turning on an “infinitesimal” tachyon. It satisfies the linearized tachyon 
EoM 



r 


2ga(A.) 

r 2 K(A*) T 


(r —>■ oo) , 


(5.4) 


where A* is the value of the coupling at the fixed point and is the IR AdS radius. 
Inserting the Ansatz r ~ r A * yields 


2ffa(A«) _ 24a(A») 

^j( A*) I / efF(A^)^q(A^) 


(5.5) 


where 

krfr(A) = V g (\) — xV f o(\) (5.6) 

in terms of the dilaton and tachyon potentials. The BF bound is thus given by 


24a(A*) < ^ 

Kff(A*)«(A*) 


(5.7) 
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When the BF bound is violated we denote 


v = ImA* 


/ 


beff ( A * ) ^ ( A * ) 


24a(A*) 


(5.8) 


In this case the asymptotic infinitesimal tachyon solution is oscillatory, 


r ~ r 2 sin(z/ log r + </>) , (r —)■ oo) . 


(5.9) 


Let us consider next specific tachyon solutions r m and where either an in¬ 
finitesimal quark mass m q or a condensate a is turned on in the UV, respectively. 
These solutions are conveniently expressed in units of Auv, and will have the asymp¬ 
totics of (5.9). We denote 



(5.10) 


(5.11) 


where the coefficients Ji* and 0* cannot be computed analytically but it is easy to 
extract them from numerical solutions. One can require that K, > 0 and — 7t/2 < 
4>i < 37 t / 2 . 

We are, however, interested in the solutions where the tachyon is small and finite. 
In this case the tachyon will eventually grow large when r ~ 1/A r ~ 1/Air, and drive 
the system away from the fixed point. The above formulas (5.10) and (5.11) then 
hold as approximations for 1/Auv -C r <C 1/Air. The solution for r > 1/Air depend 
on the details of the model in the IR. However, in order to satisfy the boundary 
conditions imposed by the good IR singularity, the tachyon must have certain fixed 
normalization and phase when expressed in IR units: 




(5.12) 


In general the IR asymptotics of the background depends on one parameter (e.g., T 0 
for potentials I [14]) but as the fixed point is approached, the dependence on this 
parameter appears only through the ratio Arv/Air whereas /Fir and 0 ir take fixed 
values (see Appendix I in [23]). 

As the fixed point is approached, the result (5.12) must match with the sum 
of (5.10) and (5.11). Therefore one finds 
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From here one can solve 


m q _ Km, sin ( 0 m - <$> a - vu) _ 2u 
Auv K m (j) a ) 

v _ Ki r sin ( 0 IR -4>m~ vu) 2 u 
Afj V K a sin (</> CT - <f> m ) 

where 

u = l og^ . (5.15) 

Air 

As u varies the equations (5.14) dehne a spiral on the (m q , cr)-plane, which has been 
studied recently at finite chemical potential in a different context (see [38]). Notice 
that v does not need to be small, ft can be verified numerically that the handedness 
of the spiral is such that its phase increases with increasing u (counter clockwise 
direction) on the (m q , <r)-plane. This means that 

sin (<f) m -<t> a )> 0 . (5.16) 


As we shall show in Appendix C.l, this is also required in order for the standard 
solution (with nonzero and nodeless tachyon) to be dominant (for x < x c so that the 
BF bound is violated and the spiral exists). 

Finally, let us point out some properties of the spiral as x —> x c from below. 
Notice from (5.8) that u = O (y/x c — x ). The approximations (5.10) and (5.11) 
are valid for 1/Auv r <C 1/Air, but they should join smoothly with the UV 
asymptotics of the tachyon in (5.1). At r ~ 1/A UV we find that 


^(rA uv ) 


-2 ™>q 


A 


K m sin (j) m i 


d 


uv 


dr L ^ 


T, 


t(^Aw) 


-2 


~ -— K m U COSCpr 


A 


uv 


(5.17) 

These estimates, and similar estimates for the solution (5.11), can be matched with 
the UV asymptotic formulas if 


K m K n 


~ sm < 


sin (j) a ~ y/x c — x , {x —> x c ) 


(5.18) 


An analogous argument shows that for (5.12) to satisfy generic IR boundary condi¬ 
tions, 


~ sin ~ yjx c - x , (x ->■ x c ) . 

Air 

We have found numerically that (with the convention Ki > 0 and —7 t/2 < </>j 
4> m = O (y/x c - x) = 4> a , 0IR - VT = O (y/x c - x) , 


(5.19) 
< 3tt/2) 

(5.20) 


(and (j) m — > 0 such that (5.16) holds) for all potentials which we have studied. 

Then the first node of the mass in Eqs. (5.14) in their regime of validity u •Cl occurs 
at 

vu = V log —^— = (j) IR - (j) a = 7T + O U/x c - x) . (5.21) 

Air 
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Figure 9: The dependence of a on the quarks mass for potentials II with SB normalized 
Wo at x = 2. Left: the spiral of physical solutions on the (m q , cr)-plane. Right: the same 
data after taking logarithms of both a and the quark mass. The red dots are numerical 
data while the blue curves are analytic fits. See text for details. 

This node is identified as the “standard” solution, where the tachyon has no nodes 
(in particular no nodes appear in the regime of validity of (5.10), (5.11), and (5.12)). 
Notice that solving Auv/Air from (5.21) results in Miransky scaling 13 of (3.5). Nodes 
at larger values of u, i.e., vu — (n + l)7r, with n — 1, 2, ... are identified as the Efimov 
solutions, where the tachyon has n nodes. 

5.2 Numerical results 

As an example, we have computed m q and a numerically for potentials II with SB 
normalization 14 for Wq at x = 2. The results are shown in Fig. 9. The blue dots are 
our data and the red curves are given by Eqs. (5.14). The blue line on the right hand 
side is a power-law fit cr ~ —Cm? q . The spiral structure is not well visible on the left 
hand plot because the distance of the curve from the origin decreases exponentially 
with increasing u in (5.14). In order to make the spiral structure visible, we have 
plotted the logarithms of (the absolute values of) a and m q on the right hand side 
of Fig. 9. Since the action is symmetric under r —> —r there is actually also another 
spiral which not shown in the left hand plot but can be obtained simply by a rotation 
of 180 degrees around the origin. 

Notice that when m q ^ 0 it is difficult to define a unambiguously in practical 
calculations, because the vacuum expectation value (vev) solution of the tachyon 
cannot be separated from the subleading terms of the source solution. This means 
that we also have to specify more carefully how K m and 4> m in (5.10) are defined in 

13 In principle it could be possible to satisfy the boundary conditions with biR = O {yjx c — x). 
We speculate that this happens in the model of [58] where no Miransky scaling was found. 

14 It is much easier to extract the Efimov spiral numerically for potentials II than for potentials I. 
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Figure 10: The dependence of a on the quarks mass for potentials II with SB normalized 
Wq. Left: results in the walking regime, x = 3.6. Right: results in the conformal window, 
x = 4. The red dots are numerical data while the blue curves are fits based on Eqs. (5.14). 
See text for details. 


our numerical analysis: we pick a reference solution which has a = 0 by definition, 
and consequently defines K m and 0 m (see also Appendix C.l). It is natural to expect 
that the <7 = 0 solution has a quark mass of O (Auv) in the QCD regime. Therefore 
the reference solution was fixed to be the standard, dominant vacuum solution with 
m q /Auv = 1- This kind choice is important in order to avoid unnatural fine tuning 
effects. 

The various coefficients in the solutions in (5.14) were determined as follows. 
Three solutions were chosen for the background for which a and m q are very small 
such that (5.9) holds as a good approximation for at least two periods of oscillation. 
The first (second) solution was tuned to have approximately zero a ( m q ) and was 
used to fit the constants A7, 0* in (5.10) [in (5.11)]. The third solution was used to 
fit the constants and 0 ir in (5.12). 

In the walking regime, extra care is needed in the choice of the reference solution 
having cr = 0. The chiral condensate is expected to have a node around such values 
of the quark mass where the normalizable and nonnormalizable terms in the tachyon 
solutions are nontrivially coupled. This happens for the standard vacuum at the 
crossover point between the regimes A and B, i.e., for for 


m q 

~ exp 

Auv 


2 AT 

V X C - X 


(5.22) 


The spiral in the walking regime (x = 3.6, close to x c ~ 3.7001) is show on the log- 
log scale in Fig. 10 (left). By studying numerically various quantities (for example the 
S-parameter, plotted below in Fig. 20) it is found that the choice m q / A uv = 3 x 10 -11 
lies at the crossover and is therefore chosen as the reference point with a = 0. 

As the critical value x c is approached, the Efimov spiral becomes “squeezed”, as 
seen by comparing plots in Fig. 9 (right) and 10 (left). In the walking regime, the 
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consecutive solutions with vanishing m q and a are close so that the spiral approaches a 
straight line in the log-log scale. This reflects our analytic results for the spiral (5.14) 
derived above in the limit x —>■ x c ~: in particular <f) m — <\> a ~ \/x c — x —* 0. 

We have also repeated the analysis in the conformal window. In this case the 
scaling of (3.10) is expected to hold for small quark mass and that a ~ m q for large 
quark mass. The results for potentials 11 with SB normalization for Wq at x = 4 are 
shown in Fig. 10 (right). The red dots are the data while the blue curves are given 
by the power laws mentioned above. 

Again we need to specify how a is extracted in the numerical analysis because 
the tachyon vev solution cannot be separated from the subleading terms of the source 
solution. In the conformal window, however, the solution is simple because the values 
of a grow fast with increasing quark mass. Therefore one can choose any reference 
solution “with a = 0” at very small m q , and the results are essentially independent 
of the choice. 

Notice that the plots in the walking regime (left) and conformal window (right) 
do not seem too different. Indeed the curves undergo a smooth transition at x = x c . 
The coefficients of the power laws are continuous at the transition, and the nodes 
where a = 0 or m q = 0 approach the origin of the spiral very fast as x —> x c from 
below. Actually the rate of the approach is given by the Miransky scaling law. 


5.3 Gell-Mann-Oakes-Renner relation 


The GOR relation can be obtained as usual by combining the results from two 
computations. The first result is the expression for the pion mass at small m q , 
which is obtained by analyzing the fluctuation equations at small m q as done in 
Appendices D and E. One finds that 15 


m\fl = 2M d N f N c tW 0 K 0 m q a 


1+0 

a 


(5.23) 


where k 0 = k(\ = 0) and Wo = V/o(A = 0). The second result is the relation 
between a and the chiral condensate as m q —» 0, which is obtained by deriving the 
renormalized on-shell action (i.e., the vacuum energy £) with respect to m q (see 
Appendix C.l): 


d£ 

(qq) = q^ = -2M 3 N f NJ 5 WoKoa 


The combination is the GOR relation: 

m lfl = ~m q {qq) 


l + O 


l + O 


"R A uv 

a 


m g A UV 

a 


(5.24) 


(5.25) 


It can be checked that the proportionality coefficient (here minus one) is correct 
for our normalization of /^, which differs by the factor Nf/2 from the standard 
normalization in chiral perturbation theory (see, e.g., [59]). 


15 


The pion decay constant A will be discussed in detail below in Sec. 7. 
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Figure 11: Numerical test of the GOR relation. Left: Both sides of the GOR relation 
plotted against the quark mass. The blue curve is rn q (qq) (in units of Ayy), while the red 
dashed curve is //m/. Right: The leading correction to the relation. The blue curve is 
the numerical data, and the thin dashed red line is a linear fit (2.5 m q / Auv). Potentials I 
with x = 1 and Wo = 3/11 were used. 

Notice that 

1. One might expect that the logarithmic terms which appear in the UV ex¬ 
pansion of the fields would result in correction that are only suppressed by 
1 / log m q in (5.23). Remarkably, the logarithmic terms completely cancel, as 
they should, since the QCD the relation holds up to linear corrections in the 
quark mass. Why this happens is shown in Appendix D. The cancellation 
of these corrections is also consistent with the combination m q (qq) being RG 
invariant [60]. 

2. The derivative in (5.24) is nontrivial, since changing the quark mass also affects 
the geometry even at m q = 0 due to the full backreaction between the flavor and 
glue sectors, which may add contributions to the derivative (see Appendix C.l). 

3. The relation is valid only in regime A: the correction terms in (5.25) become 
large at the crossover between regimes A and B. 

We have tested the GOR relation numerically for potentials I 
sides 16 of the relation are plotted as functions of the quark mass 
and good agreement is found. The subleading correction in Fig. 11 
linear in the quark mass as it should. 

6. Four-fermion operators in V-QCD 

After we have constructed the solutions on the (m q , cr)-plane it is straightforward 
to analyze the effect of multitrace deformations following the recipe of [40]. In the 

16 To be precise, we are checking equation (5.23) - the normalization factor of (5.24) was assumed 
in the computation and subleading corrections were dropped. 


at x = 1. Both 
in Fig. 11 (left) 
(right) is clearly 
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presence of such deformations, the coefficient of the source term, which was denoted 
by m q above, is no longer trivially related to the quark mass. Let us therefore denote 
this coefficient by a m in this section. Similarly, the coefficient of the normalizable 
term is denoted by /3 m . 

One can now study the following extra terms on the field theory Lagrangian: 


1 

N f N c 



d A xO + 


(-1 T9r 


n=2 


n 


d 4 xO n 


( 6 . 1 ) 


where O = qq/NfN c . We first replace O by —c a /3 m in W, where c a = 2M 3 W 0 k 0 £ 5 is 
the proportionality coefficient between cr and ( O ) in the absence of the multi-trace 
deformations 1 '. The result is the functional 

1 P ^max yi P 

wywyW'Idm] = m q c a / d A xfd m + / d A x^ . ( 6 . 2 ) 

Then the boundary conditions are obtained by setting the source term a m to the 
value [40] (see also [61,62]) 

^max 

m q +Yl , (6-3) 

n=2 


1 5W 
NfN c c a 6(3 m (x) 


and the vev is given by 


A . = „ s --<0 =. 


(6.4) 


The possible vacua can then be identified by overlapping this condition with the 
curves of regular solutions on the (a m , /3 m )-plane (shown above in Figs. 9 and 10). 

We are most interested in the case of double trace deformation, g 2 7 ^ 0 with other 
couplings equal to zero, since this operator becomes marginal at the critical point 
x = x c . Let us also set m q = 0. In this case 


92^cr(dm ■ 


(6.5) 


The overlap plot is shown for the phase with the Efimov spiral (0 < x < x c ) in 
Fig. 12. Notice that the second branch of the spiral, obtained by reflection about 
the origin, which was omitted in earlier plots (for example Fig. 9) was now also 
included 18 . For this special case where g 2 is the only nonzero coupling, as the mapping 
(ctm, f3 m ) 1 —> ( m q ,a ) in (6.3) and (6.4) implies, the change of the UV boundary 

17 For simplicity we will omit here the difficulty of defining /3 m in practice: we take it to be directly 
proportional to (qq) (so that k T = 0 in (C.23) of Appendix C.l). 

18 The solutions on the second branch are nontrivially related to those on the first branch only if 
operators with odd n have been included. 
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Figure 12: Solutions in the presence of a double trace deformation in the QCD regime. 


conditions with respect to the standard case #2 = 0 corresponds to changing the 
vertical axis to the line defined by (6.5) (examples are the red and dashed magenta 
lines in Fig. 12) while the horizontal axis is kept fixed. 

Notice that the solutions with <y 2 = 0 he on the vertical axis, and are denoted by 
the green line in Fig. 12. The green dot shows the (stable) standard solution, but 
there is also an infinite tower of additional intersection points near the origin, which 
are not visible as the spiral converges very fast towards zero. These intersection 
points give the Efimov solutions. 

In order to draw the phase diagram at nonzero g 2; we need to solve the free energy 
of each solution and find the dominant vacuum. We can start from the identity 

d£ 

= (qq) = —NfN c c a a . (6.6) 

In Appendix C.l we show that the conditions (6.3) and (6.4) are indeed consistent 
with (6.6) and that the higher order expectation values satisfy 

n(—l) n+1 d£ 


n = (O n ) = (0) n 
N f N c dg n K 1 W 

in agreement with the large N factorization of expectation values. 

Recall from Sec. 5 that the Efimov spiral can be parametrized as 

a m K m sin (0 IR - (j) a - vu) 


(6.7) 


-A-UV K m sin (0777, 0cr) 

Pm Km sin (0 m - 0 m - vu) 


—2u 


—2u 


( 6 . 8 ) 


A 3 

iV uv 


K a sin (cj) a - 
asymptotically at small values of the variable 

Auv 


u = log 


A 


IR 


(6.9) 
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Figure 13: The phase diagram of V-QCD in the presence of a four-fermion deformation (at 
zero quark mass and temperature). The blue horizontal line is a discontinuity at g- 2 = 0~. 
The dotted vertical line presents a BKT transition. 

Notice that if only g 2 is nonzero, m q and a still satisfy (6.8) (with a m ( f3 m ) replaced 
by m q (a), respectively), if the coefficients Ki and (j>i are redefined. This is seen by 
inserting the conditions (6.3) and (6.4), and is consistent with the change of boundary 
conditions simply corresponding to a new choice of axes in Fig. 12. 

By inserting this parametrization in (6.6) and integrating along the Efimov spiral 
we fold that for the solutions with m q = 0 (see Appendix C.l) 




( 6 . 10 ) 


SAAiAA sin((/> m (f) c r) 


where £ 0 is the free energy of the solution with a m = 0 = f3 m . Since sin(</> m — (j) a ) > 0, 
the dominant solution from those in the range of validity of (6.8) is that with the 
largest value of u. In the walking regime this can be seen explicitly. Namely, we 
argued in Sec. 5 that in the walking regime v ~ H\J{x c — x)/K —> 0, and that the 
solutions are found at vu ~ {n + l)7r, with n = 0,1, 2,.... Therefore 



pendix H there). 

Also quite in general the solution with largest |a| is the dominant vacuum. 
From (6.6) we see that the energy density is given in terms of the (oriented) area be¬ 
tween the spiral and the horizontal axis. For clockwise oriented spirals the minimum 
energy is reached at the solutions furthest away from the origin. 

It is then straightforward to construct the phase diagram, which is shown in 
Fig. 13 for the case of m q = 0 and nontrivial g 2 - For x < x c the diagram can be 
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found simply by analyzing the solutions in Fig. 12. When g 2 > 0 the configuration 
is qualitatively similar to that at g 2 = 0: the dominant vacuum for x c < x is the 
vacuum with the standard tachyon solution having no nodes. When g 2 < 0, the 
dominant solution is that of Fig. 12 (right) having sizeable |er|. As this solution is 
absent for g 2 > 0, there is a discontinuity (a “zeroth order” transition) at g 2 = 0—. 
It is smoothly connected, though, to the dominant solution of g 2 > 0 through the 
limit g 2 —» ±oo (where the constraint (6.5) gives a horizontal line in the plots of 
Fig. 12). There is also a subdominant solution shown in Fig. 12 (left), namely the 
continuation of the standard, dominant solution at g 2 — 0 to negative g 2 . We will 
argue below that this solution is metastable for small \g 2 \ and unstable for large \g 2 1. 

Notice that the chiral condensate is a scheme dependent quantity, and the scheme 
dependence is important, in particular, at large m q (see [46] for a discussion of the 
scheme dependence in the context of holography). Therefore the dominant solution 
g 2 < 0 appears scheme dependent at least for small \g 2 \. We have, however, found 
the behavior of generic schemes a oc m q with a negative proportionality constant for 
all choices of potentials we have tried, and independently of the precise definition of 

а. This is enough to guarantee that the phase diagram is that of Fig. 13. Actually, 
as stressed in Appendix B, the results at very large m q are essentially independent 
of the choices of the various functions in the V-QCD action. 

In the conformal window, the situation is even simpler. For g 2 > 0 there is only 
the solution with zero quark mass and the condensate. When g 2 < 0, there is again 
a solution with large |er|, analogous to that shown in Fig. 12 (right). When g 2 > 0, 
the transition at x = x c is similar as at g 2 = 0, i.e., a BKT transition, since changing 
the value of g 2 only amounts to changing the axes in Fig. 12 without affecting the 
structure of the spiral. 

One can also show that when a finite m q is turned on, the BKT transition and 
the chirally symmetric phase disappear, but the discontinuity at g 2 = 0“ remains. 

б. 1 Perturbative stability 

Finally perturbative stability of the solutions with modified boundary conditions 
could be analyzed following [23]. Here we will only discuss which of the solutions are 
expected to be unstable, and will not prove the stability of any solution. Naturally, 
the modification of the boundary conditions for the background also implies that 
the boundary conditions of the fluctuations are similarly changed as the fluctuations 
must preserve the physical value of m q . 

Let us first analyze any solutions with walking, i.e., «» 1. Recall that the stan¬ 
dard solution has been shown to be stable, while the Efimov solutions are unstable 
when g 2 = 0 [23]. The instability appeared in the scalar flavor singlet and nonsinglet 
sectors. In order to study stability of the other solutions, one should look at the 
scalar fluctuation equations. In the walking regime, they admit simple solutions in 
the UV and in the vicinity of the (approximate IR) fixed point. In fact, as argued 
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in Appendix I of [23], the fluctuation equations (for sufficiently small mass of the 
fluctuation) take the same form as the EoMs for the background. This is true both 
in the flavor singlet and in the nonsinglet sectors. Therefore, the fluctuations in the 
vicinity of the (approximate) IRFP are analogous to (3.3): 


ips(r) ^ C s r 2 sin [u log(rA uv ) + <j> s \ , 



< r < — 
A x 


( 6 . 12 ) 


where il>s is the radial wave function of any scalar fluctuation mode. The solution 
necessarily has a node if this approximation is valid for more than half a period of 
the sine function. In terms of variable u = log(Auv/Am) the node therefore appears 
for uu > 7T (where we used the fact that Air ~ A r whenever walking is present). 

Such a node of the wave function (say at zero momentum) implies an instability, 
for generic UV boundary conditions for the fluctuations. It is straightforward to 
prove this in the case of nonsinglet scalars for which the fluctuation equations can be 
cast into the Schrodinger form, and with the standard UV boundary conditions such 
that the fluctuation wave function is normalizable in the UV (see also the analysis 
in [63]). We shall sketch the proof here. 

Denote the UV normalizable (but not necessarily IR normalizable) Schrodinger 
wave function by 0 and the location of any of its nodes by r 0 . By studying the 
variation of 

/>ro(m 2 ) 

0=1 0(r)(—0"(r) + Vs{r)(f>{r) — m 2 0(r))dr (6.13) 

Jo 

with respect to the mass, we find that 


dr 0 
dm 2 


1 

0'(r o ) 2 


cr 0 


4> 2 {r)dr < 0 


(6.14) 


We see that all nodes move towards the IR as m 2 is lowered. But as m 2 —y — oo 
the solution for generic UV boundary conditions is 0 oc exp(|m|r) and has no nodes 
for r <C 1/m. Therefore all nodes must disappear either by moving to r = oo or 
by merging with other nodes. However, the Schrodinger equation does not admit 
solutions with a double node for regular potentials Vs, so merging of the nodes is 
not possible, and consequently all nodes must disappear by moving to the far IR. 
In particular, if there the wave function has a node when m 2 = 0, the node must 
move to r = oo at some negative value of m 2 . At this value, 0 is IR normalizable: 
otherwise it could not have a node in the far IR if m 2 is slightly perturbed. Tachyonic 
normalizable mode marks the presence of instability. Putting the above observations 
together, we conclude that a node of the m 2 = 0 wave function marks the presence 
of an instability. 

The above proof does not apply directly to our case because the UV boundary 
conditions are modified, and consequently the integral in (6.14) is divergent. The 
divergence can be regulated by introducing a UV cutoff at r = e, but this results in 
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extra counterterms on the right hand side of (6.14) and its negativity is no longer 
obvious. 

The proof can be fixed, however, in our case (i.e., scalar fluctuations and node 
in the region of validity of (6.12)) when x —> x c from below. This is because the 
node becomes well separated from the UV region: the location of the node satis¬ 
fies log(r 0 /Auv) 1/^x c — x as seen from (6.12), where v ~ y/x c — x. For generic 
boundary conditions, the right hand side of (6.14) is then dominated by the contribu¬ 
tions to the integral near the node. The counterterms which cancel the divergence of 
the integral are essentially independent of x for any reasonable boundary conditions, 
and therefore negligible. The rest of the proof remains unchanged. 

Let us then study the solutions with m q = 0 as x —* x c . Recall that the phase 
differences 0 ir — (f) a — tt and <j> a — (j) m will be oc yjx c — x in this limit as we argued 
in Sec. 5. The standard solution (the green dot on the vertical axis in Fig. 12) is in 
the regime of validity of the approximations leading to (6.8). It is found at 


VU = 0IR - (j) a = VT T O {y/x c ~ x) 


(6.15) 


Based on the above analysis, the onset of the perturbative instability is also expected 
at vu — Ti = O (y/x c — x). As the standard solution is stable, the critical value of u 
must be larger than that given in (6.15). Inserting this in (6.8), the critical value of g 2 
is given by the ratio am/{f3 m c a ). The factors of y/x c — x cancel leaving a O (l/A^y) 
number, which must be negative given the handedness of the spiral. The critical 
value of g 2 is therefore of the same order as the value corresponding to the dashed 
magenta line in Fig. 12. The parts of the spiral which are closer to the origin from 
the critical points (near the intersection points with magenta dots) are unstable. 

The above analysis cannot be applied in the QCD regime where x c —x is not small. 
This regime could be analyzed numerically. The natural expectation is that the 
critical value is still negative and O (l/A^y). This is supported by the fact that the 
Efimov vacua are unstable, which has been shown by a numerical computation [23]. 

To conclude this section, the phase diagram is that shown in Fig. 13: For g 2 > 0, 
the diagram is identical to that of g 2 = 0. The result is somewhat different from the 
result obtained in the gauged Nambu-Jona-Lasinio model, where the four-fermion 
operator is also slightly different, and the lower bound of the conformal window 
appears to increase with g 2 [64]. As any negative g 2 is turned on, the “standard” 
vacuum immediately becomes unstable and the dominant vacuum has much larger 
|cr|. For x < x c and g 2 < 0, the perturbation analysis suggests the the standard 
vacuum is metastable for small \g 2 \ and becomes perturbatively unstable for larger 
\g 2 \, with the critical value being O (1/Ayy). 
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7. S-parameter and current-current correlators 


The vector-vector and axial-axial correlators have the structure 


i J d A x e~ iqx (0\ T{j a ^ v \x)J b u {v \ 0 )} | 0 ) 
i J d 4 x e~ iqx (0\ T{j^ A \x)J b(A \0)}\0) 


26 


ab 


Nf 
2 6 ab 


N 


f 


{q 2 v»v - %<iv) n v (q 2 ) 
- q^h) n A (q 2 ) 
+ q^q^dq 2 ) , 


(7.1) 

(7.2) 


where 

Jl {v) = qi^ a q , Jl {A) = qi^t a q • ( 7 . 3 ) 

and t a are the generators of SU(Nf) with the normalization T rt a t b = 6 ab j 2. The 
factors 2/Nf on the right hand side were added to ensure that fly and II 4 are 
proportional to Nf. The numerical factor was chosen such that these factor are 
equal to one when Nf = 2 , which is the smallest number for which the flavor non¬ 
singlet currents are defined. This will result in the standard normalization of the 
S-parameter. The sign convention is rj^ = diag(—, +, +, +). Notice that n L vanishes 
for zero quark mass because then = 0. 

When the spectrum is discrete (i.e., 0 < x < x c or m q is finite), we may formally 
write the correlators as sums over the contributions from the meson states: 


n 4 = -4 
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q 2 


(mi P) ) 2 


(7.4) 

(7.5) 

(7.6) 


Depending on the choice of the holographic action, these series may not converge. 
This issue will be discussed below in the context of the S-parameter. The residues So 
of the “spurious” q 2 = 0 pole in (7.4) and (7.6) must identical for the pole to cancel 
in the full correlator. At zero quark mass also So must be related to the pion decay 
constant, S 0 — f 2 = (/i ) 2 , since II^ vanishes. 

The difference of the vector-vector and axial-axial correlators involves the quan¬ 
tity 

D(q 2 ) = q 2 U A (q 2 ) - q 2 U v (q 2 ) , (7.7) 

which is nontrivial only when chiral symmetry is broken. The expansion of D(q 2 ) at 
q 2 = 0 defines the S-parameter: 

zV) = s 0 -T , 2 + g 9 . + .... ( 7 . 8 ) 
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Here Sq is the same coefficient which appears in (7.4) and (7.6), and we will also 
study the higher order coefficient S 2 below. 


7.1 Correlators and the S-parameter in V-QCD 

Let us then recall how the correlators can be computed in V-QCD. As pointed out 
in Sec. 2, the vector currents are dual to the gauge fields in the DBI action (2.5). 
Following the standard approach (see [23] for additional details), we carry out the 
fluctuation analysis writing down an Ansatz which separates the spatial and radial 
dependence of the fluctuation modes in momentum space. The spatial (radial) wave 
functions of the vector, transverse axial, and longitudinal axial modes are denoted 
by V, A and V (Vv, Vbn and The radial wave functions are IR normalizable 
and satisfy the UV boundary conditions 

1 = tMcP 2 ) = Me,p 2 ) = Vr(e,p 2 ) - p(e,P 2 ) • (7.9) 


Here t/jp is the radial pion wave function. 

The terms of the on-shell V-QCD action which are quadratic in the vector fields 
can then be written as 


= -M 3 AL 


d A p 


KWP^KhP) V,e A w 2 d^ v (r,p%_ 


4 

1 


4 M 


S A = -M d N r 


(2vr] 

d ‘ P A‘MP^Al(-p) V f e A w 2 d r Mr,P 2 ) L, 


1 


S L = -M N c 


(2tt) 4 ^ 

d A p 


V a {p)V a (-p) V f e A p 2 w 2 d r 'ip L (r,p z )\ , 


(2tt) 


where 


pnv = ^ _ 


pVpV 

p2 


(7.10) 

(7.11) 

(7.12) 

(7.13) 


projects to the transverse parts of the wave functions. To compute the vector-vector 
correlators, we need the precise dictionary given in terms of the couplings to field 
theory currents: 




(7.14) 

(7.15) 


Applying the gauge/gravity correspondence with these couplings leads to the follow¬ 
ing expressions for the form factors: 


1 


<? 2 n,( 9 2 ) = ~-M 3 N,N e V,e A w 2 d^,(r,q 2 )| r=e , 


(7.16) 


where / = V, A, L. 
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Notice that for m q > 0 (so that IU is nonzero) the wave functions and i/jl 
satisfy the same fluctuation equations as q 2 —> 0. Consequently 


limq 2 U A (q 2 ) = limg 2 hU(g 2 ) (7-17) 

( j->0 <?— s-0 


and this number equals Sq of (7.4) and (7.6) which is therefore well defined. This 
equality ensures the cancellation of the “pole” at q 2 = 0 as pointed out above. 

By inserting (7.16) in the definitions (7.7) and (7.8) we obtain for the S-parameter 


5 = it M d Nf N c V f e A w 2 


d 2 


drdq 


;Mr, (f 


d 2 


drdq 


; Vv(r,« 2 


(7.18) 


- r=e, q 2 =0 


This formula is, however, not convenient for high precision numerical computations, 
because the subleading terms at r = e are only suppressed by logarithms of e in V- 
QCD, so that extremely small values of e would be needed to obtain accurate results. 
By an analysis of the fluctuation equations, it is possible to derive more convenient 
integral representations for S (see Appendix E). We find that 


D(q 2 ) = q 2 U A (q 2 ) - q 2 U v (q 2 ) = M 3 N f N c / du ^ v V f e 3A KT 2 ^ A , (7.19) 

Jo 

where u is the Schrodinger coordinate defined in Appendix E, and that 


POO 

7rM 3 N f N c J du V f e A w 2 [ip\ - -0y] g2=0 

(7.20) 

POO 

7rM 3 N f N c j ^ du V f e A w 2 [t/j 2 a - l] q2=0 . 

(7.21) 


The formula for the S-parameter is well-known in the context of simple bottom-up 
models (see, e.g., [65]). In Appendix E we write it in a form which holds for very 
generic holographic models. We used the fact that 


Vvl q 2 =0 = 1 

in order to obtain the last expression (7.21). 


(7.22) 


7.2 Numerical results for the S-parameter and f n 

Using equations (7.16), (7.19), and (7.21) it is straightforward to compute the S- 
parameter, the pion decay constant as well as the coefficients S 0 and S 2 of (7.8) 
numerically (see [23] for additional details, and [26, 66, 67] for analysis of the S- 
parameter in other holographic models). Figure 14 shows the results for S at m q = 0 
(dashed red curves) and at m q / Auv = 10” 6 (blue curves) for potentials I (left) and II 
(right). The numerical value of M 3 was fixed such that the asymptotics of the 
vector-vector correlator matches with perturbative QCD (see Eq. (C.10) in [23]). 
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Figure 14: The normalized S-parameter as a function of x for m q / Auv = 10~ 6 (blue 
solid curves) and m q = 0 (red dashed curves). Left: potentials I with Wo = 3/11. Right: 
potentials II with SB normalized Wo- 

The most striking feature in these plots is the discontinuity of the S-parameter 
in the conformal window. When x > x c , the S-parameter immediately jumps from 
zero to a O ( NfN c ) number when any finite m q is turned on. The mechanism which 
leads to this discontinuity will be discussed in detail below (from the holographic 
viewpoint), but it appears rather natural: the S-parameter is 0(NfN c ) whenever 
the geometry has the IR singularity, and vanishes only for zero quark mass in the 
conformal window where there is an IRFP instead. The result is also consistent with 
the analysis based on field theory at qualitative level [35,68]: the S-parameter is 
finite except for exactly zero mass in the conformal window. 

There is, however, one striking difference [22,23] with respect to previous results: 
the S-parameter increases with x in regime A, whereas many earlier analyses [26,69 
71] suggest that the S-parameter is suppressed in the walking regime and may even 
vanish as x —> x c . Recall, however, that the IR behavior of the potentials in the V- 
QCD action has not yet been fitted to QCD or lattice data, and such fits may affect 
the ^-dependence of the S-parameter. By analyzing the form of the fluctuation wave 
functions (see Appendix G in [23]) in the integral formula (7.21) one can indeed 
check that it is dominated in the IR (for small m q and for all values of x). This is 
consistent with the analysis of Appendix I of [23], where the same is argued to hold 
for the meson masses and decay constants. 

We have also computed the pion decay constant which is defined in terms of the 
residue of LR at the pion mass whenever m q ^ 0 (see Eq. (E.21) in Appendix E). 
The results for both potentials are given in Fig. 15, and they are also compared 19 

19 Notice that Sq is UV divergent whenever the quark mass is finite, as can be seen by inserting the 
UV expansions of the wave functions [23] in (7.16), and needs to be renormalized. The divergence 
is oc m 2 q in agreement with the one-loop field theory computation of Appendix F. At small quark 
masses it is irrelevant how the renormalization is done because the difference between all reasonable 
renormalization schemes is negligible due to the smallness of the coefficient in the divergent term. 
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Figure 15: The pion decay constant as a function of x for m q /Auv = 10~ 6 . Left: poten¬ 
tials I with Wo = 3/11. Right: potentials II with SB normalized Wo- The (normalized and 
squared) pion decay constant is shown as the blue curves, and the values of the constant 
y/So are also shown as the red dashed curves for comparison. 



Figure 16: The x dependence of the higher order coefficient S 2 for m q / Auv = ICR 6 and 
for potentials I with Wq = 3/11. Left: S 2 in units of Air. Right: The dimensionless 
product SoS 2 - 

to the constant S 0 . As expected the pion decay constant (blue curves) match with 
y/So (dashed red curves) in the QCD regime. The ratio fl/S 0 decreases fast with 
increasing x for both potentials when x > x c , suggesting that the pion decouples. 

The higher order coefficient S 2 (in units of Air) is also shown for potentials I in 
Fig. 16 (left) and the product S 0 S 2 in Fig. 16 (right). The dependence on x is similar 
as for the S-parameter when x < x c . 

We have not tried to analyze the various observables in the BZ limit x — > x-qz 
because this region is not the most interesting one from a holographic viewpoint. I 11 
general, however, the slow RG flow in the BZ limit causes that all ratios of energy 
scales to be typically oc exp[#(a;Bz — a;) -2 ]. Therefore even ratios which are expected 
to be close to one for generic values of x easily blow up in the BZ limit - in this case 
# in the above relation is small, but it is difficult to define the scales such that it 
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Figure 17: The mass dependence of the S-parameter for x = 1 (blue curve), x = 4 (dashed 
red curve), and x = 4.5 (dotted magenta curve). Potentials I with Wq = 3/11 were used. 
The thin lines are extrapolations given by fits to the asymptotic behavior. 

would be exactly zero. In view of this, it is not surprising that only the dimensionless 
quantities S and SqS -2 approach finite values in the BZ region. 

Let us then analyze the mass dependence of f n and S in more detail. Fig. 17 
shows the mass dependence of the S-parameter in the QCD-regime (x = 1, blue 
curve), in the walking regime {x = 4, dashed red curve), and in the conformal 
window (x = 4.5, dotted magenta curve). The dependence on m q is relatively mild 
for all values of x = 0 (apart from the discontinuity at m q = 0 which is only present 
in the conformal window and the fact that S varies slower as x increases which is due 
to the RG flow as discussed at the end of Appendix A). In particular the limiting 
value as m q —* oo is independent of x. In this limit the S-parameter is expected 
to approach the value N c Nf/127T from perturbative QCD (see Appendix F). Even 
though V-QCD is not expected to reproduce perturbative results in general, the 
limiting value in V-QCD is numerically close to the QCD number 1/127T ~ 0.0265. 

The dependence of f n on m q is demonstrated in Fig. 18. Again the different 
plots are in the QCD regime (x — 1, top left plot), in the walking regime (x = 4, top 
right plot), and in the conformal window (x = 4.5, bottom plot). For small m q the 
dependence is weak, but when m q / A uv 3> 1 the decay constant vanishes very fast 
with increasing m q . This signals the decoupling of the pion mode in regime C, and is 
consistent with the findings of Appendix B (see Eq. (B.45)). Actually, all low-lying 
meson states are expected to decouple for potentials I. 

7.3 Scaling of the S-parameter 

We have also analyzed the S-parameter in the limit m q —» 0. Numerical results are 
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Figure 18: The mass dependence of f n for x = 1 (top left), x = 4 (top right), and x = 4.5 
(bottom). Potentials I with Wo = 3/11 were used. 


0 S(0+)-S(m q ))/N c N f 


(S(0+)-S(mq))/N c N f 




Figure 19: The mass dependence of the S-parameter for x = 4.5 in log-log scale. The 
red dots are the data and the blue lines are power-law fits oc m^' 08 . Left: potentials I with 
Wo = 3/11. Right: potentials II with SB normalized Wq. 


shown in the conformal window (for x = 4.5) in Fig. 19, where the red dots are our 
data and the lines are given by the functions 


S(m q ) - 5(0+) 
N c N f 


Pirn? 2 , 


(7.23) 


where the parameters as well as S(e), were fitted to the data. Here it is under- 
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stood that 


(7.24) 


S(0 + )= lim S(m q ) 

TYlq —^0-b 

which is a hnite number whereas the S-parameter vanishes at zero quark mass: S'(O) = 
0. Thus there indeed is a discontinuity at m q = 0. For both potentials (3 2 — 0.08 fits 
the data very well. 

In order to understand the power law in (7.23), it is useful to first discuss in 
more detail how the discontinuity at m q = 0 arises. The mechanism is the same 
which was studied in detail in the case of m q = 0 and x — > x c in section 6 and 
Appendix I of [23]. Here instead x > x c and m q — y 0. In both cases the RG flow 
of the coupling approaches the fixed point (A = A*) but misses it finally (due to 
the hnite tachyon). For such hows, it is useful to divide the background to the UV 
and IR sections, having A < A* and A > A*, respectively. Considering hows which 
get closer and closer to the hxed point, the S-parameter can be computed more and 
more precisely in terms of the IR section. The IR part takes a hxed shape in this 
limit, explaining the hnite value of the S-parameter. For exactly zero quark mass 
the IR section of the background becomes disconnected from the UV section and is 
therefore not present in the physical vacuum solution, which now ends at the IRFP. 
This is rehected in the vanishing value of the S-parameter. 

By using similar arguments, we can also sketch how the power law in the mass 
dependence arises. It is understood to be the leading “perturbation” of the IR back¬ 
ground due to the fact that the hxed point was not reached exactly. In the conformal 
window, the how towards the hxed point is given by 

A ~ A* — C\ (rAuy) 6 , (7.25) 

where 6 is related to the dimension of the TrF 2 operator at the hxed point. It can be 
computed as the derivative of the holographic beta function at the hxed point (when 
the tachyon is set to zero exactly) [14,72], One hnds that 

/ 9IAA 2 

A FF -4 = h= W4--^-2 , (7.26) 

where V) are the coefficients of the Taylor expansion of V e g at A = A*: 

V e g = V g - xV f0 = V 0 + V 2 (A - A*) 2 + • • • . (7.27) 


Notice that V2 < 0. Flow toward the hxed point ends when r ~ 1/Air. The difference 
of A with respect to the hxed point value when this happens is given by 


A* — Air = A 


A(r = 1/Air) ~ 





(7.28) 


where (A. 14) was used to obtain the last expression. This difference controls the 
deviation of the IR section of the background from its limiting shape as m q —>- 0 
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m q 

Auv 


Figure 20: The mass dependence of the S-parameter in the walking regime in log-log 
scale. The red dots are the data and the blue lines are power-law fits. Left: potentials I 
with Wo = 3/11 and with x = 4. Right: potentials II with SB normalized Wo and x = 3.6. 


and correspondingly the deviation of the S-parameter from the limiting value 5(0 + ). 
Indeed by using the explicit expressions for the potentials one obtains 

r 

— ~ 0.0780 , (x = 4.5) , (7.29) 

a value in a very good agreement with the fit from above, shown in Fig. 19. 

Notice that the difference A* — Air, not only controls the corrections to the S- 
parameter at small mass, but also to other quantities that can be defined in terms of 
the IR section of the background. Examples are decay constants and meson masses in 
IR units, which are therefore expected to have qualitatively similar m q dependence to 
the S-parameter at small masses 20 . Indeed, the power law of (7.28) agrees with that 
found for the scaling corrections to generic correlators by analyzing the Wilsonian 
RG flow perturbatively in the vicinity of the fixed point [73]. 

The mass dependence of the S-parameter can be analyzed similarly in the walking 
regime (x —> x c ~). The above calculation is approximately valid in regime B, i.e, 
when the quark mass controls the amount of walking. One can use (3.9) together 
with (7.28) to obtain 

A - “ Air ~ (yy) 2 ■ ( exp 


2K 


\Jx c - x 


<C 


m n 


A 


< 1 


uv 


(7.30) 


Therefore one should effectively take A* —> 2 in (7.28), meaning that the power is 
continuous over the conformal transition at x — x c , as A* = 2 at the transition. For 
even smaller m q , i.e., in regime A, the mass term of the tachyon can be treated as a 
linear perturbation to the whole background, and therefore the mass dependence of 
the S-parameter is linear. Both the scaling of (7.28) and the linear dependence can 
be seen in the numerical results in the walking regime, see Fig. 20. 

20 We have found numerically that ratios of masses or decay constants typically follow the scaling 
of (7.28) more accurately than their values in IR units. 
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Figure 21: S-parameter computed from the spectral representation compared to the exact 
value. Left: the dependence of the value of S on the mass cutoff in the series. The thin 
horizontal lines are the limiting values which match with direct computation of S. Right: 
the difference between the number obtained from the series and the exact value as a function 
of the cutoff. The blue, dashed red, dotted magenta, and dotdashed green curves have 
x = 1, 3.5, 4, and 4.5, respectively, in both plots. The filled circles and boxes denote the 
masses of the lowest five vector and axial mesons, respectively. 


7.4 Convergence of the spectral representation 

Finally we will analyze the spectral representation of the S-parameter in order to 
understand better why it increases with x in the QCD and walking regimes. By 
inserting (7.4) and (7.5) in the definition for S, one obtains 


n— 1 


00\2 


c rn 

(mi V) ) 


(mi A) ) 2 / 


(7.31) 


This series may, however, be ill defined. For example, the decay constant approach 
asymptotically constant values whereas the masses have linear trajectories m 2 ~ n 
for potentials I (see Appendices E and F in [23]). One can check that (7.31) is 
convergent thanks to the asymptotic cancellation of the vector and axial terms, but 
it is not absolutely convergent, and therefore the result may depend on the ordering 
of the terms. The definition of (7.31), where the states are ordered in terms of their 
excitation numbers n, would work for potentials with linear trajectories if the slopes 
of the vector and axial spectra are the same. This is not 21 the case for potentials I, and 
consequently (7.31) is incorrect. The solution to this issue is simple: the contributions 
should be ordered by the meson masses rather than by the excitation numbers. 

It is straightforward to verify numerically that (7.31) converges towards the S- 
parameter if the terms are ordered according to the masses. It turns out to be 

21 The question whether the slopes should be the same in QCD is unsettled [74], but usually 
they are assumed to be, which can be obtained for potentials I by changing the IR asymptotics of 
w{\) [23], 
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convenient to define a mass dependent cutoff, 

°° f t f(y )\2 ( a a )\2 \ 

S(m cnt ) = 4vr ( ( V ) /cut(mi y) ) - ■ (J /cut(^ A) ) J (7.32) 

n=l \( m ™ ) ) / 

where 

fcut(m) = ^ ^1 - tanh ■ (7.33) 

A smooth cutoff function was chosen instead of a step function because it im¬ 
proves convergence drastically. The convergence also means that the value of the 
S-parameter is determined through the dynamics in the deep IR, because the same 
holds for the masses and decay constants (the argument can be made precise in the 
walking regime, see Appendix I in [23]). 

The convergence of the regulated series (7.32) towards the S-parameter is demon¬ 
strated in Fig. 21. The resolution of the cutoff function was fixed to the mass differ¬ 
ence of the two lowest vector states: 8m = ^ —m^ ] . The speed of the convergence 

is best visible from the right hand plot, which shows 

AS = S - S(m cut ) (7.34) 

as a function of the cutoff. The convergence is exponential for all values of x , but 
becomes significantly slower as x increases and one moves from the QCD regime to 
the conformal window. The slowness of the convergence is not due to changes in the 
spectra. To show this, the filled circles and boxes marking the masses of the lowest 
five vector and axial mesons, respectively, were added in each plot. It is seen that 
the spectrum changes relatively little with x in unit of Air. 

Finally let us try to extract the reason for the increase of S with increasing x in 
the region of low values of x from the plots of Fig. 21. Notice that the curves for x — 1, 
3.5 and 4 essentially overlap at small values of m cut in the left hand plot. The curves 
for x = 1 and x = 3.5 deviate from that of x = 4 as m cut increases, and after deviating 
rapidly saturate to the final value of S. Therefore it appears that contributions to the 
S-parameter are roughly mass-independent up to a saturation scale, which increases 
with x. In order to have a good estimate for the S-parameter, a growing number 
of terms need to be included in the sum with increasing x, whereas the individual 
terms in the sum are of roughly constant size. Therefore the increase of S with x in 
the QCD and walking regimes can be seen to be due to slower convergence of the 
sum. In the conformal window, i.e., for the curve with x = 4.5, something different 
happens. The convergence of the sum is even slower, but the contributions at fixed 
m cut are suppressed, resulting in the decrease of the S-parameter with increasing x. 

8. Finite temperature phase diagram 

The finite temperature phase diagram has been studied in detail for IHQCD in [75] 
and for V-QCD at zero quark mass and at small values of the quark mass in [18]. 
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Here this study is extended to large values of m q as well as very high values of x. 
The code for constructing solutions at finite temperature is available at [76]. Some 
extra tricks are necessary in order to obtain reliable results in the BZ region and at 
very large m q (see Appendix G). 

First recall the generic structure of the (x,T) phase diagram [18], which was 
already reviewed in Sec. 2. At zero quark mass, there is a first order deconfinement 
transition in the QCD and walking regimes, but there is also the possibility (depend¬ 
ing on the choice of potentials and the value of x) of a chiral symmetry restoration at 
a separate second order transition. In the conformal window, there is a continuous 
phase transition at zero temperature. When a finite quark mass is turned on, chiral 
symmetry is always broken, and the second order chiral transition will become a fast 
crossover when the quark mass is small, and completely disappear at larger quark 
masses. The system is in a tachyonic thermal gas (TG) phase at small temperatures 
for all 0 < x < xbz- As the system is heated, there is a first order deconfinement 
transition to the high temperature phase, which is implemented through a transition 
from TG phase to the black hole (BH) phase in holography. 

The existence of the deconfinement transition requires an order parameter. While 
QCD at finite N c and Nf has no order parameter related to deconfinement, the 
pressure acts as an effective order parameter at large N c . This is clear in the’t Hooft 
limit, where the number of degrees of freedom (and consequently the pressure) is 
O (Af?) in the low temperate phase and O (AT?) in the high temperature phase. In the 
Veneziano limit the number of degrees of freedom is of the same order in both phases, 
but the phase transition may still be identified as a discontinuity of the pressure. In 
fact, the pressure of the model is still exactly zero in the TG phase, because our 
approach does not capture the contributions corresponding to loops of pions (as well 
as mesons with higher masses) in this phase. Including these contributions in the 
model would affect the critical temperature, and potentially even alter the order of 
the transition [21]. 


8.1 Scaling laws at finite temperature 

The critical temperature has nontrivial dependence on rn q , which can be analyzed 
analytically. Notice that the temperature brings in an additional energy scale with 
respect to the zero temperature solutions. Another difference is that the definition of 
the standard reference scale Air cannot be extended to the BH solutions in a natural 
manner (because the geometry now ends at a horizon rather than an IR singularity). 
Therefore it is understood that Air is defined below through the TG solution (or 
equivalently through the zero temperature solution at the same values of m q and x). 

Let us first discuss the mass dependence of the critical temperature, which can 
be inferred by using the results from [18] and from Sec. 3. The temperature of the 
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black hole can be related to the metric through the formula 


1 / r h dr y 1 

Aire 3Ah \J 0 e 3A h) J 


( 8 . 1 ) 


where Ah and r> t are the values of the scale factor and the bulk coordinate at the 
horizon. In V-QCD models T(r^) has a nontrivial minimum (for tachyonic BHs) and 
the transition takes place at the scale of the minimum. Indeed, the entropy density 


s BH = l-A/Wy 34 ' 


( 8 . 2 ) 


decreases monotonically (and fast) with r^, as suggested by the UV and IR (zero 
temperature) expansions of A(r) and as can be verified numerically. Further, the 
geometry of the BH solution approaches smoothly the TG solution as —y oo, and 
therefore p TG = lini^oo p B u- Integrating p' BB {r h ) = s(r h )T'(r h ), a node p B h = 0, 
and consequently a first order phase transition, is found near the minimum of T(r^). 

In conclusion, one should locate the minimum of T/r/J in order to determine the 
scaling of T c . When the geometry is close to AdS, (8.1) implies that T ~ l/r h . This 
result holds both in the UV asymptotic region and when there is an approximate 
IRFP, i.e., walking. When the quark mass is large, there is also an approximately 
AdS region where the flavors are already decoupled but the dilaton A is still small. 
In summary, 

T ~ l/r h , (r h « -±-\ . (8.3) 

For r'h 1/Air the temperature increases with as seen by studying the IR ex¬ 
pansions (see [18]). Therefore, the minimum of T(r/ l ) takes place at rv j 1/A m . By 
continuity, (8.3) implies that 

T c ~ Air (8.4) 

for all rrtq > 0 and 0 < x < x B z- 

The scaling results in units of Arv immediately follow by using the results from 
Sec. 3: 


In regime A, T c rs./ Air rs_/ Arv for small values of x and 


T, 


A 


IR. 


A 


UV 


A 


UV 


exp 


K 


\Jx c - x 


(8.5) 


as x —» x c from below. The dependence of T c on x at m q = 0 [18] is in 
qualitative agreement with analysis based on field theory (see, e.g., [77]), and 
the dependence on m q is expected to be a linear perturbation. 


• In regime B, 

T c ^ j m q 
A uv ~ V A uv 


( 8 . 6 ) 
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when x < x c and 


T t 


1 

A* 


when x c < x < xbz- 


A 


uv 



(8.7) 


• In regime C, 


T c ~ f rn q y- bo/ ^ M 
A UV ~ V A uv/ 


( 8 . 8 ) 


In addition to the phase transition, also various crossovers can be identified as 
the maxima of the interaction measure 


e — 3 p Ts — 4p 




rp4 


(8.9) 


As it turns out, such crossovers reflect the different regions of the zero temperature 
geometry. This can be understood by approximating the horizon as a sharp cutoff 
added on the zero temperature background. Substituting an AdS metric in the 
formulas (8.1) and (8.2), one finds that the interaction measure vanishes. 

First, there is the crossover which marks the transition from the quasi-conformal 
or walking phase (with approximate IRFP) to the asymptotic UV phase [18,78]. Such 
a crossover is found whenever there is walking, i.e., in regime B, and the x —> x c edge 
of regime A. The UV asymptotics is valid for r -C 1/Auv, and the flow from the UV 
fixed point to the IRFP is characterized by A uv . Consequently, using (8.1) and (8.3), 
the crossover temperature is expected to be 


7co,qc ~ 


A 


uv 


( 8 . 10 ) 


independently of m q . 

Second, there is a crossover at large quark mass, corresponding to the transi¬ 
tion from the region where the quarks are decoupled to the UV asymptotic region. 
The decoupling of the quarks takes place at r ~ l/m q as pointed out in Sec. 3. 
Consequently, the crossover temperature is given in terms of the quark mass: 

T CO}mq ~ m q . ( 8 . 11 ) 

When T -C m q the quarks are effectively decoupled, and therefore the thermo¬ 
dynamics is the same as for pure YM (that is, the x —» 0 limit of V-QCD). Notice 
that A uv , however, is defined in terms of the UV asymptotics, i.e., effectively at 
infinite energy, and different from that of YM even as m q —y oo: for YM, Auv ~ Air 
but in the limit of large m q these scales are related through (3.14) instead. Conse¬ 
quently, in order for the thermodynamics to smoothly approach YM thermodynamics 
as m q —> oo, dimensional quantities should be expressed in units of Air or T c rather 
than in units of Auv- 
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Figure 22: The dependence of the critical transition temperatures (solid blue curves) and 
various crossover temperatures (thin dotted red curves) on x and m q . In the top left plot 
we also included the m q = 0 data, and the critical temperature of the second order chiral 
transition which is shown as the thick red dashed curve near x = x c . 

8.2 Numerical results 

Let us then illustrate the dependence of the various critical temperatures on x and 
m q numerically. The basic features of the phase diagram at small quark mass are 
demonstrated in the top-left plot of Fig. 22. The first order “deconfinement” transi¬ 
tion temperatures for potentials I at zero and at tiny (10~ 6 ) quark mass are shown as 
functions of x on the logarithmic scale. The first order transitions are shown as blue 
curves. The curves overlap at small x, but as the conformal transition is approached, 
the curves become separated. The lower curve (which overlaps with the red dashed 
curve and is therefore not well visible) is the transition temperature at m q = 0 which 
goes to zero with Miransky scaling as x —» x c . For x > x c there is no transition when 
m q = 0. When a tiny quark mass is turned on the transition (upper blue curve) is 
present for all values of x. the critical temperature T c decreases with x inside the 
conformal window 22 . 

The dashed thick red curve is the second order chiral restoration transition, which 
also shows Miransky scaling as x —» x c and is absent for x > x c . This transition only 

22 In [18] also a region where T c increased with x was seen at high x ~ 4.5 (see Fig. 27 there), but 
this effect turns out to be due to the UV cutoff being too low in the numerics. 
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Figure 23: Pressure and interaction measure as a function of temperature for small (left 
hand plot) and large (right hand plot) quark mass. On the left, the solid blue, dashed 
red, and dotted magenta curves have m q / Auv = 0, ICO 12 , and ICO 10 , respectively. On the 
right, the solid blue curves are the YM thermodynamics in our model, while the dashed 
red and dotted magenta curves have m q / Auv = 10 10 and m q = 10 s , respectively. We used 
potentials I with Wq = 3/11 and x = 4 (except for the YM curves which correspond to 
x = 0). 

exists for rn q = 0 in the walking regime. 

The dependence of T c on m q is also demonstrated numerically for V-QCD in 
Fig. 22. The top-right, bottom-left, and bottom-right plots are in the running regime 
(x — 1), walking regime (x — 4), and in the conformal window (x = 4.5), respectively, 
and T c is given by the blue curve in each plot. The power laws are in agreement with 
the above formulas. For example, in the regime C with large m q we find that the 
exponent of (8.8) is 1 — bo/b™ = 2x/ll — 0.182 at x — 1 and 0.727 at x = 4, which 
is consistent with the plots. 

The crossover between the quasiconformal and UV regions at T = T coqc is seen 
as the horizontal lines in the bottom row of Fig. 22. The ratio T COjqc /Auv is constant 
as expected, but the value of the constant deviates significantly from its expected 
value, i.e., one. This happens because Auv deviates from the scale of the UV RG 
flow at large values of x, as was explained above in Sec. 4. The crossover due to the 
decoupling of the quarks at large m q at the temperature T ~ m q is best visible in 
the top-right plot of Fig. 22. 

In addition, in large part of the parameter space there is also a separate maximum 
of the interaction measure in connection to the first order transition. This kind of 
maxima have also been included as thin red curves in Fig. 22, and can be found close 
to the blue curves denoting the transition temperatures. 

Recall that in regime C the mass gap of the mesons does not have the expected 
behavior m gap ~ 2 m q in V-QCD (if Sen-like tachyon potential is assumed). Never¬ 
theless, the crossover due to the decoupling of the quarks takes place correctly at 
T co ,m q ~ m q - This is seen as a consequence of the decoupling of the mesons having 


64 













an unphysically low mass, which was discussed above in Sec. 4. 

Finally let us study in some details of the thermodynamic functions at very 
small and very large m q . We plot the (normalized) energy density, pressure, and 
interaction measure for potentials I with x — 4 in Fig. 23 (left). The solid blue, dashed 
red, and dotted magenta curves have m q /Auv = 0, 10 -12 , and 1CT 10 , respectively. 
The value of x = 4 was chosen such that the model shows the separate second 
order chiral restoration transition, which appears as a kink in the energy density 
and the interaction measure of the functions at zero quark mass. The vertical thin 
dashed black lines mark the locations of the phase transitions. As a tiny quark mass 
(m q /Auv — 10~ 12 ) is turned on, the second order transition turns into a crossover. 
Indeed the dashed curves follow closely the solid curves except for very close to the 
kink, where the curves have a smooth behavior instead. At m q /Auv = 10~ 10 a much 
larger deviation is seen already. 

The approach to YM theory can be seen by studying the thermodynamics at 
very large quark mass. Thermodynamic functions at sizeable m q are shown for the 
same choice of potentials and compared to the results for the limit of YM theory 
(x —)■ 0) in Fig. 23. As the quark mass increases, the quarks are decoupled and the 
thermodynamics near T = T c is expected to converge to that of YM theory. Indeed, 
the (dashed red) curves for m q /Auv = 10 10 lie essentially on top of the (solid blue) 
curves of the YM thermodynamics 23 . For smaller quark mass (m q /Auv = 10 8 , dotted 
magenta curves) a much large deviation is seen already. 
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Appendix 

A. Energy scales from tachyon solutions 


Here we shall demonstrate how the various dynamical energy scales of holographic 
QCD in the Veneziano limit arise from the background solutions for the tachyon. 
We will work with the V-QCD action, but as stressed in the main text, most of the 
results are universal, i.e., independent of the details of the action. 


A.l Tachyon solutions 

In order to analyze the scales at finite quark mass it is essential to recall what the 
tachyon solution is in different regimes [14]. For the discussion here it is enough to 
keep track of the power laws in r so we will drop the logarithmic corrections. When 
the tachyon is small, it satisfies a linearized EoM. Up to irrelevant terms, we have 

t" + 3 A't' — e 2 A m 2 T ~ 0 , (A.l) 


where the scale factor has the AdS behavior near UV and IRFPs, A ~ log(£/r), and 
the tachyon bulk mass is for V-QCD 


m 


2 

T 


2 a 
n 


(A.2) 


but its precise form is not important for our arguments. The bulk mass is expected 
to take roughly constant values near the fixed points. In the deep UV one requires 
that —m 2 ^ 2 ~ 3, leading to 


r(r) 

£ 


~ m q r + or 3 , 



(A.3) 


where a is proportional to the chiral condensate. Walking takes place if the coupling 
A flows very close to an IRFP but the tachyon is nonzero (and will eventually drive 
the system away from the fixed point in the deep IR). Walking can happen either 
right below the conformal window or in the conformal window (if the quark mass is 
small but finite). Let us denote 


A*(4 - A*) 


-£ 2 m 2 


(A.4) 


where A and m T) * are the AdS radius and tachyon mass at the IRFP, respectively. 
In the vicinity of the fixed point, and if x < x c so that the (squared) tachyon mass 
is below the BF bound, £ 2 m 2 * < —4, the solution to (A.l) behaves as 

r(r) ( 1 1 A 

—— ~ C w (r A uv ) 2 sin [v log(rA uv ) + 4>] , -r— < r < — , (A.5) 

t V A uv a t y 
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where v = ImA*, If x c < x < xbz we find instead 


—^ — C m (r A uv ) a * + C a (rA uv ) 4 A * 
= C m (rAuv) 1+7 * + C a ( r Auv) 3 7 * 


Auv 


1 



1 


) • (A.6) 


In the deep IR, i.e., for r 1/A r , the tachyon and the dilaton have calculable IR 
asymptotics which depend on the details of the action [23]. These asymptotics shall 
not be needed here. 

A.2 Scaling results 

The dependence of the various scales on the quark mass can now be found by requir¬ 
ing continuity (and the continuity of the derivative) of the tachyon solution. More 
precisely, one should require that the dominant and subdominant tachyon solutions 
are both continuous, but in the cases studied here this is equal to requiring that 
continuity of the derivatives. 

A.2.1 QCD regime 

The QCD regime is defined by 0 < x < x c and x c — x > 1. Let us first consider the 
case of small mass, m q j Auv 1. In this case there is (almost) spontaneous chiral 
symmetry breaking, which in V-QCD means that the (vev term of the) tachyon grows 
large at some value of the radial coordinate (energy scale) and triggers a nontrivial 
IR geometry [14]. The fields obey their IR asymptotics exactly when the tachyon 
is large. Therefore both the UV and IR scales are comparable to the scale of the 
tachyon, Auv ~ Air ~ A r , and a ~ Ay V . In the dynamic AdS/QCD models [33] 
similar results are obtained by introducing an IR cutoff where the tachyon grows 
large. From (A.3) one indeed sees that the mass term of the tachyon is suppressed 
with respect to the vev term for r <C 1/Auv which implies that m q can be treated 
as small perturbation. 

When the quark mass is large, m q /Auv 1, the tachyon grows large at small 
r ~ 1 /m q and therefore A r ~ m q . The RG flow of coupling is determined by the 
QCD beta function for r <C 1/A r and by the YM beta function for r 1/A r as the 
growing tachyon decouples the flavor sector. Explicitly 24 , 


A 


A 



fro log(rA uv ) 
1 


1 



(A.7) 


(A.8) 


24 Precise definitions of the scales in the UV asymptotics of A would require including the higher 
order terms, but this does not affect the results. 
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Since b^/b™ < 1, requiring continuity leads to the counterintuitive result that the 
scale of the IR expansions is larger than that of the UV expansions: 
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uv 
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IR 


rrin 
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uv 


b 0 /b™-l 


m n 
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IR 


m n 
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uv 




(A.9) 


This is however consistent with the fact that decoupling the flavor degrees of freedom 


makes A to run faster than it does at small m q . 


A.2.2 Walking regime 

The walking regime is defined by x < x c and x c — x -C 1. For very small quark mass 
it is expected that the scaling is the same as for m q = 0 and that A r ~ Air - the IR 
geometry is linked to the growth of the taehyon as in the QCD regime. Requiring 
continuity 25 with generic IR boundary conditions at r ~ 1/A IR gives C w ~ A 2 R /A RV 
in (A.5). Continuity of (A.3) and (A.5) at r ~ 1/Auv further leads to a ~ A 2 R A UV , 
and from (A.3) it is seen that the quark mass term is small if 


m q 

Auv 


<C 


Air 

Kv 


~ exp 


2 K 

y/ x c - X 


(A.10) 


where we used (3.5). 

When 

m„ 

« t- 5 - « 1 , (A.11) 

A UV 

the amount of walking is controlled by the quark mass. We still require that the IR 
geometry is exactly where the taehyon is large so that A r rs-/ Air and C w A?r/A r y, 
but now continuity of (A.3) and (A.5) at r ~ 1/A UV yields 


exp 


2 K 


y/x c ~ X 


™ q _ Am_o_ 

AUV Ayy A^jy 


(A-12) 


and there is no Miransky scaling. 

When m q /A\jy 1, the coupling never flows close to the fixed point and the 
walking behavior is therefore absent. One finds the same results as in the QCD 
regime above. 


A.2.3 Conformal window 

Conformal window is the regime with x c < x < Xrz- Let us first assume that we 
are not in the perturbative BZ regime, Xrz — x > 1. When the quark mass is small, 

25 There is also the continuity of the derivative involved - this fixes the phase of the slow oscillations 
of the taehyon. Taking this into account one can derive the Miransky scaling of (3.5) as detailed 
in [14] and also here in Sec. 5. 
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m q j Auv <C 1, one again requires that A T ~ A IR and that there is walking. Continuity 
at r 1 /Aie gives 



(A.13) 


for the coefficients of (A.6). Matching the leading and subleading terms of (A.3) 
and (A.6) at r ~ 1/Auv one further obtains 




4—A* 
A* 


(A.14) 


When m q /Auv 1, the flow does not become close to the fixed point, and the 

results are the same as in the QCD and walking regimes. 

A.2.4 BZ regime 

Let us then discuss the BZ regime {x < x R z and xbz — x -C 1). One might expect 
that A r ~ m q independently of the value of m q since A* approaches one in this limit. 
We have, however, defined m q asymptotically in the UV, and the UV RG flow of the 
quark mass is singular in the BZ limit. Therefore m q and A T are not simply related 
as x —* xbz ~ see the top-left plot of Fig. 3, where A r increases with x in the BZ 
region instead of approaching m q . Therefore, we will use the scale A T in the analysis 
below. 

When the quark mass is small, the coupling flows toward the BZ fixed point when 
r < 1/A r , and like in YM when r 1/A r . Since A* is small, one can use (A.8) to 
describe the YM flow. Requiring the flow to start at A* when r ~ 1/A r gives the 
exponential scaling 



(A.15) 


In the opposite limit, i.e., large quark mass, the reasoning leading to (A.9) ap¬ 
plies, if one replaces m q by A r in the formulas. There is, however, also a subtlety 
in the definition of A uv . In [14], Auv was defined essentially as the scale where b 0 A 
becomes O (1). However in the BZ regime the coupling reaches the fixed point well 
before reaching the value l/b 0 ~ 1/(xbz — x). Consequently, Auv is exponentially 
suppressed with respect to the true characteristic scale of the RG flow in the UV. 
Let us instead denote by Auv the scale defined as [52] (roughly corresponding to the 
scale where X/b 0 becomes O (1)): 



(A.16) 


This formula is the RG flow given by the two-loop BZ beta function, and it is re¬ 
produced in V-QCD up to correction suppressed by xbz — x (for the whole flow 
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when r = 0 and in the UV for all backgrounds), whereas the geometry is AdS, 
A ~ — log (r/£), up to highly suppressed O ((xbz — X Y) corrections. Therefore (A.9) 
is a better estimate after the replacements Auv —;> Auv and m q —* A r . 

In order to take into account the extremely slow RG flow in the BZ region, the 
conditions for the validity of (A.15) and (A.9) should actually be written as 


m n 


A 


uv 


<C 1 and 


m n 


A 


uv 


> 1 


(A. 17) 


respectively. Here b\ is the NLO coefficient of the QCD beta function and b^/bi ~ 
(^bz — x ) 2 - Notice that the slowness of the RG flow is already visible at not so large 
values of x (as demonstrated by the plots in the text, see, e.g., Fig. 17), because 
we use parameters m q and A uv defined asymptotically in the UV. Indeed the equa¬ 
tions (7.28), and (A.9) depend on the ratio m q /Auv through combinations of the 
type (m q /Aijv) 0{xBZ ~ x) . 


B. Schrodinger potentials and mass scales 

Let us then analyse the behavior of mass gaps and mass splittings (separation of 
the lowest bound state masses) in V-QCD. For the flavor nonsinglet fluctuations this 
is rather straightforward, as the fluctuation equations can be transformed into the 
Schrodinger form. The singlet fluctuations are more involved, because there is in 
general nontrivial mixing between the meson and glueball states. Only in the probe 
limit x — > 0 and in the limit of large quark mass m q —>■ oo the glueballs and mesons 
are decoupled. We will restrict here to the nonsinglet states and study the scalar 
singlet states only numerically in Sec. 4.3.3. 

For every flavor nonsinglet sector (vectors, axial vectors, pseudoscalars, and 
scalars) the fluctuation equation can be written as 

-</>» + V s (u)(P(u) = m 2 (j)(u) , (B.l) 

where Vs(u) is the Schrodinger potential, and m 2 is the mass of the fluctuation. This 
form is obtained after a coordinate transformation defined by 

^ = \A + e~ 2 A KT ' 2 = G (B.2) 

dr 

and u{r = 0) = 0. The Schrodinger potential is given by 

v s (u) = =4L + H(u) , (B.3) 

where E(w) and H(u) are different functions for each sector and can be explicitly 
expressed in terms of the potentials (see Appendix A in [23]). For example, in the 
vector sector we find that 

= y/V f we A/2 , H v = 0 , (B.4) 
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and for the axial vectors = Sy but 


Ha = 


e 2 A T 2 K 
w 2 


(B-5) 


The vector/axial decay constants are given by (see Appendix F in [23]) 


f 2 = M 3 N c Nf 


(u) [d u (<l)(u)/E(u))]‘ 


mt 


(B-6) 


Notice that by using (B.3) the Schrodinger equation (B.l) may be written as 


-d„ 


( u)d u 


4>n{u) 


(u 


= [ m n - H ( U )\ = {U)M U ) 


(B.7) 


Integrating this over u and inserting in (B.6), we obtain 

M 3 N c N f 


f = 

J n 


mt 


du [ml - H(u )] E(u)4> n (u) 


(B-8) 


B.l Small quark mass 


Let us then analyze the dependence of the bound state masses on the quark mass 
in the three regimes of Fig. 2. In regime A, the quark mass is a small perturbation 
and the results from [23] can be used directly. The potential Vs can be computed 
both when r <C 1/A UV and r 1/Air by using the asymptotic expansions of the 
background. The results in the UV are independent of the potentials of the V-QCD 
action: 


14 - 


Uuv 

u 2 


u <C 


A 


uv 


(B.9) 


where Uyv = —1/4 for the pseudoscalars and 3/4 for other sectors, and we used the 
fact that u ~ r in the UV region. 

In the IR the asymptotics of the tachyon, and consequently the coordinate depen¬ 
dence of the Schrodinger potential, strongly on the choice for the potential functions 
in the V-QCD action. All regular choices considered in Appendices D and E of [23] 
lead to a confining potential, which grows as a function of u in the IR regime. 

If Ayy ~ Ai R , as is the case at small x and m q , we immediately notice that 
the Schrodinger potential has its bottom for all sectors expect for pseudoscalars at 
r ~ 1/Auv ~ 1 /Air. The value of the Schrodinger potential is O (Aj R ) near the 
bottom, as can be estimated from (B.9) above by requiring continuity at u 1/Air, 
and therefore the mass gap is ~ Air. Similarly it can be seen that the mass splittings 
and vector/axial decay constants are (9 (Air). This is rather evident as only one 
energy scale enters the definitions above, both the tachyon and A are (9(1) when 
r ~ 1/Air, and no cancellations are expected in the formulas. 
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The pseudoscalar sector is special because of the pion modes which obey the 
GOR relation, as shown in Sec. 5.3. The excited pseudoscalar states appear at mass 
scale Air, as can be seen numerically. 

In order to study regime B and the remaining part of regime A (near x = x c where 
walking is seen) we need to check the Schrodinger potential for 1/Auv -C r <C 1/Air, 
i.e., in the near-conformal region. The tachyon remains small also here and r ~ u. 
The Schrodinger potentials have been derived in [23] and read 






«m« 


A 


IR 


(B.10) 


where v w = 3/4 for vectors and axials, and depends on the anomalous dimension A* 
at the (approximate) fixed point for scalars and pseudoscalars. For scalars v w < —1/4 
when x < x C: whereas it is positive for pseudoscalars. Therefore the coefficient in the 
scalar potential is critical [28], and might potentially lead to an instability or a light 
state (see Sec. 5.3 of [23]). Numerically it is seen, however, that this is not the case 
and the spectrum of scalars is not qualitatively different from that of the vectors, for 
example. 

We see from (B.10) that Vs has similar dependence when the coupling constant 
walks as in the UV region in (B.9), for the vector and axial states. Therefore the 
bottom of the Schrodinger potential is at u ~ 1/Air. Similar arguments as above 
show that the mass gap, mass splitting, and decay constants are O (Air). As we 
pointed out above, for the scalars and pseudoscalars more careful or numerical anal¬ 
ysis is needed. The numerical result is that the masses are similarly O (Air), with the 
exception of pion masses in regime A, which obey the GOR relation. In regime B, 
pions also have masses O (Air). 

In conclusion, mass gaps, mass splittings, and decays constants are O (Air) at 
small quark mass, with the sole exception of the pions in regime A. 


B.2 Large quark mass 

Next we shall study the dependence on the (flavor nonsinglet) meson mass spectrum 
on the quark mass in the regime C (m q /Auv 1). Recall that in this limit some 
features of the bound state can be analyzed starting from field theory, because the low 
end of the spectrum becomes nonrelativistic. The expected mass gap is roughly equal 
to 2 m q , and the states can be studied by using the Schrodinger equation. For QCD 
states at large mass, one expects to have two main contributions in the Schrodinger 
potential. First, at very short distances one has Coulomb potential as perturbative 
gluon exchange dominates. Second, there is a nonperturbative confining potential 
which is expected to be linear in the distance and arise from a flux tube between 
the quarks. While it is difficult to derive such a linear potential from first principles, 
it is consistent with the observed quarkonium spectra, lattice simulations, and also 
found in holographic calculations. 
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In the limit mn q —> oo, the lowest states are therefore governed by the Coulomb 
potential, and one finds typical Hydrogen-like spectrum with negative binding ener¬ 
gies. For slightly higher states, the linear potential dominates. Assuming precisely 
linear confining potential ~ Aj R r, where Air is the scale of glueballs in QCD, it is 
straightforward to solve the Schrodinger equation. The masses of the states obey 
roughly the scaling law 

/a 4 \ V 3 

m n — 2m q ~ ^ ) n 2/3 . (B.ll) 

For charmonium and bottomonium the observed states are in the region where both 
the Coulomb exchange and nonperturbative effects are important (so that (B.ll), 
for example, is not a good approximation) and have positive binding energies. 

The behavior at large quark mass in V-QCD is somewhat dependent on the 
choice of the potentials. As pointed out above in (A.9), the scales m q and Air 
become separated as m q —> oo. For r -C 1 /m q one obtains the usual UV asymptotic 
solution, and for r 1/Air the background obeys the IR asymptotics, leading to 
the confining Schrodinger potential. In the middle, however, the behavior of the 
Schrodinger potential is nontrivial. To compute it, we first need to solve the tachyon 
from its EoM. 


B.2.1 The tachyon and the mass gap 

The tachyon background EoM may be written as [14,23] 


t" + 


3A' + A'^log(V f ,v) 


t' + e 


-2 A 
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4A' + A'^log(y fV ^) 


(r0 


/\3 


1 + e- 2A K(T') 2 d 


e~ 2A n 


dr 


log Vf — 0 . 


(B. 12) 


When l/m q <C r <C 1/Air, the tachyon has already grown large, and decoupled 
the quarks from the glue. The dilaton is still small, and its evolution is governed by 
the YM RG flow, as the quarks are decoupled. More precisely, in this regime 


logr + lo 0 4 + O ( log j rA[R j) • A y M l°g(rAi R ) + ® (log(rAm) 2 ) 

(B. 13) 

where £q = £(x = 0) is the UV AdS radius for YM, and b™ is the leading coefficient 
of the YM /3-function. The appearance of the IR scale Air in these expressions, 
which are the UV expansions for YM theory, may be surprising. In fact, it would 
perhaps be more appropriate to denote the scale of the expansions by a new quantity 
Ayy (which is not the same as Auv, the scale of the UV expansions for r -C l/m q ). 
But YM has only one energy scale, which is therefore the only scale in the model 
smaller than m q thanks to the decoupling of the quarks. That is, A/ T y rsj Air, and 
for simplicity we have already used this in (B.13). 
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Let us then insert the expansions in (B.12) in order to solve for the tachyon. 
We keep the leading behavior of all coefficients in 1/ log(rAjR). In particular, since 
A '/A! = O (log(rAiR) -2 ), the terms involving X' will be dropped. As we shall see 
below, this works at least for r close to the lower end of the scaling region. We also 
expect that the tachyon grows large, and therefore 1 <C e~ 2A k(t') 2 . The last term 
on the first line of (B.12) and the term on the second line dominate. Consequently 
the tachyon satisfies 

r\ 

4 r° r ' + ^ lo g^o, (R14) 

where we used the expansions (B.13), and k 0 = n(\ = 0). 

Let us study the asymptotics 

log Vf — —aoT Vp , (r —» oo) (B.15) 


where we require v p > 1 in order to ensure that the tachyon grows fast enough 26 . 
Then the “asymptotic” solution of (B.14) reads 27 

t — To (rrriqf , [y p = 2) (B.16) 

t ~ [£(2 - v p ) log(rm q ) + r 0 ]^ , (1 < v p < 2) (B.17) 


where 

* _ v P a pip 
4« 0 


(B.18) 


For v p > 2 there are no regular solutions. 

Let us then compute the Schrodinger potentials. First, the Schrodinger coordi¬ 
nate is given by 

^ ~ e~ A ^KT r ~ • (B.19) 

dr 

The behavior of 5 is dominated by its dependence on the tachyon in all sectors (flavor 
nonsinglet vectors, axials, scalars, and pseudoscalars): 


logS ~ 


1 


V p 


(B.20) 


Since logS increases fast enough with r, the first term in (B.3) can be approximated 
by 



d 

2 <i 2 , ,, 

d 

> log=(M) . 

+ ^ i °g=(„)= ! 

> log=w 


(B.21) 


26 The analysis can be extended to 0 < v p < 1 where some assumptions made above fail and a 
more careful analysis is needed. The tachyon still grows only logarithmically and the conclusions 
are similar to the case 1 < v p < 2. 

2 'Taking account the logarithmic flow of (B.13) would result in terms in (B.16) and (B.17) which 
have subleading (logarithmic) r dependence but are leading with respect to the constant To- For 
simplicity, we have omitted such corrections, as they do not affect our conclusions. The same applies 
to the formulas (B.29) and (B.30) below. 
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Combining these we obtain 


M ~ {r >f 

(m) \du 1 ’ 


d , , 

— iog; 

ar 


£ 2 0 ° 2 o v p r 2Vp ~ 2 

r 2 


(B.22) 


This result can be expressed in terms of the Schrodinger coordinate by using (B.16), 
(B.17) and (B.19), but the above form turns out to be more useful. The second term 
of (B.3) is negligible for the vectors and scalars, but important for the axials (and 
pseudoscalars), for which we find 


Ha — 


T 2 

wl r 2 ’ 


(B.23) 


where w 0 = w(X = 0). 

Let us first assume that the results (B.16) and (B.17) hold in the whole regime 
l/m q -C r <C 1/Air. We will later discuss when this is not the case. For 1 < v p < 2 we 
then obtain that up to logarithmic corrections, Vs ~ 1/r 2 . The potential decreases 
with r, and reaches its bottom at r ~ 1 /Air, where Vs ~ A 2 R . Therefore the 
meson mass gaps are characterized by Air, and only logarithmically enhanced with 
increasing m q . 

For v p = 2, the estimates (B.22) and (B.23) match up to the multiplicative 
coefficient. The Schrodinger potentials behave as 


Vs ~ 


( rm q ) 2 £ 
r 2 


( — < r -C 
\m q 



(B.24) 


If £ < 1, the result decreases with r. The bottom of the potential is reached at 
r ~ 1/Air, which leads to the mass gap 


^gap ~ A ir 



(B.25) 


If £ > 1 the result increases with r. The bottom of the potential is therefore at 
r ~ 1 /m q , and 

"igap ~ rn q . (B.26) 

In the marginal case £ = 1 the Schrodinger potential is flat and O ( m 2 ) in the whole 
regime. The mass gap is therefore also 


"igap ~ m q . (B.27) 

In conclusion, the mas gap has the power law expected for nonrelativistic states 
when £ > 1. However, as we shall see below, for £ > 1 the mass splitting is O (m 9 ), 
i.e., larger than that of nonrelativistic bound states. Therefore only the marginal 
case £ = 1 may potentially reproduce both realistic mass gap and splitting. 
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Let us then discuss logarithmic corrections to the critical case v p = 2. That is, 
we assume 

log Vf ~ — a 0 r 2 (logr)^ (B.28) 

as r —> oo. In this case, the tachyon solution is given by 


r ~ exp [ t 0 {rm q )Z] , 
r ~ exp j[£(l — Vi) \og{rm q ) + To] I “’k j , 


(vt = 1) 

(V£ < 1 ) 


(B.29) 

(B.30) 


and no regular solution is found when V( > 1. The leading term of the Schrodinger 
potential is now 

_ ^o a o v p T 2 (log r) 2vt 

S(w) “ 4 k 0 r 2 ' 1 j 

By similar analysis as above, the bottom of the Schrodinger potentials lies at r ~ 
1/Air and the mass gap is O (Air) when ve < 0, up to corrections which grow slower 
than any power of m q as m q —y oo. When 0 < ry < 1, the bottom of the the potential 
is at r ~ 1 /m q , and the mass gap is O ( m q ). 

Recall also that, as we pointed out above, in some cases the above analysis is not 
valid in the whole regime 1 jm q -C r 1/Air. This can happen if the terms involving 
A '(r) in (B.12) start to dominate as r approaches the value 1/Air. The growth of 
these terms requires that the factors in the square brackets depend on the tachyon, 
because A is known to obey the RG flow of YM in this regime, and remains small - 
actually r ~ 1/Air is exactly the region, where A finally reaches values O (1). Such a 
dependence on the tachyon may arise from the logarithmic derivative /A log Vf. The 
most natural Ansatz which leads to this is an exponential 


V f oc exp(—a(A)T^) (B.32) 

where the crucial point is that the factor a(A) depends on A. For the sake of generality, 
we however assume that 

^ log Vf - air dp (B.33) 

as r —> oo, where d v may differ from v p . 

From (B.12) we see that the terms which were neglected above become important 
when 

A'~\'-^logVf . (B.34) 

Inserting here (B.13) and (B.33), we obtain the condition 

[log(r Air )] 2 ~ r dp . (B.35) 

If d p > 0 and for any growing tachyon solution this condition is saturated within the 
range l/m q -C r 1 /Air, because at r ~ 1/Air the left hand side is O (1) and the 
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tachyon is already 1 (and for r ~ 1 /m q , the left hand side is sizeable whereas the 
tachyon is, by definition of m q , O (1)). 

For the sake of concreteness, let us consider the case v p = 2 so that the tachyon 
has the r-dependence of (B.16) when the additional terms oc A' are still small. 
Then (B.35) is saturated for 


m 


9 L 


log 


m n 


A 


IR 


2 


= r. 


(B.36) 


That is, the tachyon behaves as r * only up to r ~ r c , which is only enhanced with 
respect to r ~ l/m q by a logarithmic term. For r c r ■C 1/Air (he., for almost the 
whole interval l/m q -C r 1 /Air) we need to solve the tachyon EoM with different 
dominant terms. Assuming that the right hand side in (B.34) dominates over the 
left hand side, we obtain 


6 2AkX ' t ' §\ logVf ~^ logV f 

which becomes for the current case 

APh YM 

r'r^-i ~ |log(rA IR )] 2 . 

rai 

This equation is solved by 


q- 


4:d p £b 0 


,YM 


3ai 


[log(rA m )] 3 + To 


(B.37) 


(B.38) 


(B.39) 


so that the tachyon only increases logarithmically for r c r 1/Ajr. Inserting this 
expression in (B.22), and recalling that r c is close to 1 /m q , we see that the bottom of 
the Schrodinger potential is found at r ~ 1/Air, and the mass gap is again O (Air) 
up to logarithmic corrections. 

In conclusion, the phenomenologically interesting large mass gap (with power law 
dependence on m q ) is obtained only when v p = 2, 0 < vg < 1, and when // log Vf is 
suppressed (e.g., d p < 0). The last requirement means that for the critical exponential 
asymptotics log Vf ~ —ar 2 , the factor a cannot 28 depend on A. 

Interestingly we notice that potentials I, the construction of which was motivated 
by using completely independent arguments in [14,23], do produce a large mass gap. 
In this case we have 

V f ( A, r) = V f0 {\) exp(—a 0 r 2 ) (B.40) 


where ao is indeed independent of A. The value of «o was fixed in the UV (where the 
tachyon is small) to reproduce the UV dimension of the chiral condensate and quark 
mass so that 


_ v p aotl _ 3 /q 
Akq 4£ 2 


(B.41) 


28 More precisely, only dependence which is highly suppressed in the UV is allowed, for example 
terms ~ exp(—#/A). 
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The value of £ is therefore slightly below the critical value £ = 1 (because also 
(/^/d 2 < 1)- This is not surprising since holographic models are not expected to work 
perfectly in the UV region, where the coupling A is small, so fixing the parameters 
by UV arguments may lead to some tension in the model. Recall, for example, that 
the UV asymptotics of the pressure and the UV asymptotics of the correlators for 
the energy momentum tensor give slightly different numbers for the normalization of 
the glue action [79]. 

For potentials II, we used instead a(A) with nontrivial dependence on A, and 
consequently the meson mass gap will be small. 

Finally, let us point out that in all cases classified above, Ha of (B.5) contributes 
at leading order to the Schrodinger potentials near their bottoms. This is seen from 
the estimates (B.22) and (B.23) when the bottom is at r ~ 1 /Air, and by inserting 
the UV expansions in (B.5) if the bottom is at r ~ l/m q . Consequently, axial vectors 
and pseudoscalar mesons have larger mass gaps than vectors and scalars: the ratio 
of, say, the axial and the vector mass gap approaches a finite number which is larger 
than one in the limit of large quark mass 29 . 

B.2.2 Mass splittings 

Let us then discuss the mass splittings between the lowest meson states for the 
various cases described above. When the Schrodinger potential Vs(u) has a clear 
and regular minimum, the splittings can be estimated by computing the second 
derivative V/'(u 0 ) at the minimum u = u 0 . Most of the potentials discussed above 
fall into two categories, where the Schrodinger potential has a clear minimum either 
at r l/m q or at r rs./ 1/Air. 

When the minimum is at r ~ 1 /m q , we found that the mass gap was 0(m q ). 
From the definition (B.2) of the Schrodinger coordinate we see that u ~ r and the 
scale of the derivatives of all relevant functions is given by d/dr ~ m q . Therefore m q 
is the only scale that enters the analysis, and the splitting is also O ( m q ), i.e., much 
larger than expected for nonrelativistic bound states. Notice that this also includes 
some of the asymptotics with critical v p = 2, i.e., those with 0 < Vt < 1 , and those 
with Vi = 0 and £ > 1 . 

When the minimum is at r ~ 1/Air, the computation is more involved. We shall 
only discuss the case v p = 2 and V£ = 0, for which r ~ ( rm q Y and the mass gap is 
large, 



(B.42) 


Recall that the minimum is at r ~ 1/Air when 0 < £ < 1. At the minimum, all 
fields (r, A, and A) have logarithmic or power-law dependence on the coordinate r 

29 When 1 < v p < 2, the term Ha is actually leading by a logarithmic factor, so the mass gaps of 
the axials and pseudoscalars are logarithmically enhanced with respect to those of the vectors and 
the scalars. 
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so the scale of the r-derivatives is Air. The relevant quantities are, however, the 
w-derivatives of the potential, which can be estimated by using the chain rule and 
du/dr ~ r ~ (m q /Air)^. We find that 

/ , \ ( 2 -n)£ 

q'(«„)~A f R , V s n> ( u °) ~ A m" (B.43) 


where n > 2. The extent of the Schrodinger wave functions around the minimum 
is determined for the lowest fluctuation modes by Vg(u 0 ), and roughly given by 
A u = u — u 0 ~ 1/Air. Outside this region, the wave functions vanish very fast. 
The higher order derivatives in (B.43) vanish as m q —» oo, and therefore higher 
order terms in the Taylor expansion of the Schrodinger potential are suppressed (for 
u — Uq ~ 1/Air) and the potential takes the Harmonic oscillator form. In particular, 
the masses are given by 


m l = Vs(u 0 ) + 



+ o 



(B.44) 


where n = 0,1,2,.... From here we see that the mass splitting is suppressed as 
~ A m (Ai K /m q f at large m q . 

In the remaining case (v p = 2, Vi = 0, and £ = 1) the bottom of the Schrodinger 
potential, as obtained from the leading tachyon solution, is flat for l/m q <C r -C 
1/Air, suggesting a relatively small level splitting. There are, however, subleading 
logarithmic corrections to the solution due to the YM RG flow, which we have not 
computed. They will be important, and are expected to cause a minimum either at 
r ~ l/m q or at r ~ 1/Air. Then the mass splittings of the very lowest states will be 
as in one of the cases discussed above (with £ = 1 if the minimum is at r ~ 1 /Air), 
up to logarithmic corrections, unless the subleading corrections cancel miraculously. 


B.2.3 Decay constants 


Finally let us discuss the mass dependence of the (vector/axial) decay constants, 
given in (B.8). The Schrodinger wave function vanishes very rapidly whenever m\ < 
Vs(u), which limits the integral to the “classically allowed” region. 

First we notice that when the Schrodinger has its bottom at r ~ 1/Air., the decay 
constants of the lowest states are very small, because S in the above integral formula 
contains the factor y / V/ which is exponentially suppressed because the tachyon is 
also large near r ~ 1/Air. For the interesting case of v p = 2, V£ = 0, and 0 < £ < 1 
we find that 


f 2 

J n 

nJn c 


~ exp 



(B.45) 


One can also show that the pion decay constant has similar dependence on m q . 
Recall that Hy = 0 for the vectors and Ha is given in (B.5). The low-lying states 
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are therefore asymptotically decoupled. Sizeable decay constants are only found for 
highly excited states. More precisely, the suppression factor disappears when the 
classically allowed region extends to r ~ 1 /m q where the tachyon is no longer small. 
Since Vs ~ 1/u 2 in the UV, this requires m n ~ m q (whereas the lowest states had 
m n ~ AiR(m g /AiR,)^). Therefore the coupled states appear at the scale where mesons 
are expected in QCD. By using (B.8), it is possible to show that for such states the 
decay constants are f^/(N c Nj) rv-/ Air- 

When the Schrodinger potential has its bottom at r ~ 1 /m q , it is easy to see 
from the above expressions that 


f 2 

J n 

N c N f 




m„ 


(B.46) 


for the lowest states, because m q is the only scale which enters the formulas. 


C. Analysis of the free energy 

C.l Chiral condensate as the derivative of free energy 

The chiral condensate (including the sum over flavor) can be defined as 30 

d£ 1 0S on — sheii ,e 1 1 

dm q V 4 dm q 

where £ is the (zero temperature) energy density of QCD, and the subscript E denotes 
that Euclidean signature was used (so that the sign of the action is opposite with 
respect to the Minkowski signature, which was used in the main text). 

By computing the renormalized on-shell action using the identity (C.l) one can 
find the relation between the chiral condensate and the coefficient a in the vev term 
of the tachyon. This is, however, slightly complicated in V-QCD. First, since there 
is full backreaction: changing the value of the quark mass will affect the geometry, 
possibly leading to nontrivial contributions in the m q derivative of (C.l). Second, 
the counterterms needed to regularize the on-shell action depend on the quark mass, 
and may also contribute in the derivative. As it turns out, these issues can be fully 
solved in the limit of zero quark mass. 

Below we will first demonstrate how the free energy, and the backreaction in 
particular, can be analyzed by using the EoMs for the background. This will be 
compared to the direct computation of the regularized on-shell action done in Ap¬ 
pendix C.3. In the calculations below the fields shall be decomposed as 

A = A 0 + A t + A\ A = A 0 + A r + Ai t — To + ti + T q (C.2) 

30 Recall that the quark mass is defined only up to a proportionality constant in the holographic 
model, which we have set to unity for simplicity. The inverse of this constant would appear in the 
definition of the chiral condensate. 
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asymptotically in the UV (where r —y 0). Here the terms with subscript zero (one) 
are the source (vev) terms. The terms A r and A r and are sourced by the leading 
quadratic terms in the tachyon in the UV (oc Tq) and are therefore oc m q . The 
term r q is the leading nonlinear term of the tachyon. It is oc m q and computed in 
Appendix C.2. It is shown to be subleading with respect to the vev term n and 
therefore irrelevant in the calculations below. It is also argued that the same is 
true for similar terms due to quartic and higher order tachyon perturbation in the 
expansions of A and A. That is, the terms sourced by ToT q or Tq are subleading with 
respect to the vev terms A\ and Ai. The UV expansions of the various terms will be 
given in Appendix C.3. 


C.1.1 Chiral condensate at vanishing quark mass 


Let us start by computing the chiral condensate at zero quark mass. This is the 
simplest case because the square of the source term of the tachyon does not contribute 
in the on-shell action, and as we shall see, there is no issue with extracting value of 
a (the proportionality coefficient of the vev term) from the UV asymptotics. 

Consider a perturbation of a generic background (around m q = 0) which keeps 
Auv (but not the quark mass) fixed. The variation of the background solution in the 
UV is found by using the UV expansions from Appendix C.3: 


„4 r 


8A{r) ~ <5Ai(r) = 8Q — 

P 


MA(r) ~ MAr(r) = 
<5r(r 


1 + 0 


logrA 



9 

— B a a 5m a 

^4 

l + O i 

( 1 \ 

8 

32 q \ 

P 


y log r A J 


(C.3) 

(C.4) 


= 5m q r( 


8a r 3 


logrA) p 
— logrA) p 


1 + 0 


1 + 0 


log rA 

1 

log rA 


where 


B a = xw 0 n 0 e 5 


(C.5) 


(C.6) 


The chiral condensate may then be computed by studying the variation of the 
on-shell action as seen from (C.l). The computation is analogous to the holographic 
derivation of the first law of thermodynamics at finite temperature and chemical 
potential: the differential of the free energy equals the variation of the action, and 
expressing this in terms of the various UV and IR boundary terms gives the terms 
in the desired expression. In order to be as precise as possible, we will formulate the 
computation in terms of a conserved infinitesimal current, which can also be easily 
used in the computations later. It can be found as follows. 

The on-shell Lagrangian can be expressed as a total derivative [14]: 

Lon—shell oc — [—2A'e 3A ] . (C.7) 
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But any leading variation of a generic Lagrangian around its on-shell value is a total 
derivative as well, given formally by 31 


d 

dr 


sr a 


(C.8) 


where the sum is over all fields in the Lagrangian. Requiring equality of the generic 
expression and the variation of (C.7), one identifies the following infinitesimal con¬ 
served “current”: 

8e 3A \'5\ xe 3A VfKT r 5T 


J = 6e 3A A'5A - 6 e 3A 5A' - 


3 A 2 + e~ 2A K{r') 2 


(C.9) 


Indeed it is straightforward to check that J' = 0 by using the EoMs (and their 
variations). 

We will then require that the variation of the background is regular in the IR. By 
using the IR expansions of the background [14] in the expression (C.9) one sees that 
J vanishes in the IR (r —y oo) so it must vanish everywhere for regular variations: 


J = 0 . 


(C.10) 


In the UV the above expansions are inserted, leading to 

0 = lim J = —15 5Q - B a a5m q . (C.ll) 

r-> o 4 

The expression for the free energy density in terms of Q, cr, and m q can be extracted 
from the finite temperature computation in Appendix C.3 (see the expression for the 
free energy (C.63)). At small quark mass and at T = 0 (so that also (7 = 0) one 
obtains the expression 32 

£ = M 3 N 2 c (l5G + i. B a m g a) . (C.12) 

Inserting here (C.ll), 

5£ = -2M 3 N 2 c B (T a5m q . (C.13) 

Therefore the chiral condensate is given by 

(qq) = - 2M 3 N 2 B a a . (C.14) 

The coefficient in this equation can be fixed by using the asymptotics of the scalar- 
scalar correlator (see Appendix C in [23]), which leads to 

NtN 

{w) = —■ ( c - 15 ) 

31 In the case of the gravitational action there is a complication because R contains second deriva¬ 
tives. It is, however, well known that the second derivatives can be isolated in another total 
derivative term, related to the Gibbons-Hawking boundary term. 

32 We omitted contributions from the reference background used to regulate this expression, but 
as we shall argue below, this does not affect the result at m q = 0. 
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C.1.2 Chiral condensate at finite quark mass 

Let us then try to generalize the above computation to finite rn q . Recall that there 
is an issue in the definition of the vevs Q and cr. In principle, they are well defined 
as the coefficients of the vev terms A\ and T\, respectively. However such definitions 
are useless in practice, because the source terms A 0 and r 0 cannot be solved per- 
turbatively to high enough orders in order to separate them from the vev terms in 
the UV. In practice only differences of the vevs can be computed. Therefore we will 
define the values of Q and cr with respect to some reference solutions. Let us first 
discuss how we define a reference solution for all values of m q . 

First a (IR regular) solution with exactly zero tachyon (and therefore zero quark 
mass) is picked. This solution is chosen to have, by definition, Q = 0 (and trivially 
a = 0). Therefore A T and A\ in (C.2) are zero, and the solution for the scale factor 
A defines the source term Aq. Let us call this solution 1. Then solution 2 is chosen 
which has finite m q , defines the nonnormalizable term r 0 of the tachyon in the UV, 
and has a = 0 by definition. Further one requires that Q = 0 also for solution 2. 
Since Aq was already defined, this choice also defines A r . 

For 0 < x < x c , the solution 2 can be identified with the solution that was used 
to define the vanishing of a in the plots of Fig. 9. In principle any solution can be 
picked, but as we argued in Sec. 5, a choice which avoids fine tuning is to pick the 
solution 2 is the crossover between regimes A and B when 0 < x < x c and a solution 
with tiny quark mass in the conformal window. 

We have argued that “nonlinear” terms of the tachyon are suppressed with re¬ 
spect to the vev terms. Similarly, the terms sourced by quartic tachyons in the UV 
expansions for A and A are suppressed with respect to the vev terms. Therefore 
fixing the values of Q and cr for solutions 1 and 2 is enough to define the vevs for all 
backgrounds. It is also possible construct a reference background which has vanishing 
G and cr for arbitrary m q , as a (generally not IR regular) combination of the two IR 
regular solutions 1 and 2 constructed above by appropriately scaling the constructed 
A r and r 0 to have the desired quark mass. 

The above construction implies that the vevs can be written as 

G = G — Gi — rn q kA , a = a — m q k T (C.16) 

where the hatted quantities are the exact coefficients of the vev terms, and the 
quantities without hat are given by the above subtraction procedure. The coefficients 
ki and Gi are independent of m q . 

Let us then proceed with the calculation. Consider again a perturbation of a 
background with 5m q ^ 0, but now at generic value of m q . The UV expansion of the 


84 


perturbation has additional terms, in particular 


1 r 2 

5A = — — B a m q 5m q — (-logrA) 


—2 P 


1 + 0 


(: 


\ log r A 


4 


i + o 


logrA 


(C.17) 


where we chose to use the hatted vevs. The variation of b^X includes a term oc rn q dm q 
which turns out to only contribute at subleading orders in the UV, and the term 
(compare to (C.4)) 


45 - 

9 . 1 

4 


( i \1 

."if 49 - 

—B a ( o5m q + m q 5a) 

P 

i + o 

y log r A J 


(C.18) 


It is again useful to study the current (C.9). As above, J = 0 because it vanishes 
in the IR limit. In the UV one obtains 


0 = lim J 

r —>-0 


1 9 

— 15 5Q — -B a m q 5a — ~^B a a 5m q . 


(C.19) 


In particular, we expect that no terms oc m q 5m q arise in the UV because a and Q 
are exactly the coefficients of the vev terms. This identity may be rewritten as 


5 


15£ + -B a m 



—2B a a 5m q , 


(C.20) 


where the expression in the square brackets has similar structure to the regularized 
vacuum energy given above in (C.12). The natural expectation is that indeed 


S 


m 3 n 2 c 


15 Q + -B a m 



(C.21) 


but in principle the expression could also contain additional terms oc m 2 (or con¬ 
stants) which would cancel in the regularization. 

We see that 


5S = -2 M 3 N^B a a 5m q = -2 M 3 N*B a (a + k T m q ) 5m q , (C.22) 

where extra terms oc m 2 in the definition of S would effectively change the value of 
k T which remains unknown in any case 33 . Therefore the chiral condensate is given 
by 

(qq) = -2M 3 N*B a a = -2 M 3 N 2 c B a {a + k T m q ) . (C.23) 

33 A possible handle to control the value of this coefficient is to fix somehow the subleading linear 
corrections to the GOR relation discussed in Appendix D. 
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C.1.3 Free energy on the Efimov spiral 


Let us then discuss how the free energy is computed for the configurations of Fig. 9, 
which are found when 0 < x < x c . Let us assume for simplicity that we are able to 
define the vevs such that the constant k T of (C.23) vanishes. In this case the chiral 
condensate is simply given by the derivative of the free energy density with respect 
to the quark mass: 

d£ 

{qq) = ^ = -‘2M :i N^B <T cr = -N f N cC(J a . (C.24) 

Inserting the asymptotic result of (5.14) in (C.24) and integrating, one can readily 
find the free energy for asymptotically small quark mass: 

JV/JVcAJv {S ~ £o) 


where Sq is the free energy of the solution having m q = 0 = cr, and u = logAuv/Ai R . 
In general, free energy differences are given by the area between the spiral and the 
horizontal axis (see, for example, Fig. 12). 

Interestingly, the result for both zero m q and zero o simplifies to 


r K 2 

L cM v IR 


—4 U 


2 KmKfr sin(0 m ^cr)^ 


(C.25) 
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-sin(0i R -0 CT -fcn)sin(0 IR -(/) m -/cn) + - sin (^> m ~<f> a ) 
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vc„K 2 r 


—4u 


2 A UV A^jy 


QK-mK-a Sm(0 m (^ ct ) 


(C.26) 


NfN c Ai 


<£-£o) = 


vc„K 2 { 


uv 


sill [ O rT} 0o-) 


e _4 “ < 0 , 


(C.27) 


where we used the handedness of the spiral in (5.16): sin(0 m — (j) a ) > 0. Therefore 
these solutions are dominant over the solution with zero tachyon. 


C.1.4 Free energy with multi-trace deformations 

Let us then study how the free energy is computed in the presence of multi-trace 
deformations. Recall that the UV boundary conditions are 

Timax 

am = m q + ^2 Sncr^r 1 ( C - 28 ) 

n=2 

Pm = a (C.29) 

and that the identity (C.19) is interpreted as 

1 9 

0 15 5Q -B a a m 5(3 m ~^B a f3 m 5cx m 
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By using the UV boundary conditions, this identity can be rearranged to read 


^max -i 

= - c a a5m q - V —c™cr n 5g n . 

n =2 

(C.31) 

Therefore the quantity in the square brackets is identified as the new free energy (over 
NfN c ). It matches with the free energy obtained by renormalizing the action (C.12) 
np to the last term involving the couplings g n . This extra term is proportional to 


15M 3 „ i 

Q T g c a cx m <7 T / j 

n =2 


X 


n 


n 


g n c n a a n 


w-p m 


8W 

8/3 m (x) 


(C.32) 


where W is as in (6.2). Therefore the term agrees with that found in [61]. 
The condensates are given by 


(O) 

(O n ) 


1 d£ 
NfN c dm q 
(— l) n n d£ 
NfN c dm 


= -c a a 



CgO~ ) 


(C.33) 

(C.34) 


where the normalization factors were read from (6.1). Therefore we find agreement 
with the large N factorization of the expectation values. 

The free energy on the Efimov spiral asymptotically close to the origin can be 
computed as above. That is, by inserting the spiral equations (6.8) to (C.33) and 
integrating, we obtain 


NfNcA^ 


{S-So) = 


c a m q a 
2 At TV A?' 


UV ^UV 
VCgKf R 


1 t max ' \ 

E 


n =2 


2 n 


—Au 


8K m KgSm((j) m -(j ) a ) 


(C.35) 


Notice that the term involving g n in (C.35) vanishes for n = 2, which is consistent 
with nonzero gq amounting to a redefinition of the parameters of the Efimov spiral. 
The higher order terms are suppressed as u grows because a ~ exp(— 2u). Therefore 
we conclude as above that the solutions with m q = 0 but a ^ 0 dominate over the 
solution with cr = 0, and that the solution with the smallest u (for u 1 so that the 
parametrization of the spiral is accurate) has the lowest free energy. 


C.2 Contributions oc m z q in the on-shell action 

The on-shell action may involve terms oc m q which arise due to the nonlinear nature 
of the tachyon EoM. We will now check how these terms behave. We will need the 
UV asymptotics of the source and the vev terms of the tachyon which are given in 
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equations (C.48) and in (C.54) in Appendix C.3. Using results from [14], the tachyon 
EoM can be written in the UV as 


r 


U + 1(3 


2p 


log(rA) 




(C.36) 


where k 0 = k( 0). All contributions suppressed by l/log(rA ) 2 or more are hidden 
in the coefficients of the terms linear in the tachyon, as well as all contributions 
suppressed by 1 / log(rA) in the nonlinear terms (including the complete terms oc t 2 t' 
and oc r 2 (r') 3 ). 

The “nonlinear” tachyon solution is then found by writing 


r(r) = r 0 (r) + r g (r) , 


(C.37) 


where the UV expansion of To is given in (C.48) and r q is the term which needs to 
be solved. Using the fact the To solves the linear equation, one finds that 


r'g - ( 3 


2 P 


log(rA) 


n~ 

1 q — 


4^or , 3 3/to , 


~ Kom?l (— log(rA)) 3p 


(C.38) 


One finds that (dropping the terms corresponding to the source and vev terms, which 
arise as solutions to the homogeneous equation) 


Tq = k 0 m 3 

i 2(4p — 1) 


r 3 (— log(rA)) 


— 3p+l 


1 + 0 


log(rA) 


(C.39) 


which agrees with the result of [46] when p = 0. 

Notice that the nonlinear term is subleading to the vev term in (C.54) provided 

that 

P > i , (C.40) 

which is satisfied for QCD in the Veneziano limit as 


= To = 9 > _9_ 1 

P ~ bo ~ 22 - 4x ~ 22 4 ’ 


(C.41) 


where b o and 70 are the leading coefficients of the beta function and the anomalous 
dimension of the quark mass, respectively. Therefore one can conclude that only 
linear terms of the tachyon EoM are relevant in the computation of the on-shell 
action 34 . The conclusion is different from that of [46] where the quark mass did 
not run. Notice, however, that these terms are only logarithmically suppressed and 
would be important for large values of m q if a finite value of the UV cutoff was used. 

34 Notice that even though the nonlinearities in the tachyon asymptotics do not appear directly in 
the computation of the on-shcll action, nonlinear terms of the EoM are important in general because 
they affect the values of the vevs such as the chiral condensate through IR boundary conditions. 












C.3 Regularized on-shell action 

The free energy is given by the on-shell value of the action. It is a UV divergent 
quantity and needs to be renormalized 35 . We work here at finite temperature, and 
subtract the UV divergences by considering the difference with respect to a reference 
(zero temperature) background with the same values for the sources. We could also 
consider the difference between two different zero temperature backgrounds. 

It is not difficult to show that the Lagrangian is a total derivative even in the 
presence of the tachyon (as already shown at zero temperature in [14]). After a 
straightforward computation, and taking into account the Gibbons-Hawking term, 
one finds that 

Son—shell = M 3 N^f3V 3 (6A'f + f)e 3A \ r=e + counterterms . (C.42) 

When the quark mass is finite, the tachyon contributes at O (r 2 ) (terms oc m q ) 
and at O (r 4 ) (terms oc m q a ) in the UV expansions. These contributions need to 
be taken into account in the holographic renormalization procedure. The calculation 
of [75] at Nf = 0 used a reference (thermal gas) solution to subtract the divergences. 
We will use the same technique here. Denoting the fields and other variables of the 
reference solution by tildes, the following conditions need to be fulfilled at the UV 
cutoff: 


/?e A ^ = /3e A(e) v / 7M , V 3 e sA ® = U 3 e 3A(e) , (C.43) 

A(e) = A(e) , r(e) = r(e) , (C.44) 

where the possibility that the cutoffs of the two solutions are different was included, 
e / e, as in [75]. In the absence of the tachyon this was convenient since the UV 
scales A = Auv of the two solutions could be chosen to be the same, A = A. The 
fact that the present system has two scalars suggest that it is better to choose here 
e = e and satisfy the last two conditions by varying the sources A and m q . However, 
it turns out to be convenient to still require that A = A and vary the cutoff (and m q ) 
instead. This is so because the tachyon contributions are independent on whether one 
chooses to vary A or e: The variation if suppressed by O (e 4 ), i.e., e/e = 1 + O (e 4 ), 
so that the effect on the tachyon contributions will be down by O (e 6 ) and therefore 
negligible. Keeping A fixed one can maintain very close contact to the computation 
of [75], 

Turning on the quark mass will modify the UV expansions of A and A. at this 
point it is useful to write down the complete UV expansion including all relevant 
terms. One can decompose, as r —>■ 0, 

A = Aq + A t + Ai , A = Aq + A t + Ai, / = 1 + fi , t = tq + ti. (C.45) 


35 For detailed analysis of the holographic renormalization of dilaton gravity, see [60,80] 



The source terms have the expansions [14] 


A 0 = - log r + log £ + 

1 O 


bo^o 

To 


log rA 


9 log r A 
1 


+ o 


(log rA ) 5 


(log rA)' 


— = m q r(— log r A) 


-p 


1 + 0 


logrA 


(C.46) 

(C.47) 

(C.48) 


and p = 70/&0 with 70 being the leading coefficient of the anomalous dimension of the 
quark mass. The O (r 2 ) terms A r and A r were also included here which are sourced 
by the O (r 2 ) tachyon perturbation 36 . As it turns out, A r is suppressed by logarithms 
of log rA with respect to A r so that it will not enter the calculation and its expansion 
will not be needed. A r is given by 


A = ~^ B ^ m +p (-logrA) 2p 


1 + 0 


log r A 


where Wo = V/o(0), n 0 = k( 0 ), and as above 

B a = xW 0 k 0 £ 5 . 

The vev terms have the expansions (compare to [75]) 


Mr) = 


1 + 0 


log rA 


45 9 \ r 

feoAi(r) = ( ~—G - ) p 


4 r 


1 + 0 


logrA 


C r 


/l(r) = ~4 p 
n (r) 

—■— = err 


1 + 0 

logrA) p 


logrA 
1 + 0 


. 

log r A J 


(C.49) 

(C.50) 

(C.51) 

(C.52) 

(C.53) 

(C.54) 


Let us then go on with the renormalization procedure. In the expressions below 
only the differences of the vevs Q and a with respect to the values of the reference 
solutions will appear. Therefore, without loss of generality, one can set the vevs of 
the reference solution to zero. As motivated above, one can choose A = A. Then the 
relation between e and e is fixed 3 ' by the first condition in (C.44): 


e 

e 


1 + 


45 


9 


- ^Am q a ) ^(-logeA ) 2 


1 + 0 



(C.55) 


36 There are also O (r 4 ) terms which arise due to the perturbation from the quartic terms in the 
tachyon. As we argued above, these contributions are subleading with respect to those arising from 
the vev terms. 

37 It can be checked that the variation of m q , which would enter through the term A T is subleading. 
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whereas the second condition sets 


m q = m q + cre 2 (— log eA) 


2 P 


1 + 0 


log eA 


(C.56) 


It is then straightforward to calculate the renormalized pressure. Inserting the 
expression of the (unrenormalized) on-shell action from Eq. (C.42) one finds that 


-fiVk p = lim 

e ->0 . 


S(e) - S(e) 


(C.57) 


= M 3 IV 2 lim 


>Va (6 A'(e)f(e) + f(e))e 3A ^ - 6 ^V 3 A\e)e 3A ^] . (C.58) 

After eliminating /3 and V 3 by using the conditions (C.43), the expression for the 
nressure reads 

vTW - 6A'(e)/(e) - /'(e)' , 


pressure reads 

p = M 3 Nj. lim e 3A(e) 

C e—>0 

The relation (C.55) implies that 

A'(i)e- A ® = A'(e)e~ A ^ + (+ \ 

\2 o 


B a m q a | — 


1 + 0 


. log eA / 


(C.59) 


(C.60) 


where the expansion (C.46) was used. Further, notice that the variation of m q enters 
through 

1 A r / 1 \ 1 

(C.61) 


~ 1 e 4 

A T (e) = A r (e) + —B a m q a — 


1+ °{++) 


so that 


A(e) = i(e) + 


r 1 

e 4 

1 \ 1 

y + — B a m q (j 

lo 

P 

[ i + U (logeA)J 


(C.62) 


similar result can De iouna ror tne derivative of A. Inserting these relations and 
the expansion of / in (C.59) one obtains the final result 38 


p = M 3 N 1 (-C-\$g--B < 


■m„a 


D. Gell-Mann-Oakes-Renner relation 

In this Appendix we check the GOR relation explicitly. The starting point is the 
normalization integral (E.13) for the pion wave function. When m q is small, it is 
dominated 39 at small u ~ r. Also the wave function ipp is constant in the UV up to 

38 Recall that the vevs here should be interpreted as their differences with respect to the values 
of the reference solutions, because we set the vevs of the reference solution to zero. 

39 As can be verified numerically, this holds for the pion state but not for higher pseudoscalar 
states. 
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corrections O (r 2 ) (see Appendix G in [23]). Therefore integral one needs to compute 
is (the UV contribution to) 


I = dr - 


VfK,e 3A r 2 


(D.l) 


Let us first try to compute the integral by using the UV expansions of the various 
fields, given in Appendix C.3. We obtain 


I = 


WW 3 

1 


dr 


T(r)‘ 


1 + 0 


logrA 


(D-2) 


dr r 3 


W 0 k 0 £ 5 J 0 


m q r (—log rA) 


-p 


1 + 0 


log rA 


+ ot 3 (—logrA) p 


1 + 0 


log rA 


-n 2 


where A = Auv and k 0 = k( 0). As m q —» 0 this integral is dominated by the regime 
with r ~ > Jm q /a . Substituting here v ~ r(— logrA) p one obtains 


I = 


1 


W 0 K 0 O 


foo V 

dv — 


1 + 0 


i 


log vA 


'0 


( m q v + crv 3 y 


2W 0 n 0 £ 5 m q a 


1 + 0 


log(m q A 2 /a ) 


(D.3) 

Notice that in the walking regime corrections are suppressed only when m q is a small 
perturbation, i.e., (A.10) holds. Also, the result (D.3) is not valid in the conformal 
window for any m q . 

From (D.3) we see that the UV RG flow leads to correction terms which are 
suppressed only by logarithms of m q . It is however possible to show that such cor¬ 
rections vanish to all orders. First we notice that the tachyon EoM (B.12) can be 
written as 

r" + log ( e 3A V f0 K ) r + = 0 (D.4) 


K 


up to corrections suppressed by r 2 in the UV. Let us denote by T a the solution having 
m q = 0 and a — 1, and let r be a generic solution with m q ^ 0. The Wronskian 


W = tt' - 


satisfies 


with the solution 


W + ^_ log (e 3A V f0 O) W = 0 

C w 


W = TtL- TrrT = 


(D.5) 

(D.6) 

(D.7) 


e 3A V f0 K ' 

By inserting the UV expansions for all fields we see that the constant is given by 

Cw = 2WoK 0 £ 5 m q . (D.8) 
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Notice that (D.7) may be written as 

d Try 


a 


w 


dr t e 3A Vf 0 KT 2 

By integrating this identity, the UV contribution to the integral (D.l) becomes 


(D.9) 




rcut _ 1 r a (r CVLt ) 

o C w T{r cu t) 


lo e 3A V f0 KT 2 C w r 0 C w t{t cu t) ’ (D-10) 

Choosing the cutoff such that r cut \Jm q /a (but so that it is not too large so that 

our approximation are still valid) we obtain 


I = 


1 + 0 


m n 


ar 


cut 


(D.ll) 


/ 0 e 3A Vf 0 KT 2 C w a 

Our approximations will break down at r cut ~ 1/Auv, where the r-dependence of 
the tachyon changes qualitatively. In order to complete the calculation, one should 
estimate the contributions to the normalization integral (E.13) for r > 1/Auv- But 
in this regime the dependence on the tachyon is regular and one can just do a Taylor 
expansion at m q = 0. The cutoff dependence should cancel against the UV con¬ 
tribution given in (D.ll). As the dust clears, we can effectively set r cut ~ 1/Auv, 
in (D.ll), obtaining 


/ = 


Cw o’ 


1 + 0 


m 9 A UV 

cr 


The normalization condition (E.13) becomes for small quark mass 

*,.( 0 ) 


1 _ V’P,7t(0)^ 


1+0 


m q A UV 

a 


2£ 5 W 0 K 0 m q a 


1+0 


m Q A uy 

O’ 


(D.12) 


(D.13) 


where we restricted to the lowest pseudoscalar mode, the pion. Solving for ipp )7T ( 0) 
and inserting in (E.21) one obtains the GOR relation 


//m^. ~ 2M 3 NfN c £ 5 W 0 K 0 m q a = 2 M 3 N 2 c B a m q o ~ — m q (qq) 


(D.14) 


with corrections suppressed by m q K\ ]V /a. It can be checked that the proportionality 
factor is correct for our definitions of f n and (qq). 


E. Fluctuation equations, f w , and the S-parameter 

The radial wave functions for the flavor nonsinglet fluctuations satisfy the following 
equations [23]: 


d u {Cvd u 'i(jv) — Cyq' J i\)v 

(E.l) 

duiCydu^A) = Cy(H A + q 2 )i/j A 

(E.2) 

d u (Cyd u ^ L ) = CyH A (lpL ~ tpp) 

(E.3) 

H A d u il) P = -q 2 d u ^> L 

(E.4) 
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where 


— = v/l + e~ 2A n(r') 2 = G , 
dr 


H a = 


4 T 2 ne 2A 


w “ 


Gy = VfW 2 e A . 

In addition it is convenient to define the pseudoscalar wave function 

" 0 P = —Cydu^L 

which satisfies the single equation 


(E.5) 
(E.6) 

(E.7) 

(E.8) 

C v H a V f T 2 ne 3A ' ^ 

At nonzero quark mass the pion decay constant may be defined in terms of the 
pole of n L in (7.16) at q 2 = — m As the computation is similar also for higher 
modes, a generic pseudoscalar fluctuation can be considered. To compute the decay 
constants one needs to study the fluctuations for small values of 


where 


d u {Cpd u ^p) = C P (H A + q 2 )^ P 

c P = 1 


6q 2 = q 2 + m 2 n (E.10) 

where m 2 is the mass of the fluctuation mode. The wave functions are written as 


pJL = k L + Sif-ipL + O ((Sq 2 ) 2 ) 
ipp = k L ipp.n + Sq 2 p> P + O ({Sq 2 ) 2 ) 


(E. 11) 
(E. 12) 


where /cp is a normalization constant which will be fixed below. 

When 5q 2 = 0 the wave functions are normalizable. One can therefore choose 
the fields ^/.n to satisfy the usual normalization condition 


1 = 


duCp^p^ = 


du iCv '»-'- 22 


H, 


' ( 9u^pL,r 


(E. 13) 


with V'p.n defined as in (E.7) for the normalizable mode. At finite but small 6q 2 we 
impose the standard normalization condition in the UV: 


1 = V’i(o) - k L Pq 2 'ii L ( 0) ; 0 = ipp( 0) oc MO) , 


(E. 14) 


where we used the fact that the normalizable wave functions vanish in the UV. 
Expanding the fluctuation equations at small 5q 2 gives 


d u {Cvd u 'ijjL ) — CyH A {^>L — -0p) 
H A d u ^) P = mldu'tpL - d u i> Lt1l . 


(E.15) 

(E.16) 
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By using these equations and the fluctuation equations for ij} I n one finds the identity 


d u 


m n C V ( "0 L dui>L,n - -0 L,ndui>L ) + H A C V ( 0p9 u 0p, n - 0p, n 9 u 0 


— HaCv'iI’P^H’P^ ~ 1pL,n) 


(E.17) 


Integrating this over u, using the boundary conditions from above, and further using 
the fluctuation equations to simplify the result one obtains 

1 f°° 

Cvi>Ldu^L,n\ u=0 = -2/ du 0 P,n< 9 « 0 P,n . (E.18) 

m n Jo 

Inserting (E.14), integrating partially, and using again fluctuation equations this 
relation simplifies to 


C'v9u' 1 Pl,ti\ u= q 


k L Sq 


2 poo 


duCp^ 2 Pn = 


k L dq 2 


(E.19) 


By using this result and the definitions from (E.ll), lip from (7.16) can be written 


as 


^ M 3 N f N c a , , M 3 N f N c , a , 

U d q ) = 4g2 Cvd u M u =o ~ 4q 2 CvkrA^u 


I u =0 


M 3 N f N c 1 n , 

-{Cydu^L^n) 


u =0 


q 2 5q 2 

for small 5q 2 . The decay constants are given in terms of the residue at 5q 2 = 0 


(E.20) 


fi< = MW / JV c (Cya„V’i,»)1„,„ = M J jV / jV t # f .,„) 2 |„ =0 


From (E.20) we see that for q 2 ~ m 2 , as m q —y 0, 


n L (q 2 ) ^ 


2 \ _ , f>l 


f 2 

J 7r 


f 2 

J 7T 


g 2 (g 2 + m 2 ) g 2 g 2 + m 2 


(E.21) 


(E.22) 


so that So (defined as the zero momentum value of g 2 IIp(g 2 )) approaches f 2 . The 
same factor for the axial-axial correlator does not have a pion node, so that simply 


IW) - 4 

qz 


(E.23) 


in the same scaling limit. Consequently, the axial-axial correlator has the correct 
structure at small momentum: 


(qV - «v) m(« 2 ) + ?vin( 9 2 ) = ir 


qfiqu 


(E.24) 


g 2 + m 2 

That is, the longitudinal term arises from the coupling of the pion to the axial current. 
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Let us then derive a few integral representations which are useful when computing 
observables numerically. First we analyze 

D(q 2 ) = q 2 Ii A {q 2 ) - q 2 H A {q 2 ) = ^ M 3 N f N c C v [d u ^ v (u,q 2 ) - d u ^ v (u, q 2 )\ u=e , 

(E.25) 

where we inserted (7.16), and used the fact that r ~ u near the boundary. Recalling 
the normalization ?/y (u = e) = 1 = ij)v{u = e) we obtain 


C v [dui>v ~ < 9 rVvL =e = J dud u [VvCVduVt 4 - VuCVduVv] 

poo 

= du 'il> v H A C v , ijj A , 


(E.26) 


where fluctuation equations (E.l) were used at the last step. Taking e —^ 0, it follows 
that 

1 


D(q 2 ) = -M 3 N f N c 


du ^vH a Cv^a ■ 


(E.27) 


A rather similar formula can be derived for the S-parameter. Let us denote 


d l 

w ivlA 


q 2 =0 


— Vv/a • 


Then the S-parameter (7.18) can be written as 

S = t rM 3 N f N c C v d u (i> A - V 

Derivating the fluctuation equations with respect to q 2 results in 

d u {C v d u ipv) = CWv 
d u {C v d u ^ A ) = C V ^ A + C v H A ip A 


(E.28) 


(E.29) 


(E.30) 


where it is understood that ^v,A are evaluated at q 2 = 0. The combination appearing 
in the S-parameter can then be written as 


Cyd u 


ipA-ipv — J du d u vC v d u ip v - i> v C v d u i) V 

- i) A C v dui) A + i> A C v d u i> A (E.31) 

POO 

= / du C v (ip a ~ V’y) , (E.32) 


where we dropped higher order terms in e in the first step (noticing that 'ipv,A vanish 
fast in the UV) and used the fluctuation equations (E.l) as well as equations (E.30) 
in the second step. Taking e —» 0, we obtain the final result 


S = nM 3 N f N c j du C v - ^ 2 V \ q2=0 


(E.33) 


Here one could also insert that i/jy = 1 when q 2 = 0. 
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F. Free field (one-loop) computation of the S-parameter 

The free field result for the vector-vector and axial-axial correlators is given by 


d A x e~ iqx (0| T{j a ^ v / A \x)J b l , {y/A \ 0)} |0) 


N c S ab r d A k Tr [ 7 ^( 75 )(ft + m)Y ( 75 )(ft - j + m)\ 

2 J (27t) 4 ( k 2 + m 2 — ie) [(k — q ) 2 + m? — ie] ’ 

where the 75 ’s are only present in the axial-axial correlator. The loop is divergent 
but the contribution to the S-parameter will be finite. 

Doing the integral with dimensional regularization one obtains 


r' 

■ q 1 (4U 2 X 


dxx( 1 — x) 


x-log(m 2 + x(l — x)q 2 ) — 7 + log 47 T + O (e) 


n _ 2AW f 1 

n - i{q >-wl 


dx x(l — x) 4 -— 


x-log(m 2 + a;(l — x)q 2 ) — 7 + log 47 T + O (e) 


(4tt) z q 2 


x dx -log(m 2 + x(l — x)q ) — 7 + log 47 r + O (e) . (F.4) 

./n £ 


Consequently 


D(g 2 ) - D(0) = - 


2 N f N c rrr 


da; log 1 + a;(l — a;)—- , 


where the value at q 2 = 0 was subtracted in order to remove the logarithmic diver¬ 


gence. 

For the S-parameter one obtains the well-known result 


5 = -47t£>'(0) = 


N f N c f 1 


dxx( 1 — x) = 


N f N c 
12t r ’ 


whereas the corresponding finite difference reads 


D(q 2 )-D( 0) _ N f N c m 2 


27r q 2 


— / da; log l + a;(l —a;)—- 


NfN c m 


—— * /1 T 4—— arctanh- 


7r 


'1+4^ 

<7 


1 (F. 8 ) 
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in agreement with [ 68 ]. The series as q 2 —> 0 and as m 2 —> 0 are given by 


—47T 


D(q 2 )-D( 0) = N f N c 
q 2 127r 

_ N f N c 


7T 


1 q 2 1 f 

1 - - - 1 -—-)- 

10 m 2 70 m 4 

2 2 2 4 

rn m rn m 

q 1 2 q 1 q 1 q 4 


(F.9) 


4 2 

m m 

— log — + 
q 4 q 2 


G. Details on numerics 

The numerical results in this paper were computed for sets of potentials termed 
“potentials 1” and “potentials II”. They are exactly the sets given in [23], but we 
repeat their definitions here for completeness: 

• Both Potentials I & II. 

Vg{ A) = Po 
Yfo(A) = IPo 

• Potentials I. 

a(A) = a 0 , k(A) = ----. (G.2 ) 

(1 + ^a A ) 413 

• Potentials II. 

1 + aiA + i 

«(A) = a o - - -n§- , «(A) = --. (G.3) 

i + i) 0 + ^r 


1 + PiA + P 2 A : 


yi+ioga+A) 


1 + i; 


2/3 


[l + W 1 X + W 2 X 2 ' . 


(G.l) 


Here the most of the coefficients are fixed by matching to perturbative QCD. 
Exceptions include the normalizations factors Vo, which Exes the UV AdS radius, 
and Wo, which remains as a free parameter. We also set £(x = 0) = 1, and choose the 
parameter Ao, which only affects the higher order coefficients of the UV expansions, 
such that the higher order coefficients have approximately the same relative size as 
with standard scheme choices in perturbative QCD. Explicitly, the coefficients satisfy 


U 0 = 12 , 

Wi = 


V,= 
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27tt 2 ’ 
24+ (11 -2x)W 0 
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(G.4) 



Two qualitatively different choices for W$ are possible: either constant Wo, which 
satisfies 

0 < Wq < 24/11 , (G.5) 


or Wq fixed such that the pressure agrees with the Stefan-Boltzmann (SB) result at 
high temperatures [18] (without the need to include x dependence in the normal¬ 
ization of the action). The latter option is given explicitly (when £(x = 0) = 1) 
by 


12 

Wo — — 
x 


(1 + H 2 /3J 


(Stefan-Boltzmann) , 


(G.6) 


so that the AdS radius is 


£(x) 




(G.7) 


In this article we have always chosen Wo = 3/11 for potentials I, and the SB normal¬ 
ized Wo for potentials II. For the coupling of the gauge fields w(X) which is required 
for the computation of the vector correlators, we used w(X) = n(X) for potentials I 
and w(X) = 1 for potentials II. 

The numerical result were mostly computed as detailed in [14,18,23]: coupled 
ordinary differential equations were solved by shooting from the IR, either starting 
from near the IR singularity (at zero temperature) or near the horizon (at finite 
temperature). In order to obtain accurate and reliable results, some tricks had to be 
used in various cases, in particular near the BZ point x = Xbz and at large quark 
masses. 

In general, the difficulties of the numerics in V-QCD arise from two sources: the 
IR singularity and the logarithmic corrections to the UV asymptotics. The latter 
easily lead to sizeable numerical errors when one tries to extract the values of the 
sources or vevs near the boundary. To improve the accuracy of the numerical analysis, 
we did the following: 


• The scale factor A was used as the coordinate instead of r. This makes it easier 
to analyze the UV asymptotics, because the mapping from r to A is logarithmic 
and expands the details near the boundary. 

• Near the IR singularity (at zero temperature), where possible, EoMs were 
implemented such that they do not contain the extremely small factors oc 
exp(—or 2 ). In many cases, this could be done by writing all expression in 
terms of Vi 0 g(A,r) = log V/(A, r) = —or 2 + logV/ 0 (A) rather than in terms of 

Vf- 

• Near the BZ region and at very large m q , the zero temperature backgrounds 
were constructed by shooting from the IR toward the UV in four steps. Very 
close to the singularity, where the tachyon is large and completely decouples the 
flavors, YM solution was used for the geometry and the dilaton, whereas the 
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tachyon was solved by using a simplified EoM with only the dominant terms 
on top of the YM background. In the next step, the complete tachyon EoM 
together with the decoupled YM flow for the geometry was used. In the third 
step, all fields were solved from the full coupled EoMs. In the final step, the 
tachyon was again solved separately since it is decoupled from the other fields 
near the boundary. Actually f = e A r was used as the field, because it decreases 
much slower than r, and the flow can then be tracked closer to the boundary. 

The first two steps were necessary because the tachyon EoM becomes stiff in 
the IR, in particular for large m q , so that the tachyon takes larger values in 
the IR than otherwise. The last step was useful in particular at large x, where 
the RG flow of the fields is slower, and one needs to solve them closer to the 
boundary in order to control the leading logarithmic behavior of the various 
terms. 

• For the finite temperature backgrounds in the BZ region and at large rn q , the 
background was solved in two steps, which were analogous to the two last steps 
of the zero temperature construction. That is, the full system was used in the 
IR, and the tachyon was treated as a decoupled field near the UV. 

• The S-parameter was computed by using the integral formula (7.21) rather 
than UV asymptotics of the fluctuation wave functions. 

• The nonsinglet meson wave functions were computed by rewriting the fluctua¬ 
tion equations into a system of two coupled first order equations (rather than 
a second order equation). Again the equations were written in a form which 
does not contain explicitly the tiny factors exp(—ar 2 ). 
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